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PREFACE 


Hie  purpose  of  this  project  was  to  generate  and  detect 
features  in  connected  speech.  It  is  a  part  of  the  larger 
speech  recognition  problem. 

Phonetic  units,  which  we  called  observations  and 
phonets.  were  generated  from  connected  speech.  The  time 
file  was  converted  to  feature  space  using  the  Fourier 
Transform.  Phonet  occurrences  were  detected  using  Minkowski 
One  and  Two  distance  measures.  Phonet  matches  were  ranked 
in  order  from  minimum  to  maximum  distance  between  each 
input  phonet  and  each  phonet  in  the  template  file. 

An  algorithm  was  developed  to  partition  feature  space 
into  classes  of  phonetic  units.  The  central  tendency  of 
each  class  was  considered  to  be  a  template  against  which 
observations  were  to  be  compared.  The  algorithm  is  adap¬ 
tive  in  that  a  supervisory  capability  is  provided  to  decide, 
on  the  basis  of  detection  results  and  class  variability, 
if  class  descriptors  should  be  modified.  Algorithm  per¬ 
formance  in  additive,  white,  guassian  noise  was  evaluated 
and  decision  parameters  were  related  to  probability  of 
classification  error. 
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ABSTRACT 


"Phonet  Detection  in  Connected  Speech " 

This  project  generated  phonetic  units,  termed 
"phonets,"  from  digitized  speech  files.  The  time  file 
was  converted  to  feature  space  using  the  Fourier  Transform, 
and  phonet  occurrences  were  detected  using  Minkowski  One 
and  Two  distance  measures.  Phonet  matches  were  detected 
and  ranked  for  each  phonet  compared  against  a  template 
file.  Phonet  short-time  energy  was  included  in  the  output 
files.  An  algorithm  was  developed  to  partition  feature 
space  and  its  performance  was  evaluated. 


I .  INTRODUCTION 
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|  BACKGROUND: 

i 

This  project  concerns  the  generation  and  detection 
of  features  in  connected  speech.  It  is  a  part  of  a  larger 
effort  directed  towards  connected  speech  recognition. 

We  have  divided  the  speech  recognition  problem  into 
six  parts: 

(1)  Digitization  of  the  acoustic  signal, 

(2)  Feature  extraction, 

(3)  Partitioning  of  the  feature  space  into 

classes, 

'  (4)  Feature  detection, 

(5)  Word  recognition  from  detected  features,  and 

(6)  Concept  recognition  from  recognized  words. 
Falkey  (Ref  2)  and  Seelar.ct  (Ref  1)  have  worked  on  Parts 
I,  II,  and  IV  of  the  problem.  This  project  extends 
Seelandt's  work  on  feature  detection  and  offers  an  algorithm 

j  for  partitioning  feature  space  into  classes.  Montgomery 

«  (Ref  10)  has  worked  on  Part  V  of  this  problem  and  Part  VI 

i  has  yet  to  be  addressed. 

In  the  remainder  of  this  chapter,  we  mention  works  in 
*  the  literature  which  influenced  the  development  of  this 

project  and  state  its  scope. 
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Felkey  was  able  to  generate  spectrograms  from  digi¬ 
tized  speech  (Ref  2).  This  capability  was  extended  by 
Seelandt  (Ref  1)  into  a  Speech  Sound  Analysis  Machine 
( SSAM )  which  gave  an  operator  the  ability  to  listen  to 
speech  segments  chosen  from  a  spectrogram  displayed  on  a 
Tektronix  4010-1  Graphics  Terminal.  Seelandt  was  able  to 
use  his  machine  to  choose  phonetic  units  which  he  thought 
might  be  detectable  in  connected  speech 1  He  was  able  to 
construct  template  files  of  these  units  and  then  to  detect 
similar  units  in  connected  speech.  Detection  was  accom¬ 
plished  by  calculating  a  Minkowski  metric  of  order  one 
(Ref  1  and  9)  between  each  phonetic  unit  and  spectrum  of 
consecutive  time  slices  from  the  input  speech.  The  template 
which  matched  the  speech  segment  in  the  sense  of  minimum 
distance  was  deemed  the  detected  phoneme. 

Felkey  and  Seelandt  both  applied  the  Fourier  Transform 
and  extracted  features  in  the  spectral  domain.  Other 
transformations  from  measurement  space  can  be  applied. 
DeSouza  and  Thomson  derived  likelihood  ratio  statistics 
for  the  case  in  which  two  estimated  Linear  Predictive 
Coding  (LPC)  vectors  are  being  compared  (Ref  3).  They  also 
investigated  the  distribution  and  sensitivity  of  other  LPC 
distance  measures. 

Kassam  has  compiled  an  alphabetical  listing  of  papers 
from  the  engineering  literature  on  nonparametric  detection 
theory  and  application  (Ref  4).  Nonparametric  detectors 


are  based  on  statistical  hypothesis  testing  principles 
for  situations  where  parametric  atatiatical  models  cannot 
be  specified  for  the  observation  under  the  null  hypotheeis. 
The  usual  performance  characteristic  is  minimum  Type  I 
error  probability  or  false-alarm  rate.  K ass am  also  states 
that  robust  procedures,  where  not  only  the  false-alarm 
probability  but  also  the  power  of  the  test  is  considered 
in  defining  a  performance  criterion,  may  be  applicable  to 
situations  where  parametric  models  cannot  be  assumed. 

Lainiotis  (Ref  5)  has  studied  adaptive  systems  for 
detection  and  pattern  recognition  as  well  as  for  feature 
extraction.  He  states  that  for  the  supervised  learning 
case,  optimal  adaptive  systems  are  realizable  in  a 
partitioned  form  which  consists  of  a  linear,  nonadaptive 
part  made  up  of  a  bank  of  Kalman  filters,  and  a  nonlinear, 
adaptive  part  made  up  of  probability  computers.  He  offers 
relatively  tight  upper  and  lower  bounds  to  the  probability 
of  error  in  the  feature  extraction  problem  in  terms  of 
functionals  such  as  the  Bhattacharya  coefficient. 

Work  done  at  Johns  Hopkins  University  points  out  that 
the  Past  Fourier  Transform  (FFT)  is  an  extremely  robust 
procedure  for  large  time-band  width  signals  (Ref  7).  The 
report  comments  that  the  FFT  is  robust  enough  to  give  good 
performance  on  an  error  probability  basis  on  random  input 
time  seriea  with  probability  distributions  so  extreme  as 
to  be  unlikely  to  occur  in  nature.  The  price  for  this 
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robustness  is  rsducsd  power;  that  is,  the  probability  of 
■eking  a  correct  binary  decision  conditioned  on  knowledge 
of  the  signal  sent,  according  to  the  report. 

We  decided  to  continue  feature  extraction  in  the 
spectral  domain  for  the  following  reasons: 

(1)  We  could  build  directly  on  the  work  already 
accomplished  by  Seelandt  and  Felkey. 

(2)  We  could  apply  the  already  available  Eclipae 
AP/130  Array  Processor  to  the  task  of  transforming  the 
speech  time  series  to  the  spectral  domain. 

(3)  While  no  one  transformation  into  feature 
space  that  we  found  documented  in  the  literature  was  clearly 
a  best  choice  for  our  purpose,  we  deemed  it  wise  to  sacri¬ 
fice  power  for  the  robustness  of  the  FFT. 

SCOPE: 

The  goals  of  this  project  are: 

(1)  To  generate  phonetic  units,  which  we  term 
phonets  or  observations,  from  speech  files  placed  on  disk. 

(2)  To  detect  occurrence  of  members  of  one  set 
of  phonetic  units,  referred  to  as  phonets,  in  a  contiguous 
file  of  observations  generated  from  speech. 

(3)  To  apply  the  Eclipse  AP/130  Array  Processor 
to  tasks  one  and  two  for  speed. 

(4)  To  build  a  sufficient  degree  of  flexibility 
into  the  phonet  generator  and  distance  computer  to  allow 
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factors,  such  as  FFT  window  point  sisa  and  typa,  filtering, 
and  distance  rule,  to  be  varied  and  their  effect  investi¬ 
gated. 

(5)  To  identify  a  means  for  generating  an 
adequate  phonet  set  for  our  purposes. 

Tasks  one  through  four  are  accomplished  by  development  of 
the  Acoustic  Analysers  an  interactive  software  package 
which  accomplishes  those  tasks.  Task  five  is  accomplished 
by  formulation  of  an  algorithm  for  partitioning  feature 
space  into  template  classes.  Additionally,  some  character¬ 
istics  of  phonets  and  their  interaction  with  observations 
are  illustrated. 

This  work  was  not  concerned  with  acquiring  speech  or 
generating  time  files  containing  speech.  It  takes  the 
techniques  Seelandt  (Ref  1:9-12)  used  as  given  and  relies 
on  speech  files  resident  on  disk. 


II .  TECHNIQUES 


TIME  DOMAIN  PROCESSING; 

In  this  chapter,  we  discuss  techniques  used  to  imple- 
msnt  ths  Acoustic  Anslyssr.  In  this  section,  wo  discuss 
tin#  domain  processing  of  the  digitised  speech,  in  the 
following  sections,  we  describe  the  generation  of  what  we 
call  observations  and  phonets:  spectral  domain  represen¬ 
tations  of  the  fundamental  phonetic  units,  and  the  distance 
computation  between  observations  and  phonets. 

Speech  is  digitized  at  an  8KHz  sampling  rate  and- 
stored  on  disk  in  the  manner  described  by  Seelandt  (Ref  1: 
9-12)  before  being  input  to  the  Acoustic  Analyzer.  Given 
this  input  file,  the  Acoustic  Analyzer  converts  speech  to 
observations  and  then  compares  them  to  a  template  file  of 
phonets.  Time  domain  processing  is  accomplished  in  the 
segment  of  the  Acoustic  Analyzer  which  converts  speech  to 
observations  and  consists  of  the  application  of  a  window 
to  the  time  file.  The  operator  has  control  of  the  following 
window  parameters,  which  are  explained  below: 

(1)  Window  type. 

(2)  Window  size. 

(3)  Overlap  of  the  128  point  window. 


Two  typos  of  windows  oro  available:  a  Hastening  and  a 
root angular  window.  The  Hamming  window  option  waa  included 
bocauao:  (1)  Soolandt  used  it,  and  (2)  spoctral  distor¬ 

tion  caused  by  the  rectangular  window  is  an  undosirod  com¬ 
plication.  Zt  was  desirable  to  include  the  capability  to 
select  a  number  of  different  window  sizes  to  provide  flexi¬ 
bility  in  future  investigations.  It  is  clear  that  a  choice 
of  window  size  involves  a  trade-off  between  spectral 
resolution  and  time  resolution  (Ref  8:260).  It  appears 
that  several  window  sizes  will  need  to  be  tried  and  one 
best  chosen  on  the  basis  of  recognition  results.  Seelandt 
used  a  64  point  window  at  a  sample  rate  of  BKHz  for  time 
slices  SMS EC  long  (Ref  1:120).  Rabiner  and  Schafer  discuss 
the  window  effects  with  regard  to  using  time  slice  energy 
as  an  indicator  of  voiced  or  unvoiced  speech  (Ref  8:122). 
Ihey  recommend  using  a  window  duration  of  from  10-20MSEC. 

At  the  8 KHz  sampling  rate  used  by  Seelandt  to  digitize 
the  speech  files  (Ref  1:10),  a  10-20MSBC  window  is  from 
80-160  points  wide.  Also,  Rabiner  and  Schafer  show  that 
stops  and  other  short-lived  temporal  features  of  speech 
last  roughly  10MSEC  (Ref  8:38-60).  For  these  two  reasons, 
we  consider  the  128  point  window  as  our  primary  choice,  or 
point-of-departure  for  further  investigations  into  the 
question  of  optimum  window  size.  Because  of  concern  that 
a  Hamming  window  shape  coupled  with  16MSEC  window  duration 


Right  miss  mm  temporal  detail.  we  included  the  capability 
to  cause  tha  window  to  incraRant  along  tha  apaach  fila  in 
64  point  incraRants  overlapping  by  64  points.  This 
option  ia  not  available  for  tha  othar  window  sizes  and 
doss  not  significantly  add  to  processing  time. 


OBSERVATION  GENERATION; 

An  observation  is#  basically#  tha  magnitude  of  tha 
Fast  Fourier  Transform  ( FFT)  of  a  segment  of  an  input 
speech  file.  Spectral  shaping,  described  later,  is  applied 
to  the  spectral  components.  It  is  this  shaped  spectrum 
that  is  compared  to  the  phoneta  in  the  template  file:  a 
phonet  being  an  observation  assigned  to  the  template  file. 

Short-time  energy  of  the  speech  signal  is  a  parameter 
that  reflects  amplitude  variations  between  unvoiced  and 
voiced  speech  segments  (Ref  8:120).  Rabiner  and  Schafer 
define  the  short-time  energy  as: 

E  fx(m)  w  (n-m)l  2  (1) 

n  m  ■  -00  ‘ 


with: 


E_  -  short-time  energy, 
n 

£x(ra)j  -  speech  time  sequence. 

£w(n-m)}  -  window  time  sequence  of  finite  length, 
ttiey  define  the  time  dependent  Fourier  Transform  as 
(Ref  8:251): 
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selectable  rate  up  to  a  selectable  corner  frequency,  m^. 
High  frequency  components  are  emphasized  at  a  selectable 
rate  and  comer  frequency,  m2*  Emphasis  rates  are  select¬ 
able  in  units  of  dB/Octave.  Note  that  the  zero  and  L/2 
points  are  zero.  Assuming  neglible  energy  in  the  speech 
signal  at  4KHz.  and  assuming  a  DC  block,  the  observations 
computed  by  our  Acoustic  Analyzer  are  the  output  of  the 
linear  system  as  shown  in  Figure  2. 


FIGURE  2.  Linear  System  Representation  of  Spectral  computer 
In  Acoustic  Analyzer.  £x(n)J  is  the  speech 
input. 

By  Parseval's  theorem  (Ref  9:36),  the  energy  in  the 
output  sequence,  E*,  is  related  to  the  energy  in  the  input 
speech,  Er,  under  assumptions  that  there  is:  (1)  no  DC 
value  to  the  speech  input,  and  (2)  no  speech  energy  at  or 
above  4KHz.  Hence,  the  observation  energy,  E*,  computed 
by  the  Acoustic  Analyzer,  can  be  used  as  an  indication  of 
voiced/unvoiced  speech  in  the  same  manner  as  the  short-time 
energy.  This  value  can  be  used  by  a  word  recognition 


machine  to  aid  in  that  process.  The  output  file  at  this 
point  contains i  for  each  observation*  the  observation 
energy  as  the  first  component  and  spectral  components  one 
through  L/2  -1  in  the  second  through  L/2  component  posi¬ 
tions  . 

For  the  same  reason  Seelandt  did  (Ref  1:121),  we 


normalized  the  energy  in  each  observation  to  a  constant 
value.  This  provided  a  measure  of  automatic  gain  control, 
in  that  if  the  only  difference  between  two  speech  segments 
was  their  energy  content,  the  components  in  observations 
produced  from  them  would  be  of  approximately  the  same 
amplitude.  Thus,  variability  caused  by  different  speech 
amplitudes  would  be  reduced.  To  normalize  the  observation 
energy,  we  scaled  each  spectral  component  by  the  factor 

io4/VeJ  • 

We  observed  the  effectiveness  of  obtaining  automatic 
gain  control  through  energy  normalizing  by  including  the 
capability  to  scale  the  spectral  components  by  the  factor 
104/E*,  rather  than  energy  normalizing.  We  generated  two 
time  files  containing  tones  at  200Hz,  lKHz,  and  2.5KHZ, 
each  at  an  amplitude  of  50  units  in  File  A  and  200  units 
in  File  B.  Figure  3  is  an  observation  from  File  A  and 
Figure  4  is  from  File  B;  both  energy  normalized.  In  these 


and  the  next  two  figures,  the  value  plotted  in  position 
zero  along  the  horizontal  axis  is  the  energy  value,  V4- 
One  can  see  that  the  amplitudes  of  the  spectral  components 
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are  nearly  the  same  in  the  two  observations.  The  amplitudes 
of  the  spectral  components  differ  markedly  between  the  next 
two  figures  which  were  scaled  by  the  factor  104/En. 

Figure  5  was  computed  from  File  A  which  contained  tones 
at  an  amplitude  of  200  units  and  Figure  6  was  computed 
from  File  B  which  contained  tones  at  an  amplitude  of  50 
units.  One  can  also  see  that  energy  normalization  removes 
from  the  observation  spectrum  all  of  the  variability  caused 
by  speech  amplitude  differences. 

Energy  normalization  has,  at  least,  one  other  conse¬ 
quence:  it  limits  the  maximum  Euclidean  distance  between 

any  two  observations.  To  see  this,  suppose  the  observa¬ 
tions  have  been  energy  normalized  by  having  their  com¬ 
ponents  scaled  by  the  factor  K/tJ E  ,  wneie  E  is  the  energy 
in  observation  being  normalized.  We  write  each  observation 
as  a  vector  and  note  that  the  energy  in  the  difference 
between  any  two  observations  is: 


*1 


(4) 


2 

Now  the  energy  in  each  observation  is  K  so  the  Euclidean 
distance  is: 


|*1  “  *2|^  K  (5) 

So  the  Euclidean  distance  is  upper  bounded  by  VT  times 


the  scale  factor  constant,  K 
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PHONET  GENERATION: 


Phonets  are  selected  from  a  set  of  observations  com¬ 
puted  by  the  Acoustic  Analyzer  (AA)  using  the  spectrum 
option.  We  presently  have  no  systematic  means  for  selecting 
an  adequate  set  of  phonets.  Such  a  set  should  span  the 
observation  space  in  that  any  observation  can  be  seen  to 
resemble  at  least  one  phonet  in  a  detectable  way.  It  should 
also  be  discriminatory  in  that  any  observation  should 
resemble  only  one  phonet  better  than  all  others  to  a  detect¬ 
able  degree.  One  last  property  that  an  adequate  phonet  set 
should  have  is  that  of  compactness,  in  that  the  number  of 
phonets  should  not  be  so  great  as  to  impose  prohibitive 
computational  burdens.  Given  any  phonet  set,  determining 
if  it  has  these  three  basic  properties,  requires  considera¬ 
tion  of  the  distance  rule  to  be  used. 

DISTANCE  COMPUTATION: 

Seelandt  used  the  Minkowski  One  (Ml)  distance  and  a 
smallest  distance  determination  (Ref  1:53-57)  to  choose 
five  closest  phonetic  units  to  each  observation  run  through 
his  distance  computer.  He  used  this  measure  because  it  was 
convenient  to  compute.  Our  Acoustic  Analyzer  includes 
this  distance  measure  as  one  option.  We  also  included  the 
Euclidean,  or  Minkowski  Two  (M2)  distance  (Ref  9:232). 

One  option  of  our  Acoustic  Analyzer  is  a  distance  computer 
which  uses  either  the  Ml  or  M2  measure  to  calculate  a 
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symmetric  distance  matrix.  <^.  The  element  <P(i,  j)  is  the 

th 

distance  between  the  i  observation  and  the  j  phonet. 

The  other  distance  option  in  our  AA  uses  either  distance 
measure  to  choose  a  selectable  number  of  closest  phonets 
to  each  observation.  In  the  file  for  each  observation  is 
the  observation  number  as  the  first  element,  the  observa¬ 
tion  energy  as  the  second  element,  an  ordered  pair  of  (phonet 
number,  corresponding  maximum  distance),  in  the  next  two 
elements  and  the  selected  number  of  ordered  pairs  (phonet 
number,  corresponding  minimum  distance).  These  last  ordered 
pairs  are  themselves  ordered  minimum  distance  choice-to- 
maximum.  Files  for  consecutive  observations  are  contiguously 
arranged  in  the  output  file. 

We  chose  to  implement  the  M2  distance  measure  because 
it  is  the  Euclidean  distance  in  a  geometrical  sense  in  the 
absence  of  noise.  We  implemented  the  Ml  rule  to  provide 
the  capability  to  compare  results  obtained  using  the  Acoustic 
Analyzer  with  those  using  Seelandt's  Ml  implementation.  It 
is  acknowledged  that  other  distance  measures  and  more 
elaborate  decision  rules  will  be  needed  to  achieve  a  speci¬ 
fied  error  probability  in  noise.  This  implementation  assumes 
a  completely  deterministic  system  without  noise  of  any  kind. 
Noise  is  introduced  into  the  analysis  of  the  feature  space 
cluster  algorithm  described  in  the  next  chapter. 
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III.  RESULTS 


THE  ACOUSTIC  ANALYZER  (AA)  : 

One  objective  of  this  project  was  to  implement  a 
machine  like  Seelandt's  (Ref  1)  Speech  Sound  Analysis 
Machine  (SSAM)  on  the  Eclipse  AP/130  Array  Processor,  taking 
advantage  of  its  fast  computational  capability.  This  work 
resulted  in  delivery  of  a  fast,  flexible  software  package 
with  the  required  features: 

(1)  A  flexible  spectral  computing  option  which 
transforms  time  slices  of  speech  into  observations. 

(2)  A  flexible  distance  computer  which  calculates 
a  symmetric  distance  matrix  between  each  observation  and 
phonet. 

(3)  Another  distance  computer  which  chooses  a 
selectable  number  of  best  matches  from  a  distance  matrix. 

(4)  A  graphics  capability  which  allows  one  to 
view  plots  of  speech  files,  observation  files,  phonet  files, 
and  distance  files. 

Ifre  ability  to  select  parameters  such  as  those  which 
define  the  FFT  time  window,  define  processing  to  be  done 
on  the  transformed  spectral  units,  and  characterize  the 
detection  scheme,  has  been  included  in  our  Acoustic  Analyzer. 
We  chose  to  implement  a  128  point,  overlapped  FFT  window. 


! 

In  the  discussion  on  observation  generation,  we  saw  that 
normalizing  the  energy  in  each  observation  removed  much  of 
the  variability  from  speech  amplitude  variation,  in  the 
next  section,  we  will  illustrate  effects  of  other  parameter 
selections. 

PARAMETER  VARIATION; 

In  this  section,  we  illustrate  observations  and  dis¬ 
tances  computed  using  several  parameter  settings.  We 
use  only  one  speech  file,  CT56.0B,  which  is  an  adult  male 
saying  the  words  "five"  and  "six."  It  is  beyond  the  scope 
of  this  project  to  illustrate  parameter  variation  effects 
for  a  large  number  of  speech  files?  but,  there  is  reason 
to  believe  that  such  effects  would  be  similar  to  those 
illustrated  here. 

Seelandt  (Ref  1}  identified  the  following  parameters 
as  requiring  study  in  order  to  determine  optimum  settings: 

(1)  The  number  of  consecutive  time  slices  grouped 
together  to  represent  a  phoneme.  Seelandt 's  phonetic  units 
were  five  time  slices  long. 

(2)  The  FFT  window  size.  Seelandt  used  a  64 
point  window. 

(3)  Rate  at  which  the  speech  time  signal  is 
sampled  when  digitized.  Seelandt  used  speech  sampled  at 
8KHz. 

(4)  Window  shape:  rectangular,  Hamming,  or 
other.  Seelandt  used  the  Hamming  window. 
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(5)  Spectral  shaping:  pre-emphasis  to  emphasize 
frequency  components  above  a  corner  frequency,  and  de¬ 
emphasis  to  attenuate  frequency  components  below  a  corner 
frequency. 

(6)  Spectral  domain  thresholding.  Seelandt 
normalized  the  energy  in  the  spectrum  in  each  time  slice 
when  that  energy  was  above  a  threshold.  He  attenuated  the 
spectrum  in  time  slices  which  contained  less  energy  than 
the  threshold. 

PHONEME  LENGTH: 

We  decided  to  treat  phonetic  units  one  time  slice 
long,  rather  than  form  five  time  slice  phonemes,  for  the 
following  reasons.  First,  it  is  more  convenient  to  process 
single  unit  phonets  than  multiple  unit  phonemes.  Second, 
the  dynamic  behavior  of  speech  segments  would  be  visible 
with  more  resolution  in  a  plot  of  aistance-versus-phonet 
number.  Third,  until  the  dynamic  behavior  of  adjacent 
short-time  speech  units  is  better  characterized,  we  believed 
that  it  would  not  be  prudent  to  average  their  effect  in  a 
distance  measure.  We  reasoned  that  it  would  be  better  to 
preserve  the  time  slice-by-time  slice  variation.  Then  when 
short-time  speech  dynamics  are  better  understood,  one  can 
'  simply  reformat  spectral  and  distance  files  generated  by 

'  our  Acoustic  Analyzer  to  accommodate  any  desired  phoneme 

length. 
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Computing  distances  on  single  time  slice  phonets,  we 
found  that  adjacent  time  slices  were  similar  in  the  midst 
of  a  vowel  and  not  similar  in  unvpiced  fricatives.  From  a 
spectrogram  of  file  CT56.0B,  we  knew  that  the  speaker  was 
uttering  the  word  "five"  during  time  slices  40-80  and  the 
“si"  part  of  "six"  during  time  slices  90-100.  The  spectro¬ 
gram  is  provided  in  Appendix  A.  We  used  a  128  point,  non- 
overlapped  Hamming  window  with  lOdB/Octave  de-emphasis 
ending  at  300Hz  and  lOdB/Octave  pre-emphasis  starting  at 
500Hz.  We  computed  the  M2  distance  between  the  file  and 
itself  for  time  slices  41-102:  that  is,  we  ran  observations 
41-102  against  phonets  41-102  of  the  same  file.  Plots  of 
observations  41,  51,  61,  71,  81,  and  91  against  phonets 
41-102  are  shown  in  Figures  8,  10,  12,  13,  14,  and  15, 
respectively.  The  similarity  of  adjacent  observations  in 
the  vowel  regions  in  Figures  10  and  12  are  apparent,  as  is 
the  dissimilarity  of  adjacent  observations  elsewhere.  This 
effect  is  illustrated  on  an  expanded  scale  in  Figures  9  and 
11,  which  show  distances  between  the  observations  45  and  55, 
respectively,  plotted  against  phonets  41-70. 

We  have  seen  this  effect  every  time  we  have  looked  for 
it.  Adjacent  time  slices  resemble  each  other  in  both  M2 
and  Ml  space  much  more  in  regions  of  vowels  when  a  spectro¬ 
gram  shows  clear  formant  structure  than  in  regions  where 
a  spectrogram  shows  little  formant  structure.  This  formant 
structure  is  visible  in  plots  of  the  observations. 
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Copy  ovoilablo  to  DTIC  does  not 
permit  fully  legible  reproduction 


FIGURE  10.  Distance  plot  of  observation  number  51  against  41  through  102 
M2  rule,  128  point,  nonoverlapped  Hamming  window,  energy 
normalized. 


ThNCE  plot 


FIGURE  12.  Distance  plot  of  observation  61  against  41  through  102. 

M2  rule,  128  point,  nonoverlapped  Hamming  window,  energy 
normalized . 


DISTANCE  PLOT 
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FIGURE  15.  Distance  plot  of  observation  91  against  41  through  102. 

M2  rule,  128  point,  nonover lapped  Hamming  window,  energy 
normalized . 


Figures  16  and  17  are  plotted  using  the  spectral  plot 
option  of  our  Acoustic  Analyzer,  and  are  observations  45 
and  47  from  the  observation  file  above.  Though  these 
observations  are  separated  by  32MSEC,  their  resemblance 
is  clear. 

FFT  WINDOW  SIZE: 

We  have  discussed  our  choice  of  FFT  window  size  in 
Section  II  of  this  report.  A  128  point  window  provides 
more  frequency  resolution  than  does  a  64  point  window,  and 
leads  to  useful  behavior  of  short-time  energy  with  regard 
to  voiced/unvoiced  speech  (Ref  8).  We  also  found  when  we 
computed  distances  between  an  observation  file  and  itself 
in  voiced  vowel  regions  of  pronounced  formant  structure, 
that  adjacent  observations  had  similar  formant  structure 
using  a  64  point  window  as  with  a  window  size  of  128  points. 
To  illustrate,  observations  were  computed  using  a  64  point 
Hamming  window.  Observations  were  de-emphasized  at  lOdB/ 
Octave  below  300Hz  and  pre-emphasized  at  lOdB/Octave  above 
500Hz.  The  formant  structure  can  be  seen  to  develop  in 
the  sequence  of  plots  of  observation  number  86,  87,  88, 

89,  and  90;  Figures  18,  19,  20,  21,  and  22,  respectively. 
These  observations  are  in  the  same  speech  segment  as  Figures 
8-17;  there  are  twice  the  number  of  observations  in  a  time 
segment,  here,  at  64  points  per  time  slice.  The  M2  dis¬ 
tance  measure  can  be  seen  to  develop  broad  valleys  in 
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FIGURE  17.  Spectral  plot  of  observation  number  47.  One-hundred  twenty-' 
(128)  point,  nonover lapped  Hamming  window,  energy  normalized 


FIGURE  18.  Spectral  plot  of  observation  number  86.  Sixty-four  (64)  point 
Hamming  window,  energy  normalized. 
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FIGURE  20.  Spectral  plot  of  observation  number  88.  Sixty-four  (64)  point 
Hamming  window,  energy  normalized. 


FIGURE  21.  Spectral  plot  of  observation  number  89.  Sixty-four  (64)  point 
Hamming  window,  energy  normalized. 


regions  of  most  pronounced  formant  structure  in  Figures 
23 ,  24,  25,  26,  27,  and  28.  These  are  plots  of  observa¬ 
tions  number  85,  90,  115,  120,  125,  and  130,  respectively, 
all  against  phonet  numbers  81-144. 

SAMPLING  RATE,  WINDOW  SHAPE,  AND  EMPHASIS: 

We  used  speech  files  digitized  at  an  8KHz  sampling 
rate  and  12  bit  quantization  because  it  was  the  highest 
sampling  rate  available.  It  is  well  known,  of  course,  that 
speech  sampled  at  this  rate  is  perfectly  intelligible  when 
reconstructed.  Seelandt  describes  the  digitizing  aparatus 
(Ref  1).  Here  we  discuss  the  advisability  of  changing  that 
rate.  We  also  discuss  window  shape. 

The  capability  or  our  Acoustic  Analyzer  to  calculate 
FFTs  of  window  size  512  or  1024  points  allows  one  to  double 
or  quadruple  the  sampling  rate  while  maintaining  the 
behavior  of  short-time  energy  discussed  by  Rabiner  and 
Schafer  (Ref  8)  and  in  Section  II  of  this  report.  Using 
the  Array  Processor,  the  computational  burden  of  large  win¬ 
dows  is  light.  The  advantages  would  not  only  be  an  increase 
in  spectral  resolution,  but  also  a  decrease  in  spectral 
distortion  due  to  the  longer  window.  The  effects  of  window 
shape  become  more  important  as  pre-emphasis  is  applied  to 
high  frequencies,  especially  so  when  a  large  portion  of  the 
energy  in  speech  is  concentrated  in  the  lower  frequency 
range.  In  this  case,  the  sidelobes  from  the  low  frequency 


FIGURE  23.  Distance  plot  of  observation  85  against  81  through  144. 

M2  rule,  64  point  Hamming  windov,  energy  normalized. 


FIGURE  24.  Distance  plot  of  observation  90  against  81  through 
rule,  64  point  hamming  window,  energy  normalized. 
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FIGURE  26.  Distance  plot  of  observation  number  120  against  81  through  144. 
M2  rule,  64  point  Hamming  window,  energy  normalized. 
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energy  are  emphasized.  Sivian  showed  that  formant  amplitude 
rolled-off  at  approximately  lOdB  per  Octave  relative  to  the 
fundamental  (Ref  11).  Kuligowski  has  shown  the  same  (Ref 
12:20.  21).  If  pre-emphasis  is  applied  to  the  speech 
spectrum  to  compensate  for  this  roll-off,  then  sidelobes  of 
the  window  spectrum  excited  by  large  amplitude  fundamental 
and  first  formant  components  have  a  20dB  per  Octave  advantage 
over  the  high  formants:  lOdB  per  Octave  because  the  formant 
is  down  that  far  and  lOdB  per  Octave  because  the  sidelobe 
is  emphasized  that  much.  De-emphasis  of  the  fundamental 
in  the  spectral  domain  does  not  affect  these  sidelobes 
because  they  are  generated  in  the  transformation  and  de¬ 
emphasis  is  applied  afterwards.  They  can  be  lowered  by: 

(1)  High-pass  filtering  before  the  FFT  to  decrease 
the  energy  in  the  fundamental,  and  hence  in  the  sidelobes. 

(2)  Using  the  longest  possible  window  size  and 
least  distorting  window  shape  while  maintaining  the  10MSEC- 
to-20MSEC  window  duration  for  short-time  energy. 

The  first  solution  is  difficult  to  apply.  Kuligowski  has 
shown  (Ref  12:20,  21)  that  the  location  of  the  fundamental 
varies  by  as  much  as  an  octave  over  the  vowels  he  listed, 
and  depending  on  speaker  sex  and  age.  Further,  for  several 
vowel-speaker  combinations,  the  fundamental  of  one  vowel- 
speaker  combination  is  at  the  same  spectral  location  as  the 
first  formant  at  other  vowel-speaker  combinations.  An 
adaptive  filter  before  the  FFT  would  be  necessary  to 
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attenuate  the  fundamental  prior  to  transformation  without 
attenuating  the  first  formant.  The  second  solution  can  be 
implemented  with  enough  memory  and  a  fast  enough  analog-to- 
digital  converter.  At  a  1024  point  window  size,  a  sampling 
rate  of  68KHz  would  provide  a  window  duration  of  15MSEC. 

As  has  been  previously  noted,  the  FFT  is  a  robust 
feature  (Ref  7).  It  may  be  sufficiently  robust  for  an  8KHz 
sampling  rate  and  a  128  point  Hamming  window  with  the  obser¬ 
vation  pre-emphasized  at  lOdB  per  Octave.  The  measure  of 
sufficiency  is  in  the  way  the  features  cluster,  and  an 
adequate  feature-space  partitioning  algorithm  for  this 
problem  has  not  yet  been  implemented. 

It  is  clear  that  the  fundamental  should  be  de-emphasized 
in  the  observation  spectrum.  It  is  a  high  energy  feature 
which  is  highly  speaker  dependent.  The  problem  with  de¬ 
emphasizing  it  is  that  for  some  vowel-speaker  combinations, 
it  is  located  at  or  near  the  first  formant  of  other  vowel- 
speaker  combinations  (Ref  12:20,  21).  We  chose  to  de- 
emphasize  at  lOdB/Octave  ending  at  a  corner  frequency  of 
300Hz  as  a  compromise  between  allowing  too  much  fundamental 
energy  to  remain  in  the  observation,  and  reducing  the 
amount  of  energy  in  the  important  first  formant. 

THRESHOLDING: 

The  last  parameter  Seelandt  identified  for  further 
study  is  a  threshold  below  which  phonetic  units  should  be 


attenuated  when  their  energy  content  is  low.  When  the 
energy  in  a  time  slice  did  not  exceed  a  preset  threshold, 
the  spectral  components  were  attenuated  rather  than  energy 
normalized.  Seelandt  did  this  to  prevent  vectors  which 
consist  predominately  of  background  noise  from  entering 
into  the  decision  process  (Ref  1:121,  122). 

We  found  that  unvoiced  fricative  sounds,  especially 
those  in  sounds  involving  the  letter  "f,"  are  low  energy 
sounds  and  hard  to  distinguish  from  noise  in  a  spectrogram. 
Phonet  plots  did  reveal  structure  in  these  sounds  that 
might  distinguish  them  from  white  noise  (flat  spectrum), 
or  other  background  noise  models.  Distance  measures  com- 
puted  to  phonets  which  represent  background  noi=e,  coupled 
with  knowledge  of  phonet  and  observation  energy,  may  permit 
a  phonet- to-word  translation  machine,  such  as  Montgomery's 
(Ref  10),  to  distinguish  between  unvoiced  fricative  sounds 
and  noise,  and  may  aid  in  detecting  word  boundaries  as 
well.  Figures  29-55  illustrate  observations  taken  from 
regions  of  the  speech  file  CT56.SP  known  to  contain  noise 
as  well  as  the  words  "five"  and  "six."  The  plots  are  in 
dB  versus  spectral  component  number.  Component  number  zero 
is  the  observation  energy,  also  in  dB.  Figures  29-34  are 
from  a  region  where  there  was  no  speaker  sound.  Compare 
those  figures  with  Figures  35-42  which  are  from  a  region  of 
the  file  where  the  speaker  begins  the  "five"  sound.  It  is 
difficult  to  detect  where  the  "f"  sound  begins,  either  by 
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studying  the  spectral  distribution  or  by  noting  the  energy. 
One  can  see  that  as  the  characteristic  formant  structure 
associated  with  the  vowel  "i“  sound  develops ,  the  energy 
increases  from  roughly  63dB  to  roughly  lOOdB  in  the  "i" 
sound.  But  it  does  not  appear  that  energy  content  is  a 
useful  means  for  discriminating  between  noise  and  the  "f" 
sound. 

Figures  43-50  were  taken  from  the  word  "five"  from 
well  into  the  vowel  sound  "i"  through  the  "v"  sound  on 
into  a  pause  between  the  words  "five"  and  "six."  One  can 
see  the  spectral  distribution  change  from  the  "i"  formant 
structure  through  a  "v"  sound  structure  of  relatively  flat 
spectrum  but  moderate  energy,  to  a  nonspeaker  background 
noise  of  flat  spectrum  and  low  energy.  Figures  51-56  show 
a  sequence  of  observations  from  the  speakerless  background 
region  between  the  words  to  well  into  the  "s"  in  “six." 

Note  that  as  the  characteristic  upward  slope  ( lef t-to-right) 
of  the  spectral  distribution  of  the  "s"  sound  develops,  the 
energy  increases  roughly  25dB.  One  can  see  in  this  overall 
sequence  of  figures  that  energy  is  a  poor  discriminant  be¬ 
tween  noise  and  some  unvoiced  sounds.  We  conclude  that 
energy  is  an  unreliable  discriminant  between  background 
noise  and  unvoiced  speech. 

The  question  of  how  unvoiced  sounds  may  be  detected 
in  noise  needs  to  be  addressed.  It  is  beyond  the  scope  of 
this  project  to  do  so  because  it  remains  to  be  seen  if 
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FIGURE  29.  Spectral  plot  observation  number  21,  128  point  overlapping 

Hamming  window,  lOdB  de-emphasis  to  300HZ,  lOdB  pre-emphasi s  from 
500HZ,  energy  normalized. 
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FIGURE  30.  Spectral  plot  observation  number  24,  128  point  overlapping 

Hamming  window,  lOdB  de-emphasis  to  300HZ,  lOdB  pre-emphasis  from 
500HZ,  energy  normalized. 
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FIGURE  31.  Spectral  plot  observation  number  31,  128  point  overlapping 

Hamming  window,  lOdB  de-empliasis  to  300JIZ,  XOdB  prc-cmphasis  from 
500HZ,  energy  normalized. 
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FIGURE  32.  Spectral  plot  observation  number  34,  128  point  overlapping 

Hamming  window,  lOdB  de-emphasis  to  300HZ,  lOdB  pre-emphasis 
from  500HZ ,  energy  normalized. 
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FIGURE  36.  Spectral  plot  observation  number  85,  128  point  overlapping 

Hamming  window,  lOdB  de-emphasis  to  300IIZ,  lOdB  pre-emphasis 
from  500HZ,  energy  normalized. 
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FIGURE  37.  Spectral  plot  observation  number  87,  128  point  overlapping 

Hamming  window,  lOdB  de-emphasis  to  300HZ,  lOdB  pre-emphasis 
from  500HZ,  energy  normalized. 
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FIGURE  38.  Spectral  plot  observation  number  88,  128  point  overlapping 
Hamming  vjindcv?,  IQdB  de-emphasis  to  300HZ,  lOdB  pre-empha 
from  500HZ,  energy  normalized. 
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FIGURE  39.  Spectral  plot  observation  number  90,  128  point  overlapping 

Hamming  window,  lOdb  de-emphasis  to  300HZ,  lOdB  pre-emphasis 
from  500HZ,  energy  normalized. 
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FIGURE  AO.  Spectral  plot  observation  number  95,  128  point  overlapping 

Hamming  window,  lOdB  de-emphasis  to  300HZ,  lOdB  pre-emphasis 
from  500HZ,  energy  normalized. 
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FIGURE  43.  Spectral  plot  observation  number  120,  128  point  overlapping 
Hamming  window,  lOdB  de-emphasis  to  300FZ,  lOdB  pre-emphasis 
from  500HZ,  energy  normalized. 
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FIGURE  44.  Spectral  plot  observation  number  130,  128  point  overlappin 
Hamming  window,  lUdB  de-emphasis  to  300HZ,  lOdB  pre-emphas 
from  500HZ ,  energy  normalized. 
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FIGURE  45.  Spectral  plot  observation  number  135,  128  point  overlapping 
Hamming  window,  lOdB  de-emphasis  to  300HZ,  lOdB  pre-empliasi 
from  500HZ,  energy  normalized. 
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FIGURE  48.  Spectral  plot  observation  number  150,  128  point  overlapping 
Hamming  window,  10dB  de-emphasis  to  300HZ,  lOdB  pre-emphasis 
from  500HZ,  energy  normalized. 
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FIGURE  50.  Spectral  plot  observation  number  170,  128  point  overlapping 
Hamming  window,  lOdB  de-emphasis  to  JOUHZ.,  lOdB  pre-emphasi 

from  500HZ,  energy  normalized. 
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FIGURE  5?.  Spectral  plot  observation  number  195,  128  point  overlapping 
Hamming  window,  1  Gun  de-emphasis  Lu  3GGHZ,  iGuo  pi'e-erapnasL 
from  500HZ,  energy  normalized. 
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FIGURE  54.  Spectral  plot  observation  number  200,  128  point  overlapping 
Hamming  window,  lOdB  de-emphasis  to  JUUHZ,  iUdB  pre-emphasis 
from  500HZ,  energy  normalized. 
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FIGURE  55.  Spectral  plot  observation  number  201,  128  point  cverlappin 
Hamming  window,  lOdB  de-emphasis  to  300HZ,  lOdB  pre-empnas 
from  500HZ.  enercv  normalized. 
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FIGURE  56.  Spectral  plot  observation  number  202,  128  point  over 
Hamming  window,  lOdB  de-emphasis  to  300HZ,  lOdB  pre- 
from  500HZ,  energy  normalized. 


feature  space  will  cluster  for  unvoiced  sounds  and  noise 
sources.  It  can  be  argued,  however,  that  success  of  a 
scheme  for  detecting  phonetic  units  in  noise  will  depend 
not  so  much  on  an  energy  threshold,  but  on  how  well  phonets 
representing  noise  and  speech  separate,  and  on  noise  can¬ 
cellation  techniques  employed.  Several  clustering  algorithms 
(Ref  13)  will  need  to  be  studied,  and  perhaps  implemented, 
to  quantify  phonet  detection  in  noise.  A  cluster  algorithm 
is  presented  in  the  next  section.  But  it  seems  too  early, 
at  this  point,  to  settle  on  a  threshold  criterion  such  as 
Seelandt ' s. 

A  reliable  discriminant  between  unvoiced  speech  and 
noise  may  be  one  based  on  clustering  properties  of  obser¬ 
vations  about  selected  phonets.  We  would  not  expect  dis¬ 
tances  from  a  white  noise  background  model  to  cluster,  but 
nonwhite  noise  may. 

The  distribution  of  similar  observations  in  several 
speech  segments  is  plotted  in  Figures  57-68.  In  each  plot, 

M2  distances  are  plotted  from  a  particular  observation  to 
each  observation  in  the  speech  file.  The  smallest  five 
and  the  largest  distances  are  plotted  along  with  the  energy 
of  the  particular  observation.  The  position  along  the 
horizontal  axis  at  which  a  distance  is  plotted  corresponds 
to  the  number  of  the  observation  to  which  that  distance  was 
computed.  One  of  the  minimum  distances  will  be  zero 
corresponding  to  a  perfect  match.  The  energy  value  is 
plotted  in  position  zero. 
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Figures  57  and  58  plot  the  benavior  of  two  observations 


from  the  nonspeaker  region  before  the  word  "five"  and 
illustrate  the  lack  of  clustering  characteristics  of  low 
energy  speakerless  regions.  Figures  63-65  plot  distance 
choices  from  observations  in  the  speakerless  segment  between 
the  words  "five"  and  "six"  and  again  show  the  typical  lack 
of  clustering. 

Figures  59  and  60  plot  distance  choices  for  observa¬ 
tions  in  the  "f"  segment  and  in  the  “fi"  segment,  respec¬ 
tively.  Observations  from  the  "f"  sound  of  this  speaker 
did  not  cluster  well:  but,  as  the  voiced  "i"  sound  came  in, 
distance  choices  clustered  much  tighter.  Characteristic 
of  the  tight  cluster  in  high  energy  vowel  sounds  is  the  plot 
in  Figure  62.  We  were  surprised  to  find  that  distance 
choices  from  observations  in  "s"  segments  clustered  tightly. 
Figures  66,  67,  and  68  plot  choices  from  observations  in 
the  "s"  sound  in  the  word  "six."  From  a  spectrogram  of 
the  speech  segment,  we  could  see  that  these  observations 
precede  a  strong  "i"  influence,  and  spectral  plots  from  the 
region.  Figures  54,  55,  and  56,  do  not  show  the  characteris¬ 
tic  "i"  structure.  These  are  "s"  sounds  of  moderate  energy 
content  which  cluster. 

Cm1  tl'l**  Mt<J| »>l  ,  w4*  h, <'*«•>  | •*«-*<«*•  1 1- w(t  l  1 1 1- 

Acgugtic  Analyzer  and  illustrated  observations  and  tiis  • 
tdnCod  CoillpUtcd  udillg  vulUc3  We  have  assigned  tu  .jix 
parameters.  We  have  discussed  these  parameters  in  both 
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FIGURE  64.  Five-Best  Distances.  Plot  of  observation  energy  in  ordinate 
position  zero,  maximum  and  f ive-smallest  distances  from 
observation  number  181,  M2  rule. 
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FIGURE  67.  Five-Best  Distances.  Plot  of  observation  energy  in  ordinate 
position  zero,  maximum  and  fiV6*Suisllcst  distances  from 
observation  number  200  ,  M2  rule. 
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this  and  tha  previous  chaptsr.  We  saw  that  adjacent 
observations  nay  lie  approximately  the  same  distance  from 
neighboring  observations  in  segments  of  speech  which  con¬ 
tain  high- to- moderate  energy  and  pronounced  formant 
structure*  We  saw  this  behavior  for  voiced  vowels  and  for 
the  "i*  sound*  We  conclude  from  this  observed  behavior 
that  observations  generated  from  speech  may  cluster  (Ref 
13)  in  feature  space.  That  is*  it  should  be  possible  to 
classify  observations  without  a  priori  knowledge  based  on  a 
distance  measure  computed  between  them.  In  the  next  section, 
we  describe  an  algorithm  designed  to  do  just  that. 

CLUSTER  ALGORITHM: 

In  the  previous  section,  we  illustrated  observations 
and  their  behavior  under  a  distance  measure.  In  this  sec¬ 
tion,  we  describe  an  adaptive  algorithm  designed  to  cluster 
feature  space  into  classes  of  observations.  Classification 
is  made  on  the  basis  of  a  decision  rule  without  the  benefit 
of  a  priori  knowledge.  The  algorithm  is  adaptive  in  that  a 
supervisor  controls  the  modification  of  class  descriptions 
in  response  to  the  classification  results. 

Pal,  et.  al.,  described  an  adaptive  model  (ilef  6)  for 
computer  recognition  of  vowel  sounds  using  the  first  three 
formants  as  features.  Their  method  uses  a  training  pro¬ 
cedure  for  self-supervised  learning  and  a  maximum  value  of 
fuzzy  membership  function  as  the  basis  of  recognition.  In 
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the  remainder  of  this  auction,  wo  will  present  an  algorithm 
baa ad  on  tho  abova  model*  but  roodifiod  to  partition  tha 
apaco  of  obaarvationa  in  noiao  into  clusters. 

Tha  nodal  for  our  clustor  algorithm  ia  ahown  in.  Figure 
69.  It  ia  rodrawn  from  Rafarenca  6.  It  usaa  a  claaaifiar 
which  calculataa  a  aufficiant  atatiatic  betwaan  a  aat  of 
phonata,  R,  and  an  input  observation,  x.  The  obaervation 
ia  aaaignad  to  tho  phonot  claaa  if  the  aufficient  atatiatic 
ia  amall  enough;  othorwiae,  tha  obaarvation  defines  a  now 
claaa . 
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FIGURE  69.  Model  of  an  adaptive  recognition  scheme  adapted 
from  reference  6.  The  classifier  assigns  the 
input  vector  to  the  class  to  which  it  conforms 
most  closely*  The  supervisor  decides  if  tha 
class  descriptors  should  be  modified. 
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Each  sample,  X^,  In  the  sample  space.' A*  is  an  obser¬ 
vation  computed  by  our  Acoustic  Analyzer.  Ws  aaaums  that 
thsss  samples  have  a  central  tendency  about  which  they 
cluster;  this  central  tendency  is  approximated  by  the 
ensemble  average  of  all  observations  which  cluster  tightly 
about  it* 

The  sample  space  is  partitioned  into  a  sequence  of 
classes  £a^  j  .  Bach  class  is  described  by  its  central 
tendencyf  which  we  call  a  phonetr  by  its  variance,  which 
is  the  second  central  moment  of  all  observations  which 
cluster  tightly  about  the  central  tendency;  by  an  inner 
cone  parameter  which  quantifies  the  adjective  "tight"  as 
applied  to  a  cluster;  and  by  an  outer  zone  parameter  which 
determines  the  outer  boundary  of  the  class.  These  class 
descriptors  are  calculated  iteratively  as  the  sample 
space  is  partitioned. 

In  this  model,  the  supervisor  assumes  that  observations 
have  a  central  tendency  about  which  they  cluster,  it  also 
assumes  that  if  an  observation  lies  close  enough  to  a 
central  tendency,  then  the  probability  of  misclassification 
is  tolerable* 

Accordingly,  an  annulus  is  constructed  about  each 
central  tendency  and  each  observation  is  compared  to  the 
boundaries  of  this  annulus.  If  the  observation  is  close 
enough  to  the  central  tendency  to  fall  within  the  inner 
radius,  then  it  is  permitted  to  modify  the  description  of 
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that  class.  If  it  does  not  fall  within  the  outer  radius 
of  any  annulus,  then  a  new  class  is  created  to  which  this 
observation  is  assigned. 

In  operation,  a  tentative  classification  of  the  ith 

observation,  Xi,  to  the  jth  class,  A  ^ ,  is  based  on  the 

minimum  value  of  the  sufficient  statistic  D(i,  j):  that 

is,  tentatively  assign  x.  to  class  A.  if  D(i,  j)  =  Min 

'  3  K 

D(i,  K). 

The  classification  is  made  permanent  on  the  basis  of 
the  two  zone  parameters ,  \  1  and  X  2 ,  according  to  the 
following  rule: 

(1)  If  D(i,  j)  rS  \r  then  make  the  tentative 
classification  permanent  and  modify  the  class  descriptors. 

(2)  If  D(i,  j)  —  X2*  then  make  the  tenta¬ 

tive  classification  permanent  but  do  not  modify  the  class 
descriptors. 

(3)  If  X2  ^  D(i,  j),  then  delete  the  tentative 
classification,  creat~  a  new  class,  and  assign  as  the 
phonet  descriptor  of  the  class. 

These  zone  parameters  X 1  and  X  2  are  called  the  inner 
and  outer  zone  parameters,  respectively.  They  control  the 
classification  process  and  are  related  to  detection  perform¬ 
ance. 

The  zone  parameters  can  be  set  in  a  number  of  ways. 

They  could  be  set  arbitrarily  and  varied  until  optimum 
values  are  found  experimentally.  Or  they  could  be  related 

( 
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i  probability  of  misclassification  and  noise  energy 
i  speech  signal. 

Phe  other  two  class  descriptors,  the  phonet  and 
ice  sectors,  can  be  calculated  iteratively.  Suppose 
lence  of  observations  *  Let  X^  *  Xj  denote  the 

)lication  of  these  two  vectors  component-by-component; 
is,  let: 


xt  •  -  <Xil  Xjl-  Xi2  Xj2 . XiN  V  (7> 

»t  Xi/Xj  denote  the  division  of  these  two  vectors 
lent-by-component.  That  is,  let: 

V*J  "  <Xil/Xjl'  Xi2/Xj2 . Xit/XjN>  ,8) 

jiven  the  first  t  observations  of  the  sequence  ^XiJ  , 

Lrst  noncentral  moment  is: 

1 

Bt  *  t  I 

i=l 


(9) 


oond  noncentral  moment  is: 

St  *£  *i  •  Si' 

i=l 


(10) 


ie  second  central  moment  is: 

-t  =  t  -t  ~  (-t  *  -t5'  (11) 

ssume  another  observation  is  to  alter  these  class 


Lptors.  Then: 

t 


-t+1  t+1  -t  +  t+1  -t+1 


-t+1 


Ct  ♦  (Xt+1  •  W 


(12) 

(13) 
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(14) 


-t+1  t+1  ^t+l 


Bt+i *  • 


These  aquations  provide  the  means  for  rapid,  iterative 
calculation  of  the  class  descriptors  using  the  Array 
Processor. 

The  algorithm  given  by  Pal,  et.  al.  (Ref  6),  did  not 
specify  how  to  partition  the  feature  space  into  classes. 

It  assumed  class  partitions  and  their  descriptors  as  initial 
conditions.  That  is,  an  initial  set  of  class  descriptors 
was  resident  in  memory  at  the  start  of  their  algorithm. 

The  algorithm  given  below  in  step-by-step  format  does  not 
assume  that  the  sample  space  is  partitioned  initially. 
Rather,  it  partitions  feature  space  into  classes  whose 
boundaries  are  determined  by  noise  in  the  system,  the 
sufficient  statistic,  and  class  membership  variability. 


STEP 


ACTION 

Choose  a  variance  vector,  and  zone  param¬ 
eters  \  1  and  \  2*  0ne  way  to  do  this  i®  to 
let  V„  be  the  noise  variance  and  v  *  V„. 

Then  set  A2  *"  by  some  multiple.  Assume 
zero-mean  noise. 

Read  the  first  observation,  X^. 

Assign  to  class  A^. 

Describe  A^  by: 


1*  R«  m  X«»  C. 


Vo  +  <*1 


*!>•  Vj 


Read  the  second  observation,  X2* 

Compute  the  supervisory  parameter  between 
X2  and  R^»  that  is,  compute  D(2,l). 
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ACTION 

It  D<2,1)  2  \  x 
Then:  assign  x2  to 

modify  description  of 
else:  assign  X2  to  a2 
describe  A2 

Read  the  third  observation,  X3 • 
Choose  M£n  0(3,  K)  for  K  -  1,  2 

Say  the  minimum  is  at  K  *  1. 

If  D(3,  1)  2  \1 

Then:  assign  X3  to  class  A1 

modify  description  of  A1 
else:  assign  x3  to  class  a3 
describe  A3 


Read  mth  observation  X  . 

— m 

Choose  M£n  D  (M,  K)  for  K  *  1,  2, 

Say  the  minimum  is  at  K  *  1. 

If  D(M,  l)f\1 

Then:  assign  X  to  class  A, 

modify  description  of  A^ 
else:  assign  to  class  AR+1 
describe  class  A  , , 


Continue  in  thia  way  until  Mu  classes  have  been  des¬ 
cribed.  Mu  ia  the  largest  number  of  classes  tolerable  and 
is  the  smallest  number  of  classes  tolerable. 

At  this  point.  Mu  classes  have  been  described.  Bach 
class  was  created  because  the  observation  currently  under 
consideration  was  not  described  adequately  by  any  of  the 
classes  then  in  existence.  Adequacy  was  defined  in  terms  of 
class  variance.  Now.  at  most,  Mu  classes  are  to  be  permitted 
and  one  cannot  be  sure  at  this  point  that  future  observa¬ 
tions  will  be  described  adequately  by  the  current  set  of  Mu 
classes.  If  a  future  observation  is  not  described  adequately, 
a  new  class  will  need  to  be  created.  Since  processing  limi¬ 
tations  will  intrinsically  set  some  upper  limit  on  Mu>  the  number 
of  classes  must  be  reduced  to  M^,  some  number  less  than  Mu. 

This  reduction  is  to  be  accomplished  in  such  a  way  as  to 
increase  class  variance  for  some  of  the  classes  retained  so 
that  each  class  describes  adequately  more  of  the  feature 
space.  This  is  done  by  combining  the  most  similar  classes 
in  pairs;  that  is,  by  treating  the  phonets  as  observations 
and  computing  the  distance  measure  on  the  phonet  set. 

STEP  ACTION 


M+l  Read  the  first  phonet  R^. 

Compute  the  distance  between  R^  and  every 


ACTIOH 

other  phonet,  that  ia,  D  (1.  K)  for  K  »  2, 

...»  K,.  Chooaa  Min  D(l»  X). 

^  Ml 

Say  the  minimum  ia  at  K  ■  1. 

If  D(l,  1)  S  \ 

Thant  modify  deacription  of  claaa  A^  uaing 
aa  an  ohaervation 

elae:  do  not  modify  A^  deacriptora 
Delete  claaa  A^. 

Read  next  phonet.  Rj. 

Chooaa  Min  D(2»  K)  for  K  *  3»  ....  M 
M2 

Say  minimum  at  K  *  1. 

If  D<2.  l)i\j 

Then:  modify  A^  deacriptora  uaing  R2 
elae:  do  not  modify  A^. 

Delete  claaa  A2. 


Read  phonet  R^. 

Chooaa  Min  D(r,  K)  for  K  ■  r+1.  ...»  M 
Mr 

Say  miminum  at  K  ■  1 
If  D(r,  1)  2  Xx 

Then:  modify  deacription  of  A2  uaing 
elae:  do  not  modify  A^ 

Delete  Af. 


Continue  on  in  this  way  until  the  number  of  classes 
is  reduced  to  Renumber  the  classes  through  A^. 

Once  the  phonet  space  is  reduced  to  classes,  then 
classification  proceeds  using  the  rule  D(i,  j)  as  was  done 
above  when  the  phonet  vectors  were  classified.  How,  how¬ 
ever,  the  outer  zone  parameter  is  applied  to  the  detection 
problem  to  determine  if  a  new  class  should  be  added  to  the 
..t  [Aj]  . 

STEP  ACTION 


t  Read  observation  xn* 

Choose  Min  D(n,  K)  K  =  1,  ...,  m**=m 
K 

Say  minimum  at  K  =  1. 

If  D(n,  1)  «*  \x 

Then:  assign  to  class  Aj^ 

modify  class  descriptors  for  A^ 
else:  If  D(n,  1 )  —  X-2 

then:  assign  Xn  to  class 

do  not  modify  descriptors 
else:  assign  Xn  to  ^ 

describe  class  AM+1 


Continue  in  thi»  way  until  M+l  *  Mu»  When  M+l  =  Mu# 
than  reduce  tha  phonat  apaca  until  there  are  classes. 
Then  classify  more  observations... 


Det 

Supervisor 

\  r  a  i 

—  ■  > 

FIGURE  70.  Model  of  the  Supervisory  Detection  Process. 


A  model  of  the  supervisory  detection  process  is  shown 
in  Figure  70.  The  additive  noise  vector  is  assumed  to  have 
uncorrelated  components,  each  distributed  with  zero-mean 
and  variance  Vi.  Thus  X ^  +  1^,  its  mean  is  E^J  = 

X.  =  S^,  and  its  variance  is: 


/ 


V 


\ 


(15) 


m 


J 


Wozencraft  and  Jacobs  have  shown  that  the  optimum 


decision  rule  is  to  choose  tha  class,  a^,  to  which  the 
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observation.  Xj/  most  probably  belongs  (Ref  14:212-214). 
That  is.  choose 

Aj:  Max  pjs.j  PN  (X|S.)  (16) 

where  the  set  of  a  priori  probabilities  are  ^p|s^||  and 
PN  (X|S^)  is  assumed  to  be  the  joint  -  M  guassian  density 
of  X  conditioned  on  S^.  Melsa  and  Cohn  give  the  density 
with  mean  and  variance  V  (Ref  15:69).  and  develop  the 
multiple  decision  problem  (Ref  15:113-116).  We  parallel 
their  development. 

The  joint  -  M  guassian  density  is: 

PN  <25  |S± )  =  exp  {-4  (X  -  si)T  v-1  (x  -  s.)j  (17) 


V 


M 

7 T  2  7Tv. 

j=l  3 


The  optimum  decision  rule,  equation  16,  is  met  when  the 
logarithm  of  its  terms  is  maximum,  so  we  express  our  deci¬ 
sion  rule  as: 

A j :  Max  U.  -  h  (X  -  Sl)T  V  ” 1  (X  -  Sj)  (18) 

where:  M 

UA  =*  In  p{s.J  -  H  Y.  ln  (27Tvi)  (19) 

i=l 


Assuming  equally  likely  a  priori  probabilities, 
V-  i,  j  and  our  decision  rule  reduces  to: 


A.:  Min  (X.  -  S.  )T  V-1  (X  -  S 

j  i  _1  -1 

We  define  the  sufficient  statistic: 


then 

(20) 
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(21) 


i 

i 

i 


i 

it. 


i  * 


! 

* 

l 


* 


Li 


( 


D(i,  j)  ■  (XL  -  Sj)7  v-1  {XL  -  Sj) 

and  writ*  the  decision  rule  as: 

A,:  Min  D(i,  K).  (22) 

3  K 

Our  algorithm  approximates  as  the  ensemble  average 
of  all  observations  which  it  is  sufficiently  certain 
belongs  to  the  class  A^,  and  the  variance  as  the  ensemble 
variance  of  those  observations. 

The  probability  of  error  is  the  probability  that  there 
exists  at  least  one  class.  ,  given  S^,  for  which  D(i,  K) 
D(i,  j).  That  is.  we  let  denote  the  event: 

E^  «  {w:  D(i,  K)  •<  D(i,  j)  | (23) 

where  the  statistics  are.  themselves,  random  variables  and 
their  dependence  on  the  sample  point,  w.  is  suppressed  for 
convenience. 

We  apply  the  union  bound  as  given  by  Melsa  and  Cohn 
(Ref  15:115)  to  upper  bound  the  probability  of  error.  Say 
there  are  L  classes.  Then  we  say: 

«,{E|Sj}«I  p{ekISj}'  (24) 

K»1 

K*j 

The  task  is  to  find  P ^ |s J  V  K  /  j,  KiL. 

The  event  E^  occurs  whenever,  given  S^, 

D(i,  K)«  D(i,  j)  (25) 

(XA  -  SK)T  V'1  ~  SK)  -*  (X±  -  Sj)TV_1(Xi-Sj)  (26 
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SiT  V_1  xi  -  ^  V'1  *i  -  XiT  V-1  s*  ♦  s*  V-1  SK 
■*  xj  v"1  Xi  -  s*  v'1  XA  -  X?  V-1  s .  +  s*  v"*1  sj 
sj  V-1  SK  -  s]  V-1  s.^  2*1  V"1  <SK  -  S,) 


2xT  v”1  (sk  -  v’1  s£  -  Sj  v'1  Sj. 


(27) 


Using  Xj  «  S j  +  Nj<  we  have: 


2(Sj  +  Nj )TV~1  (SK  -  Sj)»sJv-1SK  -  SjV”1^ 


V  S  . 


„T  -1 


2Sjv'i  (h  -  s -  ^v*1  <4  -  Sjt  -  s jV~2Sj 

2NjV_1  (SK  -  Sj)»  S^V” XSK  -  2S  jV_1SK  +  sTv_1Sj  (28) 


Denote  the  left  aide  of  equation  28  by  the  random  variable 


Yv  and  the  right  side  by  the  nonrandom  variable  c 

IN  i\ 

mean  of  YKis  zero  and  its  variance  is : 


Var{vK}  -  e[y2] 


The 

(29) 

(30) 


2  E[(l  "ii  (SiK  “  Sjj))  ] 
i=l  Vi 

2  z[Z  nii(Sik~Sii)  (f  n11(slK~S11})]  (31) 
v,  1=1  V, 


i=l  i  1 

Because  the  noise  components  are  uncorrelated,  equation  31 
reduces  to: 


Var  |yk|  -  2*£  (S 
i=l 


iK  *  Si3> 


(32) 


So  we  write  equation  28  in  terms  of  the  zero  mean,  unit 


variance  random  variable: 


atS.K  -  Sjf 


ZK1 


(33) 
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where  |Sg  -  S ^ |  is  the  Euclidean  norm  between  the  phonets 
of  the  pair  of  claaaea  Ag  and  A  ^ ,  the  phonet  being  approxi¬ 
mated  by  the  mean  deacriptor.  The  random  variable  ^  ia 
distributed  guaasian  with  zero-mean  and  unit  variance.  So 
ita  density  is: 

pa  iv  *  -/Tj r  *xp  {~H  ’k] 

Then: 

i  oo  r  - 

PK  |Sj]  -  VTF /  l-*  *K  ]  d*K 


(34) 


(35) 


K1 


=  0(2^) 


(36) 


The  error  probability  conditioned  on  S..  is,  from 


equation  24, 

p{e|s3}«  X!  oiz^) 


(37) 


K=1 

K*j 


To  bound  the  overall  average  error  probability,  equation 
37  ia  averaged  over  all  classes,  so: 

pH-[iE  <3s> 

1*1  K=1 

This  average  can  be  further  upper-bounded  in  terms  of  the 
minimum  Z^.  That  is,  denote: 


Z  .  -  Min 

min 


in  I  Min  z,,.  j . 
1  [k/1 


Then: 


P  {e}  -  (L  -  1)  0  (Zmin). 


(39) 


(40) 
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The  Q-function  can  be  approximated  (Ref  15:258)  by 


Q(X)- 


2  exp  ( -H 


{e}«  l  (L  -  11  (-  *»  2min} 


2  V  (z„ 


zmin  +  yin  » 


One  might  relate  the  probability  of  classification 

error  to  the  maximum  number  of  classes,  My,  and  the  outer 

supervisory  parameter  in  the  following  way.  One  could 

consider  the  minimum  value  of  ZR,  the  right  side  of  equation 

28,  to  be  the  signal-to-noise  ratio  in  detecting  the  closest 

two  classes.  The  signal-to-noise  ratio  at  the  input  to  the 

detector  could  be  measured.  It  would  be: 

SNR  =  i*i  I*1  ~  |2  (43) 

Var  {N  j 

For  a  specified  p|e|  in  the  detector,  then, 

Zmin  =  Q-1  (  PM_)  (44) 

\(Mu -1)J 

from  equation  40  can  be  thought  of  as  the  minimum  distance 
normalized  to  the  input  variance'  at  the  detector.  Scaling 
this  distance  by  the  SNR  at  the  detector  input,  one  gets 
the  outer  supervisor  parameter  from: 


\2  Z  . 
Au  mu 


SNR. 


Once  the  feature  space  is  partitioned  into  classes  of 
phonetic  units,  then  those  phonetic  units  can  be  recognized 


in  speech.  Our  Acoustic  Analyser*  Seelandt's  Speech  Sound 
Analysis  Machine  (Ref  1)  or  some  other  distance  computer* 
can  toe  used  to  accomplish  the  recognition.  The  constraint, 
of  course,  is  that  the  phonetic  units  used  to  partition 
feature  space  must  toe  the  same  as  those  generated  from 
speech  by  the  distance  computer. 

In  this  chapter,  we  have  presented  the  results  of  this 
project.  We  presented  the  Acoustic  Analyzer  and  illustrated 
the  behavior  of  observations  under  a  minimum  distance  deci¬ 
sion  criterion.  Then  Conclusions  and  Recommendations  are 
presented  in  Chapters  IV  and  V,  respectively. 
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IV.  CONCLUSIONS 


We  have  implemented  Seelandt's  Speech  Sound  Analysis 
Machine  (SSAM)  on  the  Eclipse  S250  Array  Processor.  We 
were  able  to  generate  phonetic  units,  which  we  called 
observations  and  phonets,  from  speech  files.  We  calculated 
distances  between  observations  and  phonets,  and  we  were 
able  to  output  the  best  phonet  choices  along  with  other 
parameters  for  input  to  Montgomery's  word  recognizer  (Ref 
10).  We  also  provided  a  graphics  capability  for  plotting 
phonetic  units,  distances,  and  time  files. 

We  also  reviewed  the  parameters  Seelandt  identified  for 
further  study  (Ref  1)  and  provided  a  basis  for  choosing  param¬ 
eter  values.  With  the  values  we  chose,  we  observed  that  adja¬ 
cent  phonetic  units  in  voiced  vowel  sounds  resemble  each  other 
in  both  distance  rules.  We  found  that  while  adjacent  phonetic 
units  in  speaker-absent  background  noise  and  in  "r"  sounds  did 
not  resemble  each  other,  they  did  in  "s”  sounds.  We  conclude 
that  there  seems  to  be  a  sufficient  tendency  for  phonetic 
units  to  cluster  to  warrant  further  work  along  these  lines. 

We  also  adapted  an  algorithm  for  self-supervised  vowel 
recognition,  which  used  the  first  three  formants  as  features 
in  the  absence  of  noise  to  our  problem  of  phonetic  unit  recog¬ 
nition.  The  algorithm  is  explained  and  we  showed  how  the 
S250  Array  Processor  and  our  Acoustic  Analyzer  could  be  used 
to  implement  it. 


RECOMMENDATIONS 


We  were  unable  to  accomplish  much  important  work  that 
s  to  be  done  before  a  phonet  detector  built  along  the 
5  of  our  Acoustic  Analyzer  can  be  optimized.  One  major 
rt  is  needed  in  implementing  an  algorithm  for  parti tion- 
the  feature  space.  The  algorithm  by  Pal,  et.  al.  (Ref 
that  we  considered  in  this  project  is  one  of  many 
ads.  While  we  feel  it  has  merit:  we  found  many  others 
also  warrant  careful  consideration.  The  recent  book 
evijver  and  Kittler  (Ref  9)  discusses  the  topic  in 
il.  The  text  by  Fukunaga  (Ref  13)  also  treats  the 
ect  and  presents  algorithms  which  could  be  applied  to 
problem.  An  algorithm  should  be  implemented. 

Once  a  method  for  partitioning  feature  space  is  imple- 
ed,  the  parameters  Seelandt  discussed  and  that  were 
ted  in  this  project,  should  be  investigated  experiment- 
.  That  is,  their  effect  on  recognition  of  phonetic 
s  chosen  by  a  feature  space  partitioning  algorithm  should 
etermined. 

Several  features  need  to  be  added  to  our  Acoustic 
yzer.  For  one,  the  capability  to  overlap  FFT  windows 
izes  other  than  just  128  points  needs  to  be  incorporated, 
routine  S64,  which  allows  the  necessary  window  sizes  but 
not  permit  window  overlapping,  can  serve  as  a  departure 


point  for  a  general  purpose  FFT  routine  that  incorporates 
the  overlap  feature,  additional  window  shapes,  and  more 
sophisticated  spectrum  shaping  options. 

Also  needed  is  the  capability  for  the  distance  routines, 
DSTA  and  DSTN,  to  do  a  cross-correlation  over  shifts  of  only 
a  few  components  between  observation  and  phonet  spectrum. 

In  this  way,  the  number  of  phonet  clusters  in  feature  space 
may  be  reduced  if  several  phonet  clusters  are  only  shifted 
versions  of  others.  Dete  ♦■ion  could  proceed  conditionally 
on  frequency  shift,  and  then  choices  between  the  shifts 
could  be  made.  But  before  this  is  implemented,  it  may  be 
valuable  to  implement  first  a  partitioning  algorithm.  Once 
the  feature  space  is  partitioned,  the  cluster  characteristics 
may  suggest  whether  spectral  shifts  may  offer  improvement, 
how  far  to  shift,  and  how-  much  improvement.  Answers  to  these 
questions  depend  heavily  on  the  way  that  feature  space  is 
partitioned. 


Ill 


! 


BIBLIOGRAPHY 


1.  Seelandt,  Karl  G.  Computer  Analysis  and  Recognition 

of  Phoneme  Sounds  in  Connected  Speech,  MS  Thesis  GE/EE/81D-53. 
School  of  Engineering,  Air  Force  Institute  of  Technology, 
Wright- Patterson  APB,  Ohio,  December  1981. 

2.  Felkey,  Mark  A.  Automatic  Recognition  of  Phonemes  in 
Continuous  Speech,  MS  Thesis  GE/EE/80D-20,  School  of 
engineering.  Air  Force  Institute  of  Technology,  Wright- 
Patterson  AFB,  Ohio,  December  1980. 

3.  De  Souza,  Peter  and  Thomson,  Peter  J.  "LPC  Distance 
Measures  and  Statistical  Tests  with  Particular  Reference 
to  the  Likelihood  Ratio,"  IEEE  Transactions  on  Acoustics, 
Speech,  and  Signal  Processing,  Vol  ASSP-30,  No.  2,  April 
1982. 

4.  Kassam,  Salem  A.  A  Bibliography  on  Nonparametric 
Detection,  AD  A080187,  prepared  for  the  Air  Force  Office 
of  Scientific  Research,  AFOSR-TR-80-0025,  1979. 

5.  La.niotis,  D.  G.  Adoptive  Systems,  AD  776015,  prepared 
for  the  Air  Force  Office  of  Scientific  Research,  APOSR— 

'’4-0309,  February  1974. 

6.  Pal,  S.  K.,  Datta,  A.  K.,  and  Dutta  Majumder,  D. 

“A  Self-Supervised  Vowel  Recognition  System, "  Pattern 
Recognition,  Pergamon  Press  Ltd,  Great  Britain,  Vol  12, 

27-34. 

7.  Nonparametric/Robust  Methods  and  Non-Guassian  Models  in 
Communication  Theory ,  Johns  Hopkins  University,  AD-763407, 
prepared  for  the  Office  of  Naval  Research,  1971. 

8.  Rabiner,  L.  R.  and  Schafer,  R.  W.  Digital  Processing  of 
Speech  Signals,  Prentice  Hall,  1978. 

9.  Devijver,  R.  P.  and  Kittler,  J.  Pattern  Recognition; 

A  Statistical  Approach,  Prentice  Hall,  1982. 

10.  Montgomery,  Gerald.  Isolated  Word  Recognition  Using 
Fuzzy  Set  Theory,  MS  Thesis  GE/EE/82D-74,  School  of  Engineer 
ing,  Air  Force  Institute  of  Technology,  Wright-Patterson 
AFB,  Ohio,  December  1982. 


11.  Sivian,  L.  J.  "Speech  and  its  Power  Measurement, "  the 
Bell  System  Technical  Journal,  Vol  8,  pp.  646-661,  1929. 


I 


12.  Kuligowski,  T.  J.  An  Automatic  Speech  Teaching 
Responder.  MS  Diesis,  United  States  Naval  Postgraduate 
School,  AD  851813,  June,  1968. 

13.  Fukunaga,  K.  Introduction  to  Statistical  Pattern 
Recognition,  Academic  Press  Inc.,  1972. 

14.  Wo2encraft,  J.  M. ,  and  Jacobs,  I.  M.  Principles  of 
Communication  Engineering,  John  Wiley  and  Sons,  New  York, 
1965. 

15.  Melsa,  J.  L. ,  and  Cohn,  D.  L.  Decision  and  Estimation 
Theory,  McGraw-Hill,  Inc.,  1978. 


113 


APPENDIX-  A 

SPECTROGRAM  OP  PILE  CT56.SP 
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APPENDIX  B 

SOFTWARE  DOCUMENTATION 


STRUCTURE  OF  THE  ACOUSTIC  ANALYZER; 

The  Acoustic  Analyzer  consists  of  one  executive  routine, 
15  subroutines,  and  calls  to  two  library  files,  ibe 
executive  routine,  DRVR.SV,  and  subroutines  CHKO,  FCTR, 

GFDB,  along  with  subroutines  from  the  library  files,  are 
loaded  into  one  save  file,  DRVR.SV.  An  overlay  node  is 
created  in  the  save  file  into  which  segments  of  the  overlay 
file,  DRVR.OL,  are  loaded  when  called.  The  macro  file  which 
calls  the  loader  utility  is  included  in  this  documentation. 

Figure  7l  illustrates  the  software  structure.  As 
called,  routines  DRSQ,  DSTN,  DSTA,  or  PLTO  are  loaded  into 
the  executive  routine  at  the  overlay  node,  along  with  their 
supporting  subroutines.  The  other  routines:  CHKO,  FCTR, 

GFDB,  and  the  libraries  are  accessible  to  the  contents  of 
the  overlay  node. 

Below  is  a  list  of  each  routine  in  the  Acoustic  Analyzer, 
along  with  a  brief  statement  of  its  purpose: 

(1)  DRVR  is  the  executive.  It  allows  the  operator 
to  select  either  the  spectral  option,  either  of  the  two 
distance  options,  or  the  plot  options. 

(2)  DRSQ  calls  either  S128  or  S64  to  compute  the 
Fast  Fourier  Transform  ( FFT)  of  an  input  file. 

(3)  S128  computes  FFTs  of  an  input  file  using  a 
128  point  window  which  overlaps  adjacent  windows  by  64  points. 
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FIGURE  71.  Structure  of  the  Acoustic  Analyzer.  Options 

available  from  DRVR  are  program  segments  resi¬ 
dent  on  disk  and  loaded  into  the  overlay  node 
on  call. 


(4)  S64  computes  FFTs  of  an  input  file  using  a 
nonoverlapped  window  of  any  length  which  is  an  integral 
multiple  of  two  in  the  range  16-1024  points,  inclusive. 

(5)  DSTA  computes  distances  between  observations 
and  phonets,  and  stores  the  symmetric  distance  matrix  on 
disk. 

(6)  DSTN  computes  the  distance  matrix  that  DSTA 
computes,  but  chooses  a  selected  number  of  best  template 
matches  for  output  to  disk. 

(7)  CHOOS  is  the  subroutine  that  DSTN  calls  to 
find  best  matches  between  observations  and  phonets. 
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(8)  CHKJ  is  called  by  DSTN  to  output  data  when 
distance-file  elements  are  accumulated. 

(9)  PLTO  allows  the  operator  to  call  options 
h  either  display  spectrum,  display  distances  output 
STA,  display  distances  output  by  DSTN,  or  display  the 
ents  of  integer  disk  files. 

(10)  PLTA  displays  distance  files  computed  by 

• 

(11)  PLTS  displays  spectral  files  computed  by 

(12)  PLTN  displays  distance  files  computed  by 

I. 

(13)  PLTT  displays  any  segment  of  512  integers 
ewer  from  a  disk  file. 

(14)  CHKO  is  called  by  DSTA  and  DSTN  to  type  a 
.ion  statement  to  screen  if  a  file  unit  number  is  in  use. 

(15)  FCTR  is  called  by  GFDB  to  factor  an  integer 
i  a  linear  combination  of  four  terms. 

(16)  GFDB  is  called  by  PLTA,  PLTS,  PLTN,  and  PLTT 
:ompute  the  location  of  a  specified  data  point  in  a  file. 
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PILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OF  LAST  REVISION: 


LDDRVR.MC 

Command  Line  Interpreter 
September  21,  1982 
D.  Martin 
Acoustic  Analysis 
LDDRVR 

September  21,  1982 


PURPOSE: 

This  macro  file  loads  program  DRVR.SV.  An  overlay 
node  is  established  and  an  overlay  file,  DRVR.OL,  and  load 
map,  DRVR.LM,  are  created. 


DESCRIPTION: 


Location :  DP4 : BRATCHET 


Size:  LDDRVR.MC 

DRVR.SV 
DRVR. FR 
DRVR.  RB 
DRVR.OL 
DRVR.LM 


155  bytes 
51200  bytes 
2613  bytes 
3054  bytes 
40960  bytes 
9491  bytes 


PROGRAM  USE: 

This  macro  uses  the  libraries  GRPH.LB  and  APS. LB;  a 
link  to  these  are  required.  It  loads  the  main  program 
DRVR  and  the  subroutines  DRSQ,  S128,  S64,  DSTA,  DSTN, 
PLTA,  PLTO,  PLTS ,  PLTT,  PLTN,  CHKO,  CHKJ,  CHOOS,  FCTR, 
and  GFDB. 
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DELETE/ V  DRVR.LH 

RLDR/C/R/E/P  2000/N  DRVR  CDRSQ  S124  864 ,DSTA,DSTN  CHDOS  CHIJ, 
PLTO  PLTA  PITS  PLTN  PLTT1  CHIO  FCTR  CFDB  1 
DRVR.LH/L  GRPH.LB  APS. LI  IFLIBI 


FILE:  DRVR 

LANGUAGE:  FORTRAN  5 

DATE:  September  21,  1982 

AUTHOR:  D.  Martin 

SUBJECT:  Acoustic  Analysis 

CALLING  SEQUENCE:  DRVR 

DATE  OF  LAST  REVISION:  September  21,  1982 


PURPOSE: 

This  program  drives  the  acoustic  analysis  package 
developed  by  Captain  Dan  Martin  for  his  thesis,  when  this 
program  is  executed,  a  menu  is  displayed  from  which  the 
operator  chooses  an  analysis  option:  to  get  spectrum,  to 
get  distances,  to  get  best  distances,  or  to  plot  results. 


DESCRIPTION: 

Location:  DP4 : BRATCHET 

Size:  DRVR. S V 

DRVR. FR 
DRVR. RB 
DRVR. OL 
DRVR. LM 


PROGRAM  USE: 

CAUTION.  The  Array  Processor  must  be  properly 
initialized  prior  to  any  attempt  to  run  this  routine.  A 
macro  file,  LDAPS. MC,  in  DP4 : BRATCHET  will  accomplish  this 
initialization.  Also,  at  least  14  blocks  of  extended 
memory  must  be  allocated  to  the  ground  on  which  this 
routine  is  run.  Further,  DRVR.OL,  GRPH.LB,  and  APS. LB 


51200  bytes 
2613  bytes 
3054  bytes 
40960  bytes 
9491  bytes 


need  be  accessible  to  the  directory  in  which  this  routine 
is  run. 

This  program  is  not  called  by  any  routine  and  calls 
routines  DRSQ,  DSTA,  DSTN,  and  PLTO.  When  executed,  this 
program  presents  a  menu  from  which  the  operator  selects 
one  of  the  following  options: 

(1)  To  compute  observation  spectrum  from  a 
speech  file. 

(2)  To  compute  distances  between  observations  and 
phonets.  Phonets  are  observations  set  aside  to  represent 
recognizable  units  of  speech. 

(3)  To  compute  distances  between  observations  and 
phonets  as  in  the  previous  option,  but  to  choose  a  selected 
number  of  best  matches  based  on  those  distances. 

(4)  To  plot  files  on  the  Tektronix  4010-1  Graphics 

Terminal. 

PURPOSE 

Number  of  available  ex¬ 
tended  memory  blocks. 

Number  of  extended  memory 
blocks  over  and  above  those 
used  in  the  logical  address 
space  window. 

Sampling  rate  for  DRSQ. 

Vector  in  array  processor 
memory. 

Switch  to  select  options. 


LIST  OF  VARIABLES: 
VARIABLE  TOPE 

XMEM  Integer 

INITOEM  Integer 

SFREQ  Real 

DTARAY  Real  Array 

SKIP  Integer 


%■- — _•  •- 
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SWITCH  SETTINGS: 

SWITCH  SETTING 

SKIP  =  1  to  compute  observations 

=  2  to  compute  distances 

=  3  to  select  best  matches 

=  4  to  plot  files 

=  5  to  stop 

The  remainder  of  this  section  is  devoted  to  a  discussion 
of:  (1)  our  use  of  the  S250  Array  Processor,  (2)  extended 

memory  management,  and  (3)  modification  of  the  Acoustic 
Analyzer  (AA). 

The  Eclipse  S250  Array  Processor  combines  an  Eclipse 
S250  CPU  with  a  floating-point  Array  Processor.  This 
allows  high-speed  computation  on  real  and  complex  vectors. 

We  used  FORTRAN  5  as  the  application  programming  language 
and  Data  General's  Array  Processor  Software  (APS)  running 
under  Data  General's  mapped  RDOS  operating  system.  The 
AP/S250  contains  both  an  independent  pipelined  floating-point 
multiplier  and  add/subtract  compare  units  that  operate 
simultaneously.  Several  multiply  or  add/subtract  opera¬ 
tions  overlap  during  given  time  periods. 

AP  memory  consists  of  an  AP-mounted,  64  bit  by  IK 
(8kByte),  20MSEC  bipolar  RAM.  This  memory  is  dual-ported 
which  allowed  us  to  address  it  by  two  methods: 

(1)  We  addressed  APM  as  part  of  main  memory 
in  which  case  we  mapped  APM  into  our  32K  word  logical 
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address  space.  In  this  way,  data  transfer  between  APM  and 
logical  address  space  where  unnecessary. 

(2)  We  addressed  APM  as  local  memory  by  using  AP 
instructions. 

We  used  the  two  categories  of  APS  routines:  support 
routines  and  interface  routines.  With  the  support  routines, 
we  handled  the  following  tasks: 

(1)  AP  initialization. 

(2)  Setting  up  control  blocks  and  default  values. 

(3)  AP  memory  management. 

(4)  Miscellaneous  utility  tasks. 

With  the  interface  routines,  we  performed  specific  Array 
Processing  operations,  such  as  to  multiply  two  arrays,  for 
example.  We  used  first  level  interface  routines,  referred 
to  as  V- routines.  Before  executing  a  V- routine,  we  defined 
all  data  and  control  words  required  by  the  APS  for  that  v- 
routine.  This  setup  is  accomplished  by  calling  support 
routines  as  specified  in  the  V-routine  description.  Besides 
control  block  words,  ve  used  a  number  of  special  APS  param¬ 
eters  to  make  our  application  program  independent  of  actual 
main  memory  addresses  or  APM  offsets. 

We  managed  memory  in  our  Acoustic  Analyzer  (AA)  using 
memory  mapping  and  static  memory  allocation.  We  chose  to 
declare  a  nine-memory  block  window  in  logical  address  space 
into  which  we  would  map  APM  and  other  arrays  which  we  used 
as  ports  to  extended  memory.  These  port  arrays:  IDSP, 
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I DOB  in  DRSQ,  and  DIST,  FON,  and  OBS  in  DSTA  and  DSTN,  were 
mapped  through  extended  memory  to  access  data  with  a  mini¬ 
mum  of  disk  accesses.  This  scheme  is  illustrated  in  Figure 
7  for  DSTA  and  DSTN.  In  extended  memory,  the  first  four 
memory  blocks  are  dedicated  to  the  AP;  this  fact  cannot  be 
altered  by  the  application  program.  Our  window  in  logical 
address  space  consists  of  the  four  real  arrays:  'WKOBS,  WKFON, 

OBS,  and  PON;  and  the  one  integer  array  DIST.  Each  array  is 
1024  elements  long.  Array  Processor  memory  is  mapped  into 
WKOBS  and  WKFON.  Arrays  OBS,  FON,  and  DIST  are  mapped 
through  extended  memory  to  access  observation  spectrum, 
phonet  spectrum,  or  to  store  distances.  From  these  arrays, 
data  are  parcelled  to  WKOBS  and  WKFON  for  processing.  The 
scheme  for  DRSQ  is  the  same  except  that  the  entire  window 
is  not  used. 


FIGURE  7.  Memory  Map  for  DSTA  and  DSTN 
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We  have  modified  the  Acoustic  Analyzer  for  use  by  two 
students  also  doing  their  projects  in  the  Speech  Laboratory. 
We  allowed  them  to  run  the  spectral  option  and  the  distance 
option  non-interactively  by  setting  the  switches  in  copies 
of  the  object  code  to  their  specifications.  They  now  have 
access  to  the  original  interactive  package  with  graphics 
capability  as  well  as  a  tailored,  non-interactive  package 
capable  of  being  run  by  MACRO  files  under  RDOS.  They  use 
the  output  of  the  tailored  Acoustic  Analyzer  as  input  to 
Montgomery's  word  recognition  machine  (Ref  10)  which  is  also 
implemented  on  the  Eclipse  S/250  in  the  Speech  Laboratory. 

RELATED  PROGRAMS; 

PLTA,  PLTS,  PLTO.  PLTO,  PLTT,  DRSQ,  S128,  S64,  DSTA, 
DSTN,  CHKJ,  CHKO ,  CH00S,  F<  ""R,  GFDB,  GFPH2. 
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ROUTINE  DRVR.FR  NOTE*  THIS  ROUTINE  18  FOR  FORTRAN  8  !  ! 
THIS  ROUTINE  DRVES  DRSQ,S12B,S64,DSTA,DSTN,CH00S 
CHI J ,PLT0 ,PLTA ,PLTS , PITH ,PLTT ,CHK0 ,FCTR ,GFDB . 

BY i  CAPT  DAN  HARTIH 

BATE*  9/27/82 

SUBJ :  ACOUTSIC  ANALYSER 

THIS  ROUTINE  DRIVES  THE  ACOUSTIC  ANALYSIS  PACIA6E 
I  DEVELOPED.  THE  OPERATOR  IS  PRQNPED  TO  SPECIFY! 

1)  AHOUNT  OF  EXTENDED  HEHORY  TO  USE, 

2)  TASK  TO  DO 

•GET  SPECTRUR 
•GET  DISTANCES 
•GET  H-BEST  DISTANCES 
•GET  PLOTS 
•STOP 

FOR  DETAILS,  SEE  THE  USERS  NANUAL  OR  NY  THESIS. 


INCLUDE  “ARRATP:F5APS.FR“ 

EXTERNAL  ODRSO,ODSTA,ODSTN,OPLTO 

PARAMETER  LH*512 

PARAMETER  NDP«1024 

PARAMETER  HLM*LH/2 

REAL  DTARAY ,B,PREN,UIND,SFREO,IDOB 

INTEGER  FI LEI  ,IOPTN,IPRE,IDB,IHAM,FLEN,SKIP,INITNEM 

INTEGER  ISTART,ILAST, COUNT, FILE2, SMOOTH, HRTD, TEST, FLAG 

INTEGER  NBLKS,1FRQ( 10) ,IAMP( 10) ,XMEM,URTDA,CB1 , DUMMY 

INTEGER  CNTBR,1DSP ,SXHEM,GSPCT,IVAL 

COMMON  /APN/  DTARAY (NDP) ,B(NDP) 

COMMON  /  VALS  /  IDSP(NDP)  ,1D0B(NDP) , I H AN , I PRE , U IND ( LH ) 

COMMON  /  VALS  /  PREH(LM) ,1NITNEN,SXHEH, DUMMY (306) 

COMMON  /  VALT  /  SKIP, COUNT, SMOOTH, MRTD,NBLKS,XHEH,SFREO,FLEN 
COMMON  /  VALT  /  FILE1 (13) ,FILE2( 13) ,FLAG ,CNTBR,ISXHEH 

COMMON  /  VALT  /  JAOUT,LAStCD 

COMMON  /  VALV  /  CB1(0:CBHAX) 

COMMON  /  VALU  /  IVAK205B) 


30  CALL  VMEN(XMEM,IER) 
lNITMEH*XNEH-9 
XNEN'XNEH- 1 

35  TYPE“*YOU  HAVE*,XNEN,*  IK  BLOCKS  OF  EXT  MEMORY  ACCESSABLE’ 

40  ACCEPT“»*HOy  MANY  EXTENDED  MEMORY  IK  BLOCKS  HILL  YOU  USE?  ” 

IF(IXREH.LE«XHEH)COTO  45 

TYPE“YOU  DON'T  HAVE  THAT  MUCH  MEMORY  AVAILABLE." 

GOTO  35 

45  XHEH(IXNEH*10 

IF(XNEH.CT.O)GOTO  50 
lXNiNal*XNEN 

TYPE'* •••HEED  SPECIFY  AT  LEAST", IXNEM,"  MORE  BLOCKS!" 

GOTO  30 

50  XHEHiIXNEM*?  ;LOSE  9  BLKS  TO  THE  HINDOH 

ISXMEM*XMEM 
SFREO*IOOO.O 


IXNEM 
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CALL  RESET 

CALL  0VOPN(20,"DRVR.0L",IER> 

CALL  CHECK IER) 

INITIALIZE  ARRAY  PROCESSOR  AND  HEHORY  HAP 
CALL  APIHIT(NIL,DTARAY,T,INITHBB,IER) 

CALL  APHAP(ITARAY,0,4,IER> 

TYPB"**HHAT  HILL  YOU  HAVE  HE  DO?" 
TYPE"***»CHOOSE  OPTION****" 

TYPE"*" 

TYPE"!  TO  CET  SPECTRUH." 

TYPE"2  TO  CET  DISTANCES." 

TYPE’3  TO  CET  H-BEST  DISTANCES." 

TYPE*4  TO  CET  PLOTS. " 

TYPE"5  TO  STOP." 

ACCEPT"*»ENTER  OPTION  DESIRED)  ",SKP 

IP(SKIP.NE.l)COTO  110 
CALL  0VL0D(QDRSQ,-1,1ER) 

CALL  CHECKIER) 

CALL  DRSQ 
COTO  00 

1F(SK1P.NE.2)CQT0  120 
CALL  0VL0D(0DSTA/*1/1ER) 

CALL  CHECX(IER) 

CALL  DSTA 
GOTO  00 

1F(SKIP,HE.3)C0T0  130 
CALL  OVLODtODSTNf *lf 1ER) 

CALL  CHECKIER) 

CALL  DSTN 
COTO  00 

IF(SKIP.NE.4)C0T0  140 
CALL  0VL0D(0PLT0;*1 , IER) 

CALL  CHECKIER) 

CALL  PLTO 
COTO  00 

IF(SKIP.EQ.5)G0T0  1000 

TYPE"ERROR)  COULDN'T  FIND  YOUR  OPTION." 

TYPE" YOU  SELECTED  OPTION)  ",SKIP 

COTO  00 

CALL  RESET 

STOP 
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FILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OF  LAST  REVISION: 

PURPOSE: 

This  program  calls  subroutines  S128  and  S64  to  compute 
the  Fast  Fourier  Transform  ( FFT)  of  an  input  file.  It  uses 
the  Eclipse  AP/S250  Array  Processor.  Although  flexible, 
these  routines  were  designed  to  transform  speech  files  pro¬ 
duced  by  the  Cromemco  analog-to-digital  system. 

DESCRIPTION: 

Location:  DP4 : BRATCHET 

Size:  DRSQ.FR 

DRSQ.RB 

PROGRAM  USE: 

This  program  is  called  by  DRVR  and  calls  S128,  S64,  and 
CHKO.  This  is  a  flexible  routine  in  that  several  options 
are  available  with  regard  to  windowing  of  input  files, 
emphasis  and  scaling  of  output  files.  A  block  diagram  is 
included  at  Figure  72. 

When  this  routine  is  called,  it  prompts  the  operator 

i  for  the  following  values  and  options: 

(1)  To  pre-emphasize  high  frequencies  or  not? 


11400  bytes 
15840  bytes 


DRSQ 

FORTRAN  5 

September  21,  1982 
D.  Martin 
Acoustic  Analysis 
DRSQ 

September  21,  1982 


i 
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Enter 


A 


"Header  X  No 
v(  58)=1 


Calculate 
/indow  and 
emphasis 
vectors 


Fill  header 
for  spectrum 
file 


Pest  *  1 


Construct 
test  disk, 
file 


,/More  Yes 

^speech  filfP 
,X\to  do^  i 


Initialize 
APS  and  APM 


fReturr 


Read 

speech 

file 


FIGURE  72. 


Flowchart  of  DRSQ 
131 
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(2)  If  pre-emphasis  is  selected,  at  what  rate  in 
dB/Octave,  and  corner  frequency? 

(3)  To  de-emphasize  low  frequencies  or  not? 

(4)  If  de-emphasis  is  selected,  at  what  rate  in 
dB/Octave,  and  corner  frequency? 

(5)  What  FFT  window  length?  Any  multiple  of  two 
between  16  and  1024,  inclusive,  can  be  selected. 

(6)  Choose  Hamming  or  rectangular  window.  If  a 
rectangular  window  is  selected,  no  window  vector  is  com¬ 
puted.  If  a  Hamming  window  is  selected,  the  window  vector 
is  calculated  prior  to  the  S128  or  S64  call. 

(7)  If  a  window  length  of  128  points  is  selected, 
the  program  will  prompt  for  a  specification  to  overlap  or 
not.  If  overlap  is  selected,  the  window  will  be  moved 
along  the  input  file  in  64  point  increments. 

(8)  To  either  normalize  the  observation  energy 
to  preset  value  or  to  divide  each  spectral  component  by  the 
observation  energy.  The  observation  energy  is  the  sum  of 
the  squares  of  each  spectral  component:  one  through 
N/2-1. 

(9)  Specify  whether  to  write  spectrum  to  disk 
only,  to  the  terminal  as  well,  or  to  the  line  printer  as 
veil. 

(10)  Specify  if  a  test  disk  file  is  to  be  created 
and  processed,  or  a  speech  file  is  to  be  read  from  disk  and 
processed.  DRSQ  can  generate  a  disk  file  with  a  signal 
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i  is  a  sum  of  up  to  ten  tones  at  as  many  amplitudes, 
iperator  is  prompted  for  these  values  if  the  test  option 
tlected. 

(11)  Specify  speech  file  name. 

(12)  Specify  a  disk  file  name  in  which  to  store 
tbservations.  If  this  file  does  not  exist,  it  will  be 
:ed. 

(13)  Specify  the  portion  of  the  speech  file  to 
»ss.  Ibis  is  done  by  specifying  the  first  and  last 
blocks  to  process . 

A  header  is  prepared  and  written  to  disk  block  zero, 
e  1  contains  the  header  assignments  for  the  nontest 
on.  Table  2  for  the  test  option. 

Extended  memory  blocks  zero  through  eight  are  unavail- 
for  use;  they  are  used  by  the  window  declared  in 
sal  address  space.  The  speech  file  is  read  into  all 
Lable  extended  memory  blocks  and  array  IDSP  is  mapped 
jgh  extended  memory  accessing  data.  These  data  are 
elled  from  IDSP  to  Array  Processor  memory  in  FFT  time 
ss.  After  the  FFT  computation,  the  real  spectrum  is 
sf erred  through  real  array  IDOB  to  extended  memory, 
r  extended  memory  contents  have  been  replaced  with 
rvations,  its  contents  are  written  to  disk. 


HEADER  ELEMENT 

CONTENT 

1-13 

Observation  file  name. 

14-26 

Speech  file  name. 

27 

IPRE,  pre-emphasis  switch. 

28 

IDB,  pre-emphasis  rate  (dB/Octave). 

29 

FREQ,  frequency  at  which  pre¬ 
emphasis  starts. 

30 

FLEN,  FFT  window  size. 

31 

I HAM,  window  switch. 

32 

Unused  (set  to  zero). 

33 

TEST,  test  option  switch. 

34-54 

Unused  (set  to  zero). 

55 

COUNT,  number  of  first  time  slice 
to  do. 

56 

LTSTDO,  number  of  last  time  slice 
to  do. 

57 

HL,  number  of  elements,  per 
observation. 

58 

Overlapping  switch. 

59 

Number  of  disk  blocks  in  observa¬ 
tion  file. 

60 

IDPRE,  de-emphasis  switch. 

61 

DDP,  de-emphasis  rate  (dB/Octave). 

62 

DFREQ ,  frequency  at  which  de¬ 
emphasis  ends. 

63-256 

Unused  (set  to  zero). 

TABLE  1.  Header  Assignments  for  Spectral  File  Computed  by 
DRSQ  Calling  S128  or  S64. 
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Test  file  header  -  but  the  header  isn't  written  to  disk. 


HEADER  ELEMENT 


1-13 


14-26 


35-44 

45-54 


55-57 


59-256 


CONTENT 


Test  file  name. 

Observation  file  name. 

Switch:  1  =  pre-emp’nasize/O  = 

don't  pre-emphasize. 

Pre-emphasis  slope. 

Pre-emphasis  corner  frequency. 

Number  of  time  points-per-FFT. 

Switch:  1  =  Hamming  window/O  = 

rectangular  window. 

Unused. 

Switch:  1  =  create  test  file/O 

don't  create  test  file. 

Number  of  tones  in  test  file. 

Tone  frequencies  (HZ). 

Tone  amplitudes. 

Unused. 

Switch:  1  =  overlapping  time 

slice/O  =  nonoverlapping. 

Unused. 


TABLE  2.  Header  Assignments  for  Test  Option. 


**nnK*i*x**t? .k-,'  ti 
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LIST  OF  VARIABLES 


VARIABLE 

TOPE 

PURPOSE 

XMEM 

Integer 

Number  of  extended  memory 
blocks  to  use. 

IOPTION 

Integer 

Switch  to  set  parameter 
options. 

IPRE 

Integer 

Pre-emphasis  switch. 

IDB 

Integer 

dB  per  Octave  pre¬ 
emphasis. 

FREQ 

Integer 

Starting  pre-emphasis 
frequency. 

IDPRE 

Integer 

De-emphasis  switch. 

DDB 

Integer 

dB  per  Octave  de-emphasis 

DFREQ 

Integer 

Ending  de-emphasis 
frequency. 

FLEN 

Integer 

Number  of  points  in  FFT 
window. 

I  HAM 

Integer 

Switch  to  choose  window. 

HEADER 

Integer  Array 

Holds  header  values. 

NORM 

Integer 

Switch  to  select  energy 
normalization. 

WRTD 

Integer 

Switch  to  select  output 
option. 

TEST 

Integer 

Switch  to  select  input 
option. 

FI  LEI 

Integer  Array 

To  hold  input  file  name. 

FILE2 

Integer  Array 

To  hold  output  file  name. 

I START 

Integer 

First  input  disk  block. 

I  LAST 

Integer 

Last  input  disk  block. 

LTSTDO 

Integer 

Last  time  slice  to  do. 

VARIABLE 

TYPE 

PURPOSE 

COUNT 

Integer 

First  time  slice  to  do. 

WIND 

Real  Array 

Holds  window  values. 

PREM 

Real  Array 

Holds  emphasis  values. 

FREQ1 

Integer 

First  component  to  pre¬ 
emphasize. 

DFREQ 

Integer 

Last  component  to  de- 
emphasize. 

IDSP 

Integer  Array 

Holds  time  file. 

DTARAY 

Real  Array 

First  Array  Processor 
Memory  ( APM)  vector. 

COUNTDOWN 

Integer 

Counts  time  slices  done. 

WRTDA 

Integer 

Output  option  switch. 

NBLKS 

Integer 

Number  of  disk  blocks  to 
do. 

LASTCD 

Integer 

Counts  time  slices  to  do 

SWITCH  SETTINGS: 

SWITCH 

SETTING 

IDPRE 

= 

1  to  de-emphasize 

= 

0  to  not  de-emphasize 

I  HAM 

= 

1  for  Hamming  window 

= 

0  for  rectangular 

window 

HEADER( 58) 

S 

1  to  overlap  time 

slices 

= 

0  to  not  overlap 

NORM 

1  to  normalize  energy  in  time  slice 

= 

0  to  divide  spectrum  by  time  slice  energy 

WRTD 

s 

0  to  write  output 

to  disk  only 

10  to  write  spectrum  to  terminal  as  well 
as  disk 


12  to  write  spectrum  to  line  printer  as 
well  as  disk 


SWITCH 


SETTING 


TEST 

IOPTION 

IPRE 

SKIP 


1  to  create  an  integer  disk  file  of  tones 
to  simulate  speech 
0  to  process  speech 

0  for  ability,  to  set  parameters 
20-27  are  option  choices 

1  to  pre-emphasize 
0  to  not  pre-emphasize 

0  to  not  execute  APS  and  APM  initialization 
1  never 


RELATED  PROGRAMS: 

PLTA,  PLTN,  PLTS ,  PLTO,  DSTA,  DSTN,  CHOOS,  CHKJ,  CHKO, 
S128,  S64,  GFDB,  GRPH2 ,  FCTR,  PLTT,  DRVR. 
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UUUUUUUUUUOUOUOUUUOUUUUUUUUU 


CIHHI  SUBROUTINE  DRSO  NOTE:  THIS  SUBROUTINE  IS  POR  FORTRAN  5  !! 

THIS  SUBROUTINE  DRIVES  S12B,  S64.  IT  IS  CALLED  BY  DRVR.SV. 

BY:  CAPT  DAN  NARTIN 
DATS:  9/21/02 
SUBJ:  ACOUSTIC  ANALYSIS 

THIS  ROUTINE  SETS  UP  NECESSARY  FILES  AND  SWITCHES  FOR 
SUBROUTINES  SI20  AND  S44,  WHICH  COHPUTE  FOURIER  TRANSFORnS 
OF  SPEECH  FILES.  THE  OPERATOR  IS  PROMPTED  FOR  NECESSARY 
INFORMATION. 

DRSO  DOES  THE  FOLLOWING: 

1)  GETS  PARAMETERS  AND  OPTIONS 

2)  GETS  SPEECH  FILE 

3)  GETS  OBSERVATION  (SPECTRAL)  FILE 

4)  GETS  SPEECH  FILE  SEGMENT,  SPECIFIED  AS 
1ST  AND  LAST  DISK  BLOCKS 

5)  COMPUTES  WINDOU  VECTOR 
4)  COMPUTES  EMPHASIS  VECTOR 
7)  CONSTRUCTS  SPECTRAL  FILE  HEADER 
0)  CONSTRUCTS  TEST  SPEECH  FILE  IF  THAT  OPTION 
NAS  SELECTED 

9)  SETS  UP  WINDOW  AHD  MEMORY  HAP  FOR 
ARRAY  PROCESSOR 

10)  READS  SPEECH  FILE  INTO  EXTENDED  MEMORY 

11)  CALLS  EITHER  S12A  OR  S64  TO  GET  SPECTRUM 

12)  WRITES  SPECTRUM  FROM  EXTENDED  MEMORY  TO 
THE  OBSERVATION  FILE 

13)  READS  MORE  SPEECH  IF  THERE  IS  MORE  TO  DO 

14)  RETURNS  TO  DRVR 

FOR  MORE  INFORMATION,  SEE  THE  USERS  MANUAL  OR  MY  THESIS. 

SUBROUTINE  DRSO 
OVERLAY  ODRSO 

INCLUDE  "ARRAYP:F5APS.FR“ 

PARAMETER  LN*512 
PARAMETER  NDP'1024 
PARAMETER  HLH*LH/2 
REAL  DTARAY  ,B,PREM,WIND , SFREQ , I  DOB 

INTEGER  FILE1,I0PTN,IPRE,IDPRE, IDB, IHAH,FLEN, SKIP, HEADER ( 256 ), IDSP 
INTEGER  ISTART,1LAST, COUNT, FILE2, SMOOTH, WRTD, TEST, DUMMY, FREQ 
INTEGER  NBLKS,IFRQ(10),IAMP(10),XNEN,WRTDA,FLAG,CB1,CNTBR 
INTEGER  INITHEH ,SXHEH ,GSPCT ,LTSTDO , COUNTDOWN , NORM ,DFREQ ,DDB 
COMMON  / APH/  DTARAY(NDP) ,B(NDP) 

COMMON  /  VALS  /  IDSP(NDP) ,1D0B(NDP) , IHAM , IPRE , WIND ( LM ) 

COMMON  /  VALS  /  PREH(LM) , INITHEH, SXMEM,DUMMY(304) , COUNTDOWN 
COMMON  /  VALS  /  LTSTDO 

COMMON  /  VALT  /  SKIP, COUNT, SMOOTH, WRTD, NBLKS,XMEH, SFREQ, FLEN 
COMMON  /  VALT  /  FILEK13) , F I LE2 (13)  ,FLAC,CNTBR,1SXREH 
COMMON  /  VALT  /  JAOUT, LASTED 
COMMON  /  VALV  /  CB1(0:CBNAX) 

COMMON  /  VALU  /  NORN 


C 

4 


FLAC*0 


Define  variables  used 


F 


XHEN'ISXHEH 

£••••11  GET  SPECTRAL  OPTIONS 

S  TYPE*  Several  sets  of  defaults  exist  for  the* 

TYPE"  following  paraseters:  “ 

TYPE"  1)  PREEHPHASIS  <YES/N0>“  t 

TYPE"  2)  PREEHPHASIS  SLOPE  (dB/octave)" 

TYPE"  3)  PREEHPHASIS  STARTING  FREQUENCY- 

TYPE"  4)  DEEHPHASIS  (YES/HQ)" 

TYPE"  5)  DEEHPHASIS  SLOPE- 

TYPE"  6)  LAST  FREQUENCY- 

TYPE"  7)  FFT  VECTOR  LENGTH- 

TYPE"  A)  WINDOW  TYPE- 

TYPE"  They  ares  “ 


TYPE- 

OPTION 

NUMBER" 

TYPE" 

•21 

•22 

•23 

•24 

•  25 

•  26- 

TYPE"PREEHPH 

Y 

N 

Y 

N 

Y 

NO" 

TYPE-SLOPE 

10 

NA 

10 

NA 

10 

NA" 

TYPE"START 

500 

NA 

500 

NA 

500 

NA" 

TYPE"DEEHPH 

Y 

N 

Y 

N 

N 

N" 

TYPE-SLOPE 

10 

NA 

10 

NA 

NA 

NA¬ 

TYPE-LAST 

300 

NA 

300 

NA 

NA 

NA- 

TYPE-FFT  LEN 

64 

64 

126 

126 

128 

64" 

TYPE-UINDOU 

HH 

HH 

HR 

HH 

HH 

RECT“ 

TYPE"  If  you  choose  to  enter  values  for  these  paraseters  fros" 
TYPE"  the  Keyboard^  enter  nuseral  zero*  If  you  choose  a  default" 
ACCEPT"  option,  enter  the  option  desired:  ",I0PTN 

IF ( IOPTN.GT. 20. AND.IOPTN>LT >27) GOTO  2050 


ACCEPT-PREEHPHASIZE  HIGH  FREO  ?  ( 1  =  YES/0  =  M0> :  “,IPRE 

IF! I PRE.EQ.l ) ACCEPT "ENTER  dB/OCTAVE:  ",IDB 

IF ( IPRE.EQ. 1) ACCEPT" ENTER  STARTING  FREQ:  " ,FREQ 

ACCEPT" DEEHPHASIZE  LOU  FREQ?  4 1 *YES/0-N0) :  ",IDPRE 

IF( IDPRE.CQ.l ) ACCEPT “ENTER  DB/OCTAVE:  ",DDB 

IF ( IDPRE.EQ.  1 ) ACCEPT "ENTER  LAST  FREQ:  "  ,DFREQ 

TYPE"ENTER  FFT  VECTOR  LENGTH  (SAHPLES/VECTOR) .  LENGTH  RUST" 

ACCEPT"BE  A  HULT1PLE  OF  2  AND  IN  THE  RANCE  CIA, 10241:  "  ,FLEN 

TYPE-UHAT  TYPE  OF  UINDOU  IS  DESIRED?" 

ACCEPT-HAHHING  /  RECTANGULAR  (1/0):  “,IHAH 
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GOTO  SO 

20S0  IF( I0PTN.EQ.21IC0T0  21 

IF(I0PTN.E0.22)G0T0  22 
1 F  4  20PTN » EO « 23 ) COTO  23 
IF(I0PTN.EQ.24)G0T0  24 
IF(I0PTN.EQ.2S)C0T0  25 
IP< IOPTN.EQ.24)COTO  26 
COTO  SO 

21  IPRE‘1  ;preesphasize 


> 


ioo«io 

FREQ*500 
IDPRE* 1 
DS  B  *10 
DFREQ*300 
FLEN*84 
I HAH* 1 
GOTO  50 

22  XPRE*0 
IDB*0 
FREQ *500 
IDPRE z  0 
001*0 
OFREQ  =  0 
FLEN*84 
IHAH*1 
GOTO  50 

23  IPRE*1 
IDB*10 
FREQ=500 
IDPRE*1 
000*10 
DFREQ*300 
FLEN* 128 
IHAIt* 1 
GOTO  50 

24  IPRE*0 
100*0 
FREQ*SOO 
IDPRE*0 
000*0 
DFREQ*0 
FLEN* 128 
IHAH*1 
COTO  50 

25  IPRE*1 
100*10 
FREQ>500 
I0PRE*0 
000*0 
DFREQ*0 
FLEN*128 
1HAN* 1 
COTO  SO 

28  I PR  E  *  0 

100*0 
FREQ*500 
I0PRE*0 
000*10 
DFREQ*300 
FLEH*64 
IHAN*0 
GOTO  50 


{dO/octavt 

;start  prctaphasis  hart 


I 
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Ct* t ttt  SET  PARAMETERS 
SO  CONTIMUE 

SO  415  JS > 1 / 256 
41S  HEADER! J5)*0 

IF(FLEN.NE. 128)C0T0  52 
TYPE“»*ENTER  T  FOR  OVERLAPPING" 

ACCEPT***  'O'  FOR  NON-OVERLAPPING:  "  ,HEADER<58) 

52  CONTIMUE 

TYPE"*»*ENTER  '1'  TO  NORMALIZE  ENERGY  IN  FILE  TO  UNITY* 

ACCEPT*  OR  'O'  TO  DIVIDE  SPECTRUM  BY  FILE  ENERCY:  “,NORM 
PARAMETER  L’FLEN 
PARAMETER  HL*L/2 
C«****«  GET  OUTPUT  OPTION 

TYPE“*ENTER  OUTPUT  OPTION:  * 

TYPE"*»0  TO  URITE  SPECTRUM  TO  DISK  (AND  NO  OTHER  OUTPUT),* 
TYPE*»»lO  TO  WRITE  SPECTRUM  TO  SCREEN  AS  WELL  AS  DISK,* 
ACCEPT***12  TO  PRINT  SPECTRUM  AS  NELL  AS  WRITE  TO  DISK.  *,WRTD 
C»*»**»  GET  TEST  OPTION 

TYPE*»CET  TIME-FILE  TYPE* 

TYPE"»*ENTER  A  '1'  TO  CREATE  A  DISK  FILE  TO  SIMULATE  SPEECH:  * 
ACCEPT“*»ENTER  A  'O'  TO  CRUNCH  SPEECH:  “JEST 
IF (TEST . EQ . 0 )COTO  185 
ACCEPT**ENTER  TEST  DISK  FILENAME:  * 

READ ( 1 1 , 10 ) F I  LEI < 1 > 

8  CALL  0PEN(5,FLLE1,2,IER) 

IF(IER. HE. 1)TYPE* ERROR  ON  OPEN  OF  TEST  DISK  FILE,  IER=  ",IER 
CALL  CHKO(IER)  ;SEE  IF  UNIT  •  IN  USE 

1FUER.NE.  13)C0T0  9 

TYPE"FILE  DOES  NOT  EXIST.  WILL  CREATE  IT  FOR  YOU.* 

CALL  CFILW( FILE1 ,2,IER1 ) 

IF (IER1.NE.1) TYPE* ERROR  ON  FILE  CREATION.  IER1*  *,IER1 
CALL  CHECK* 1ER1) 

GOTO  8 

9  CONTINUE 
GOTO  190 

C**«*»»  GET  SPEECH  FILENAME  (IF  NOT  TEST) 

185  ACCEPT  "ENTER  NAME  OF  SPEECHFILE:  * 

READ (11 , 10  >  FILE1 4 1 ) 

CALL  OPEN ( 5 , F I LE 1 , 1 , XER1 ) 

IF ( 1ER 1 .NE. 1) TYPE* ERROR  ON  OPEN  OF  SPEECHFILE,  IER1*  *,IER1 
CALL  CHKO(IERl)  ?SEE  IF  UNIT  •  OPEN 

190  CONTINUE 

C«****i  GET  OBSERVATION  FILENAME 

ACCEPT'ENTER  NAME  OF  OUTPUT  DISK  FILE:  * 

REAKll  ,10)FILE2(1) 

10  F0RHAT(S13) 

12  CALL  0PEN(4,FILE2,2,IER) 

IF(1ER.NE.1)TYPE“ERR0R  ON  OPEN  OF  OBSERVATION  FILE,  IER*  *,IER 
CALL  CHKO(IER)  ;SEE  IF  UNIT  •  OPEN 

IF(IER«NE«13)C0T0  15 

TYPE’FILE  DOES  NOT  EXIST.  WILL  CREATE  IT  FOR  YOU.*' 

CALL  CFI LW( FILE2 ,2 , 1 ER1 ) 

IF (1ER1.NE.1) TYPE* ERROR  ON  FILE  CREATION.  IER1 «  *,IER1 
CALL  CHECK(IERl) 
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GOTO  12 
IS  CONTINUE 

C»««»»»  SPECIFY  PORTION  OF  SPEECH  FILE  TO  GET  SPECTRUH  OF,  OR  LENGTH 
C  OF  TEST  INPUT  FILE 

TYPE"»*SPECIFY  SPEECH  FILE  SECHENT  “ 

TYPE" (MULTIPLE  OF  4  BISK  BLOCKS  IF  TEST  FILE) “ 

IF (TEST .EQ.O)COTO  II 

TYPE“TEST  TIME-DOMAIN  FILE  HILL  START  AT  BLK  «0" 

ISTART«0 
GOTO  19 

11  ACCEPT**»»ENTER  FIRST  BLOCK  TO  BE  REAR:  “ # I  START 

19  ACCEPT“»**ENTER  LAST  BLOCK:  “ , HAST 

IF(HEADER(S4) .EQ»0)LTSTDQ>254*( 1+ILAST)/L 

{LAST  TINE-SLICE  TO  DO 

IF(HEADER(54) .EQ<0)CQUNT1ISTART*254/L+1 

;  1  ST  TINE-SLICE  TO  DO 

IF (HEADER (54) • NE • 0 ) LTSTDO  =  256» ( l  +  ILAST) /HL 

; LAST  TINE-SLICE  TO  DO 

IF(HEADER(54).NE.0)C0UNT*ISTART»254/HL+1 

;  1ST  TINE-SLICE  TO  DO 

C*Mt»«  COMPUTE  HINDOU  VECTOR 

TYPE“CONPUTINC  HINDOU  VECTOR  (IF  REQUESTED)” 

IF ( IHAN ,NE . 1 ) GOTO  210 
DO  200  J«1,L 

HINDI J) *  0.54  -  (0 .46*C0S( 4.2431 453072X J-l >/L) ) 

200  CONTINUE 

C««M»»  CALCULATE  EMPHASIS  VECTOR 

210  TYPE”CALCULATING  EMPHASIS  VECTOR  (IF  REQUESTED)” 

IFdPRE.NE.l .AND. IDPRE.NE.DCOTO  404 
DO  405  IN-1 ,HL 
405  PREH(IN)d.O 

IF(IPRE.NE.1)G0T0  295  {TO  NOT  PREENPHASIZE 
FREQ1MFREQ/SFREQ)*L+1  {FIRST  FREQ  PREENPH 
DO  400  IN-FREQ1,HL 

400  PREM(IN) * ( IN/FREQ1 )**( 0.1440944*2 *  I DB) 

295  CONTINUE 

IFdDPRE.NE.DCOTO  350  {TO  NOT  DEEHPHASIZE 
FREQ1* (DFREQ/SFREQ)»L  {LAST  FREQ 
DO  300  IN«1,FREQ1 

300  PREM( IN) « ( 1N/FREQ1 ) »• ( 0 . 1440944*2*DDB) 

350  CONTINUE 

404  CONTINUE 

CiiiiH  CONSTRUCT  HEADER  FOR  OUTPUT  FILE 

DO  410  J5«l,13 
410  HEADER  («I5)*F1LE2(J5) 

DO  412  JS|14,24 
J4*J5-13 

HEADER (J5)*FILE1(J4) 

HEADER (27>*IPRE 
HEADER (24) > I DB 
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HEADER  C  29  >  sPREO 
HEADER(30)*FLEN 
HEADER(31 >  *IHAH 
HEADER ( 33) *TEST 
HEADER(55)*C0UHT 
HEADER! 56) «LTSTDO 
HEADER  <  57 ) *HL 
HEADER! 60)* I DPR E 
HEADERUl)  *DDB 
HEADER(62)*DFREQ 
IF< IDPRE.EQ* 1 ) IPRE>1 
A  URTDA*WRTD 

1F(URTD>EQ.0)URTD*10 
UR ITE (URTD ,950) HEADER ( 1 ) 

URITE(WRTD,952)HEADER(14) 

WRITE ( HRTD  ,954 ) 

UR ITE (URTD,  956) ( HEADER (J5) , J5*27 ,31 ) 

URITE(URTD,955) 

WRITE ( URTD ,956) (HEADER! J5> ,  J5*60,62) 

UR1TE (URTD ,956) HEADER (33) 

416  IFdEST.EQ.l  )C0T0  430 

420  C0T0  500 

Ciiim  FILL  TEST  DISK  FILE 
430  IS*1 

DO  432  JJ1*1,10 
IFRQ(JJl) *0 
432  IAHP( JJ1 )  *0 

ACCEPT"«ENTER  1ST  FREQ  IH  HZ:  “ ,IFRQ( IS) 

445  ACCEPT" *EHTER  AMPLITUDE:  “,1AHP(IS) 

ACCEPT"»»»«»UANT  ANOTHER  TONE?  ( l-YES/O-NO) :  ",IT 

I F ( IT .EQ.O)COTO  455 

IS*IS*1 

IFUS.GT.  10) GOTO  455 

ACCEPT"*ENTER  ANOTHER  FREQ  IN  HZ:  ",IFRQ(IS> 

GOTO  445 

455  IF( IS.CT.10)IS>10 
HEADER  4  34 ) *  I S 

DO  456  J5*35 ,44 
J4*J5-34 

456  HEADER (JS)sIFRQ(J4) 

DO  457  J5*45 ,54 
J4»J5*44 

457  HEADER! J5 ) *IAHP(J4) 

WRITE (URTD, 962) HEADER (34) 

WRITE (URTD ,960) ( HEADER (J5) , J5>35 ,HEADER(34)*34 ) 

UR ITE (URTD, 960)  (HEADER! J5)  ,J5* 45 , HEADER ( 34 ) +4 4) 

45?  DO  495  JI«0,ILAST,4 

DO  460  IT*1  ,NDP 
460  IDSP( IT) *0*0 

DO  470  IB-1  ,IS 
DO  450  IT«1,NDP 

450  IDSP( IT) -IAHP ( IQ ) «COS ( 6 . 2 A3 1653072* I  FRO ( I Q ) ■ ( I T - 1 ) /SFREQ ) ♦ I  DSP ( I T ) 

470  CONTINUE 

CALL  WRBLX(5,JI,IDSP,4,IER> 


*■ 


495  CONTINUE 


500  CALL  FGTIHE(IHOUROf IHIHO  * ISECO) 


SIIP'O 

C««M«I  INITIALIZE  APS  AND  APH  HAP  FOR  ERDB 
IF< SKIP .EQ . 0) GOTO  502 
SKIP'O 

CALL  APIHITLNIL/DTARAY ,9(INITNEH,IER) 

CALL  API1APjDTARAY,0,4,IER> 

502  JI*1  {DISK  BLOCK  TO  WRITE  TO 

JIH'i  {DISK  BLOCK  TO  WRITE  TO 

503  IF<HEADER<50).EQ.0)XHEH'XREH>4  {GET  a  AVAIL  GTR  BLKS  EXT  HER 
I F ( HEADER ( Sft ) . NE .0)XHEH-X(1EH*2 

504  COUNTDOWN1 HEADER (56) -HEADER <55)^1 
TYPE-COUNTDOWN'" (COUNTDOWN 
IF(WRTDA.EQ.Q)WRTD:0 

505  NBLKS'ILAST-ISTART+1 
IF(NBLKS.LE.XHEH)COTO  510 

NBLKS'XHEH  {a  QTR  BLKS  CAN  DO 

510  IF<NEADER<5fl).NE.0)G0T0  512 

C0UNTc256>ISTART/L  {INITIALIZE  TIRE  SLICE  COUNT 

CALL  ERDB(5,ISTART,34, NBLKS, ICNT,IER) 

GOTO  530 

512  C0UNT«256>ISTART/HL 

CALL  ERDB ( 5 , 1ST ART , 36+XHEH , NBLKS , I CNT ,  I ER ) 

530  CONTINUE 

IF(IER.NE.9)C0T0  050 

Ti'PE“READ  END  OF  FILE!  SUCCESSFULLY  TRANSFERRED“/ICNT OTR  BLOCKS." 

NBLKS'ICNT 

HAST  =  ISTART+HBLKS- 1 

TYPE” PROCEEDING  UITH  NBLKS,"  OTR  BLOCKS  TRANSFERRED. “ 

I  ER  *  1 

050  IF(IER.EQ.1)G0T0  060 

TYPEHERROR  ON  ERDB,IER=  ",IER 
GOTO  1000 

060  ISTART'ISTART+NBLKS 

2000  F0RHAT(5C12.4) 

IF (FLEN.EQ. 120 .AND. HEADER (50) .NE .0 )COTO  090 
LASTCD 'COUNTDOWN -NBLKS >256/ L 
CALL  S64 

CALL  EWRB(4, JIH, 36, NBLKS, ICNT,IER> 

GOTO  904 

090  LASTCD* COUNTDOWN * NBLKSa 25 6/HL 

CALL  S120 

902  CALL  EWRB<4,JI,36,NBLXS>2,ICNT,IER> 

904  IF  HER . NE. 1) TYPE" ERROR  ON  EWRB,  1ER=  “,IER 

IF ( COUNTDOWN *LE.O. AND. I ER .EQ. 9) NBLKS* ICNT 

905  J1«JI+NBLXS>2 
JIH* JIHfNDLKS 

TYPE"**»COUNTDOWN*“ (COUNTDOWN 
IF (COUNTDOWN. LE.O) GOTO  911 

910  IFUSTART.LE.ILASTICOTO  505  {MORE  SPEECHFILE  TO  DO 
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911 


CONTINUE 

IF<FLEN.E0.120.AND.HEADER(5B).NE.0)C0T0  912 
HEADER! 59) < JIH  ;«  OB'S  IN  FILE 

GOTO  915 
912  HEADER ( 59 ) s  J 1 

915  CALL  URBLK (4,0 /HEADER, 1, I ER)  ;URITE  HEADER  TO  1ST  QBLK 

CALL  RESET 
URTDSURTDA 

CALL  FCT1HE(IH0UR,IM1N,1SEC) 

ISECT  *  ISEC  -  ISEC0-f60*(  IMIH  -  IHIHO  )  +  3600»  ( IHOUR  -  IHOURO ) 

TYPE'TIME  TO  GET  SPECTRUM  IS", ISECT,"  SECONDS" 


950  FORMAT!"****  THIS  FILE'S  HARE  IS:  ",S13> 

952  FORMAT!"  TIME  DOMAIN  FILE  NAME:  “,S13) 

954  FORMAT ! 7X, " 1PRE" ,6X,“IDB" ,6X , “FREQ" ,6X,"FLEN“  ,6X,"IHAM“ ) 

955  FORMAT (7X ,“IDPRE"  ,7X ,“DDB“ ,4X,“DFREQ‘’ ) 

956  FORMAT ! 61 10 ) 

950  FORMAT!"  -  TEST  3  “,14“ 

940  FORMAT!"  ",1017) 

962  F0RHAT12X ,14,"  TONES  IH  TEST  FILE.  HERE  ARE  FREQS  AND  AMPLITUDES") 
1000  CALL  OVEXIT ( ODRSQ , 1ER) 

CALL  CHECK! IER) 

RETURN 

END 


i  ' 


< 
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S:  S128  and  S64 

3UAGE:  FORTRAN  5 

E:  September  21,  1982 

MR:  D.  Martin 

JECT:  Fourier  Transform 

LING  SEQUENCE:  S128  or  S64 

E  OF  LAST  REVISION:  September  21,  1982 

POSE: 

Subroutine  S128,  called  by  DRSQ,  computes  overlapping 
point  Fast  Fourier  Transforms  ( FFTs )  on  an  integer 
e  residing  in  extended  memory.  Subroutine  S64  does  not 
rlap  the  input  window  and  will  do  FFTs  of  window  sizes 
:h  are  multiples  of  two  in  the  range  16-1024,  inclusive. 

sription: 

Location:  DP4: BRATCHET 

Size:  S128.FR 

S128.RB 
S64.FR 
S64.RB 

SRAM  USE: 

These  programs  are  called  by  DRSQ  and  call  no  sub¬ 
bines.  They  use  the  Array  Processor  and  extended  memory 
naximize  the  speed  with  which  they  calculate  FFTs.  These 
routines  are  similar  and  the  flowchart,  Figure  73,  appli 
x>th.  The  two  differences  between  them  are:  (1)  the 
low  size  of  S128  is  fixed  at  128  points  while  that  of  S64 
(electable,  (2}  the  window  of  S128  is  moved  along  the 


4866  bytes 
5378  bytes 
3910  bytes 
4276  bytes 


input  file  in  increments  of  64  points  at  each  FFT  calcula¬ 
tion  while  the  window  in  S64  is  moved  along  the  input  file 
in  increments  of  its  window  size. 

These  differences  require  that  extended  memory  be 
managed  differently  in  the  two  routines.  For  S64,  an  input 
window  of  N  integer  points  returns  N/2  real  components  which 
replace  the  input  time  slice-by-time  slice.  The  spectrum 
is  stored  on  disk  in  contiguous  units  called  observations. 
The  first  element  is  the  observation  energy,  E,  computed  as 
the  sum  of  the  squares  of  spectral  components  one  through 
N/2-1.  The  Array  Processor  FFT  routine  returns  N/2+1  com¬ 
ponents;  but  prior  to  the  energy  computation,  the  zero  and 
N/2  components  are  set  to  zero.  Each  observation,  then, 
consists  of  an  energy  value  as  its  first  element,  and 
spectral  components  one  through  N/2-1  in  the  remaining  N/2-1 
elements.  The  Kth  element  of  the  observation  is  the  Kth 
spectral  component  of  the  N-  point  time  slice. 

The  observations  are  written  to  disk  starting  at  disk 
block  number  one.  Each  N-  point  time  slice  requires  N 
words  of  memory.  Taking  advantage  of  this,  S64  stores  the 
spectrum  back  into  the  storage  locations  previously 
occupied  by  the  window.  In  this  sense,  the  FFTs  are  com¬ 
puted  in-place. 

%  But  because  S128  overlaps  its  window  by  64  points  each 

FFT  calculation,  the  in-place  scheme  will  not  work.  The 
scheme  S128  uses  has  DRSQ  read  the  input  file  into  the 
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upper  half  of  available  extended  memory.  Then  as  the  input 
window  is  moved  along  the  upper  half,  the  output  is  stored 
in  the  lower  half  of  extended  memory,  starting  at  the  first 
available  location.  Hie  input  and  output  storage  locations 
overlap  only  at  the  last  locations  at  the  top  of  extended 
memory. 

Both  routines  transform  an  integer  input  file,  such  as 
that  produced  by  the  Cromemco  analog- to-digital  conversion 
system  in  the  Speech  Laboratory,  into  a  real  disk  file. 

Disk  block  number  zero  is  a  header  containing  information 
which  describes  and  identifies  the  file.  This  header  is 
the  same  for  both  S128  and  S64,  and  its  assignments  are 
listed  in  Table  3. 

LIST  OF  VARIABLES: 


VARIABLE 

TYPE 

PURPOSE 

JAOUT 

Integer 

First  element  of  output 
vector. 

DUMMY ( 1 ) 

Integer 

Last  memory  block  read 
from. 

J6 

Integer 

First  memory  block  read 
from. 

J6S 

Integer 

Signals  first  call. 

XMEM 

Integer 

Extended  memory  available 
for  input. 

DTARAY 

Real  Array 

Array  Processor  memory. 

I  DSP 

Integer  Array 

Input  array  mapped  through 
extended  memory. 

WIND 

Real  Array 

Contains  window  to  be 
applied  to  input. 

VARIABLE 

TYPE 

PURPOSE 

PREM 

Real  Array 

Contains  emphasis  vector 
applied  to  spectrum. 

NORM 

Integer 

Switch  to  select  energy 
normalization. 

IPRE 

Integer 

Switch  to  select  emphasis 

I  HAM 

Integer 

Switch  to  select  window 
type. 

I  DOB 

Real  Array 

Output  array  mapped 
through  extended  memory. 

SWITCH  SETTINGS; 


SWITCH 

SETTING 

I  HAM 

=: 

1 

for  Hamming  window 

0 

for  rectangular  window 

IPRE 

= 

1 

to  emphasize  spectrum 

0 

to  not 

NORM 

1 

to  normalize  spectral  components  to  observa¬ 
tion  energy 

0 

to  divide  each  spectral  component  by  the 
observation  energy 

RELATED  PROGRAMS: 

PLTA,  PLTS ,  PLTN,  PLTT,  PLTO,  DRVR,  DSTA,  DSTN,  DRSQ, 


CHKJ,  CHKO ,  CHOOS,  FCTR,  GFDB ,  GRPH2 


-»..un  i  . . . 


HEADER  ELEMENT 

CONTENT 

1-13 

Observation  file  name. 

14-26 

Speech  file  name. 

27 

IPRE,  pre-emphasis  switch. 

28 

IDB,  pre-emphasis  rate  (dB/Octave). 

29 

FREQ ,  frequency  at  which  pre¬ 
emphasis  starts. 

30 

FLEN,  FFT  window  size. 

31 

IHAM,  window  switch. 

32 

Unused  (set  to  zero). 

33 

TEST,  test  option  switch. 

34-54 

Unused  (set  to  zero). 

55 

COUNT,  number  of  first  time  slice 
to  do. 

56 

LTSTDO,  number  of  last  time  slice 
to  do. 

57 

HL,  number  of  elements,  per 
observation. 

58 

Overlapping  switch. 

59 

Number  of  disk  blocks  in  observa¬ 
tion  file. 

60 

IDPRE,  de-emphasis  switch. 

61 

DDP,  de-emphasis  rate  (dB/Octave). 

62 

DFREQ,  frequency  at  which  de-emphas: 
ends . 

63-256 

Unused  (set  to  zero). 

TABLE  3.  Header  Assignments  for  Spectral  File  Computed  by 
DRSQ  Calling  S128  or  S64. 
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8UIR0UTIHE  S12B  NOTE:  THIS  ROUTINE  IS  FOR  FORTRAN  5  !! 

THIS  SUBROUTINE  CALLS  ARRAY  PROCESSOR  ROUTINES  TO  CALCULATE  12B 
POINT  FFT'S.  IT  OVERLAPS  THE  TINE  WINDOWS  BY  44  POINTS. 

BY x  CAPT  DAN  NARTIN 
DATE:  9/21/62 
SUBFt.  ACOUSTIC  ANALYSIS 

THIS  ROUTINE  IS  CALLED  BY  DRSQ  WHICH  SETS  UP  THE  NECESSARY 
FILE  INPUT/OUTPUT  FILES  AND  OPTIONS.  THE  INPUT  TO  THIS  ROUTINE 
RESIDES  IN  EXTENDED  HEHORY  AND  THE  OUTPUT  IS  DEPOSITED  THERE. 
S12B  DOES  THE  FOLLOWING: 

1)  INITIALIZES  PARAHETERS 

2)  HAPS  OUTPUT  ARRAY  IDOB 

3)  HAPS  INPUT  ARRAY  IDSP 

4)  LOADS  TINE  SLICE  INTO  DTARAY 

5)  WINDOWS  TINE  SLICE 
4)  CALCULATES  FFT 
7)  ZEROS-OUT  1ST  AND  LAST  ELEMENTS 

6)  EMPHASIZES  SPECTRUM 

9)  NORMALIZES  SPECTRUM 

10)  LOADS  SPECTRUM  INTO  IDOB 

11)  REMAPS  IDOB  IF  NECESSARY 

12)  REMAPS  IDSP  IF  NECESSARY 
FOR  MORE  INFORMATION,  SEE  THE  USERS  MANUAL  OR  MY  THESIS. 
SUBROUTINE  S12A 

INCLUDE  ”ARRAYP:F5APS.FR~  ;  APS  PARAMETER  FILE 


C»«**»*  Set  up  variables  to  be  used. 

PARAMETER  LR’64 
PARAMETER  TLH*12ft 
PARAMETER  NDP«1024 
PARAMETER  HAXLH‘512 
REAL  DTARAY | B|PREN|WIHD,SFREQ 
REAL  SUME, FACTOR, I DOB 

INTEGER  CBt, DUMMY, FLEN,FILE1,FIIE2, FLAG, CNTBR, IDSP 
INTEGER  IHAN,IPRE, SKIP, COUNT, SMOOTH, WRTD,NBLKS,XNEH 
INTEGER  INITMEN,SXMEH,CSPCT,LTSTDO, COUNTDOWN, NORM 

COMMON  /APH/  DTARAY  CHOP) ,B(HDP)  ;AP  MEMORY 

COMMON  /  VALS  /  IDSP(NDP) ,IDOB(NDP),IHAH,IPRE,WIND(HAXLM) 
COMMON  /  VALS  /  PREN(HAXLH) , IN ITMEM ,SXMEM,DUMMY(304) , COUNTDOWN 
COMMON  /  VALS  /  LTSTDO 

COMMON  /  VALT  /  SKIP, COUNT, SMOOTH, WRTD,NBLKS,XRER,SFREQ,FLEN 
COMMON  /  VALT  /  FILE1 ( 13) ,F1LE2( 13) , FLAG, CNTBR ,ISXHEM 
COMMON  /  VALT  /  JAOUT ,LASTCD 
COMMON  /  VALV  /  CB1 (OtCBMAX) 

COMMON  /  VALU  /  NORM 


<:••••••  INITIALIZE  APS  AND  APN  MAP 

IF(SK!P.EQ.0)60T0  230 

CALL  APINITt NIL, DTARAY ,9,INITMEM,IER) 

CALL  APMAP( DTARAY ,0,4,1ER) 

<  230  CONTINUE 

JAOUT'l  ; 1  ST  ELEMENT  OF  OUTPUT  VECTOR 
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t 


240 

CM 

c»» 

c 

c 


Cm 

260 

300 


DUHHYt  1 ) <9*(XNEH+NBLKS)/4 

J4*XHEH/4*9 

J4S*J4 

J5«9 

CALL  APHAPdDOB ,  J5 , -2 ,  IER ) 
CONTINUE 

CALL  APNAPdDSP ,J6,-1 ,IER) 
ITERATE  TO  LINE  700 


{LAST  KEN  BLK  READ  FROH 
B 1 ST  HEN  BLK  READ  FROH 
{SIGNALS  1ST  ENTRY 

.{HAP  OUTPUT  ARRAY 

{NAP  INPUT  ARRAY 


THE  FIRST  ELEHENT  OF  DATA  DOESN'T  RESIDE  ON  A  IX  BOUNDARY; 
NECESSARILLY,  IN  THE  FIRST  NAP.  TO  ACCOUNT 
FOR  THAT  POSSIBILITY: 

IF<J6.NE,J4S)J1*1  {GET  1ST  ELEHENT  OF  IDSP 

IF( J4.EQ . J4S) J1*1+NDP<( FL0AT(XHEH)/4- lFlX(XHEH/4) ) 

J4*J4+1 

TRANSFER  126  PNTS  FROH  IDSP 
DO  300  J2*Jl;Jl4TLH-l 
J3*J2-(J1-1) 

DTARAY (J3)* IDSP (J2) 


CkftMM  UIHDOU  VECTOR  DTARAY 
305  1F(IHAN.NE.1)G0T0  310 
CALL  APSETL(TLH,IER) 

CALL  CBSET(CB1 ,CBAXR ,B ,CBAAHH , WIND, IER) 

CALL  VLDR(CBl)  {LOAD  WINDOW  VECTOR 

CALL  CBSETICBl ,CBAXR, DTARAY ,CBAYR,B,CBAZR, DTARAY, IER) 
CALL  VHRA(CBl)  {NULTIPLY  VECTORS 


310  CONTINUE 

C»M»*i  CALCULATE  FFT  OF  VECTOR  DTARAY 
CALL  APSETL(LH;IER) 

CALL  CBSET(CB1,CBAXC, DTARAY, CBCW,CWDFT, IER) 

CALL  VFFTC(CBl) 

CALL  VBRC(CBl) 

CALL  VFFTR(CBl) 

C>»MM  COHPLEX  DTARAY  CONTAINS  1ST  LH+l  SPECTRUH  COEFFICIENTS 
C  NEED  ZERO-OUT  1ST  ELEHENT,  I.E.,  ZERO-HERTZ  TERN 

C  ALSO  NEED  ZERO-OUT  45TH  ELEHENT,  I.E.,  4 KHZ  TERH 

DTARAYU)«0,0 
DTARAY  <  2 )  *0 . 0 

CALL  CBSETICBl, CBCW,CUSTD,CBAZR, DTARAY, IER) 

CALL  VSHAICBl ) 


400  lFdPRE.NE.DGOTO  404 

CALL  CBSET(CB1 ,CBAXR,B,CBAANH,PREH,IER) 

CALL  VLDR(CBl)  {LOAD  EHPHASIS  VECTOR 

CALL  CBSET(CB1,CBAZR, DTARAY, CBAYR, DTARAY, IER) 

CALL  VHRA(CBl)  {NULTIPLY  VECTORS 

Cmmm  get  energy  in  speechfile 
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CALL  CBSET!CB1,CBAAHH, SURE, CBAXR, DTARAY, IER) 

CALL  VSER(CBl)  {SUM  VECTOR  ELEMENTS 

SORT  SPECTRAL  COMPONENTS 
BO  410  J4«1,LM 
DTARAY! J4) *SQRT! DTARAY (J4) ) 

NORMALIZE  ENERCY  IN  SPECTRAL  FILE  TO  A  CONSTANT 
OR  DO  SOMETHING  ELSE 
IF!N0RM.E0.0)C0T0  410 
FACTOR* 10000/ SORT (SURE) 

GOTO  420 

FACTOR* 10000000/SUME 
CONTINUE 

CALL  CBSET!CB1,CBSCR, FACTOR, CBAZR,B,IER> 

CALL  VMRS(CBl)  {MULTIPLY  VECTOR  BY  SCALER 

LOAD  DTARAY  WITH  NEXT  120  PNTS 
J1*JWLH  {MOVE  OVER  64  PNTS 

IF< J1 • E42 • 1025 ) GOTO  429  {NEED  TO  REMAP  TO  NEXT  EXT  MEM  BLI 
IF ( J 1  •  EQ «961) GOTO  422  {ONLY  64  PNTS  LEFT 

J3*0 

DO  421  J2*J1,J1*127 
J3* J3+1 

DTARAY (J3) *IDSP!J2) 

GOTO  429 
J3*0 

DO  423  J2* J1 ,1024  {LOAD  LAST  64  PNTS  OF  CURRENT 

J3*J3*1  {MEMORY  BLOCK 

DTARAY! J3)*IDSP!J2) 

IF ( J6>GT< DUMMY ( 1) )COTO  425  {LAST  MEM  BLOCK 
CALL  APHAP!IDSP,J6,~1,IER) 

DO  424  J2*65,120  {LOAD  1ST  64  PNTS  OF  NEXT 

J3«J2*64  {MEMORY  BLOCK 

DTARAY! J2) *IDSP(J3) 

CALL  APHAP( IDSP ,J6-1,-1, 1 ER > 

GOTO  429 

DO  427  <13*63,120  {NEED  FUDGE  LAST  64  PNTS 

DTARAY  (J3>*DTARAY!44)»EXP(64-J3) 

LOAD  NORMALIZED  SPECTRUM  INTO  IDOB 
C0UNT*C0UNT+1 
COUNTDONN*COUNTDONN-1 
IDOB! JAOUT) *SQRT(SUME) 

<13*1 

DO  430  J2*JA0UT+1 ,JAQUT+63  {LOADING  64  PNTS 
J3*J3»1 

IDOB! J2)*B(J3) 

JAOUT* JAOUT+LH 
IF! JAOUT.NE. 1025) GOTO  435 
JAOUT*  1 
J5*J5t2 

CALL  APMAP! IDOB, J5,* 2,1 ER) 

CONTINUE 
FORMAT (5C12  *4) 

IF  <  COUNTDOWN • LE  >  LASTCD ) GOTO  1000 
CONTINUE 

IF!J1.LE.961)G0T0  303 


« 
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SUBROUTINE  "44  NOTE*  THIS  ROUTINE  IS  FOR  FORTRAN  5  II 
THIS  SUIROUTINB  CALLS  ARRAY  PROCESSOR  ROUTINES  TO  CALCULATE 
FFT'S.  THB  UNION  SIZE  CAN  IB  ANY  NULTIPLE  OF  2  IN  THE  RAHCE 
114,10243.  ALTHOUGH  RESIGNS!  FOR  SPEECH  ANALYSIS,  THIS 
ROUTINE  IS  FLEXIBLE. 

BY l  CAPT  BAN  NARTIN 

BATE*  9/21/12 

SUBJ.  ACOUSTIC  ANALYSIS 

THIS  ROUTINE  IS  CALLEI  BY  IRS8  WHICH  SETS  UP  THE  NECESSARY  FILE 
INPUT/OUTPUT  OPTIONS.  THE  INPUT  TO  THIS  ROUTINE  RESIDES  IN 
EXTENIEO  NEHORY  ANB  THE  OUTPUT  IS  DEPOSITEI  THERE. 

S44  DOES  THE  FOLLOWING t 

1)  INITIALIZES  PARARETERS 

2)  HAPS  OUTPUT  ARRAY  IDOS 

3)  HAPS  INPUT  ARRAY  IDSP 

4)  LOADS  TIHE  SLICE  INTO  DTARAY 

5)  WINDOWS  TIHE  SLICE 

4)  CALCULATES  FFT 

7)  ZEROS-OUT  1ST  AND  LAST  ELEHENTS 

5)  EHPHAS1ZES  SPECTRUH 

9)  NORHALIZES  SPECTRUH 

10)  LOADS  SPECTRUH  INTO  IDO* 

11)  REHAPS  I DOB  IF  NECESSARY 

12)  REHAPS  IDSP  IF  NECESSARY 

FOR  HORE  INFORHATION ,  SEE  HY  THESIS  OR  THE  USERS  HANUAL. 


SUBROUTINE  $44 

INCLUDE  ’ARRAYFiFSAPS.fr*  ;  APS  PARAHETER  FILE 

Cttttta  Set  up  variables  to  be  used. 

PARAHETER  LNSFLEN 
PARAHETER  HLH-LH/2 
PARAHEU#  XLH’HlH-1 
PARAMETER  NDP«1024 
PARAHETER  HAXLH'512 
REAL  DTARAY ,B,PREH,WIND,SFREQ 
REAL  SURE, FACTOR, IDQB 

INTEGER  CB1 ,FLEN,FILE1 ,FILE2,FLAC,CNTBR,IDSP,INITNEH,SXHEH 

INTEGER  1HAH,IPRE,SIIP, COUNT, SHOOTH,WRTD,NBLIS,XHEH,DUHHY 

INTEGER  CSPCT,LTSTDO, COUNTDOWN, NORN 

COHHON  /APR/  DTARAY (NDP), B(NDP)  ;AP  HEHORY 

COHHON  /  VALS  /  I  DSP (NDP) , I  DOB (NDP) ,IHAH,IPRE,UIND(HAXLH) 

COHHON  /  VALS  /  PREH(HAXLH) ,INITHEH,SXHEH,DUHHY (304) , COUNTDOWN 

COHHON  /  VALS  /  LTSTDO 

COHHON  /  VALT  /  SKIP, COUNT ,SHOOTH,WRTD, NILES, XnEH,SFREQ,FLEH 
COHHON  /  VALT  /  FILEK13)  ,F1LE2 ( 13)  ,FLAC,CNTBR,1SXHEH 
COHHON  /  VALT  /  JAOUT,LASTCD 
COHHON  /  VALV  /  CB1<0:CBHAX> 

COHHON  /  VALU  /  NORH 

Cttitsi  INITIALIZE  APS  AND  APR  HAP 
IF(SKIP.EO.O)GOTO  250 
CALL  APINITtHIL, DTARAY, 9, INITHEHflER) 
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CALL  APHAP( BTARAY ,0,4, IER) 

J5*9  {HEM  BLOCK 

CONTINUE 

CALL  APHAPdBSP,J3,-l / 1 ER > 

CALL  APHAPdBOB,  J5,- 1 ,  IER ) 


{INPUT  ARRAY 
{OUTPUT  ARRAY 


JA0UT*1  {POINTER  IN  OUTPUT  VECTOR 

ITERATE  BO INC  FPT'S  BY  TINE-SLICE 
BO  700  J1*1,NDP,LH 

INPUT  SEQUENCE  IS  INTEGER,  OUTPUT  IS  REAL.  SO  OUTPUT 
TAKES  PLACE  OF  INPUT  ELEHENT  FOR  ELEHENT 
BO  300  J2«Jl,JWLH-l 
J3*J2-(J1-1) 

BTARAY (J3)*1BSP(J2) 


UIHBOU  VECTOR  BTARAY 
IFdHAH.NE.DGOTO  310 


CALL  APSETL(LH,IER) 

CALL  CBSET (CB1 ,CBAXR,B ,CBAARN,U1ND , IER ) 

CALL  VLDR(CBl)  {LOAD  UINDOU  VECTOR 

CALL  CBSET (CB1  ,CBAXR , BTARAY ,CBAYR ,B ,CBAZR ,DTARAY , I ER ) 
CALL  VRRA(CBl)  {MULTIPLY  VECTORS 

CONTINUE 

CALCULATE  FFT  OF  VECTOR  BTARAY 
CALL  APSETL(HLN,IER) 

CALL  CBSET ( CB 1 ,CBAXC,DTAR AY ,CBCU ,CUDFT , IER) 

CALL  VFFTC(CBl) 

CALL  VBRC(CBl) 

CALL  VFFTR(CBl) 

BTARAY  CONTAINS  1ST  HLH+1  SPECTRUM  COEFFICIENTS 


NEED  ZERO-OUT  1ST  ELEHENT,  I.E.,  ZERO-HERTZ  TERN 
ALSO  NEED  ZERO-OUT  33RD  ELEMENT,  I.E.,  4KHZ  TERN 
BTARAYt 1 > s  0 • 0 
BTARAY ( 2 ) *  0 • 0 

CALL  CBSET(CB1,CBCU,CNSTD,CBAZR, BTARAY, IER) 

CALL  VSHA(CBl)  {GET  SQUARE  MAGNITUDE 

IF(IPRE.NE.1)COTO  406 

CALL  CBSET ( CB1 , CBAXft ,B,CBAAMN,PREM , IER ) 

CALL  VLBR(CBl)  {LOAD  ARRAY 

CALL  CBSET(CB1,CBAZR, BTARAY, CBAYR ,DTARAY , IER) 

CALL  VNRA(CBl)  {MULTIPLY  VECTORS 


CALL  CBSET(CB1,CBAANN,SUNE,CBAXR, BTARAY, IER) 
CALL  VSER(CBl)  {SUN  ELEMENTS 

SORT  SPECTRAL  COMPONENTS 
BO  420  J4*2,HLM 
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420  DTARAY(J4)*SQRT(DTARAY(J4) ) 

£••»•••  NORMALIZE  ENERGY  IN  SPECTRAL  FILE  TO  A  CONSTANT 
C  OR  00  SOMETHING  ELSE 

IF(N0RH*E0«0)G0T0  42? 

FACTOR <10000/ SORT (SOME) 

GOTO  423 

422  FACTOR* 10000000/SUNE 

425  CONTINUE 

CALL  CBSETICBl fCBSCR * FACTOR f CD AZRtDTARAY , IER) 

CALL  VMRS(CBl)  {MULTIPLY  BY  SCALER 

Chhh  LOAD  NORMALIZED  SPECTRUM  INTO  IDOB 
COUNTDOWN*COUNTDOUN* 1 
C0UNT*C0UNT+1 

IDOB( JAOUT) *SQRT(SUHE)  {LOAD  ENERGY 
J3*l 

DO  430  J2*JAOUT+1|JAOUT+KLH  {LOADING  ONLY  HLH  POINTS 
J3*J3+1 

430  IDOBC J2) *DTARAY( J3> 

440  F0RMAT(SG12i4) 

IF< COUNTDOWN. LE.LASTCD) GOTO  1000 
700  JAOUT*JAOUT+HLH 

800  J5*J5+1 

C0T0240 

1000  RETURN 

END 


{ 
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FILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OF  LAST  REVISION: 

PURPOSE: 

These  programs  compute  distances  between  observations 
and  phonets.  DSTA  fills  the  output  disk  file  with  all  the 
distances  between  the  specified  observations  and  phonets. 
DSTN  fills  the  output  disk  file  with  one  unit  for  each 
specified  observation.  Each  unit  contains  the  observation 
number,  observation  energy,  worst  phonet  match,  maximum 
distance,  and  a  selected  number  of  ordered  pairs:  (phonet 
number,  next  best  distance). 

DESCRIPTION: 

Location:  DP4: BRATCHET 

Size:  DSTN. FR 

DSTN.RB 
DSTA.  FR 
DSTA. RB 

PROGRAM  USE: 

These  programs  are  called  by  DRVR  and  call  CHOOS,  GFDB, 
CHKO,  and  CHKJ.  They  use  the  Array  Processor  to  compute 
distances  between  observations  and  phonets  from  files 
calculated  by  DRSQ.  The  distinction  between  observation 


12912  bytes 
17402  bytes 
12699  bytes 
17488  bytes 


DSTA  and  DSTN 
FORTRAN  5 

September  21,  1982 
D.  Martin 
Acoustic  Analysis 
DSTA  or  DSTN 
September  21,  1982 


Le  and  phonet  file  is  arbitrary  in  format;  any  file  can 
run  against  itself.  The  distinction  in  the  way  these 
Jtines  handle  files  is  that  distances  to  all  phonets 
a  computed  for  each  observation. 

The  main  difference  between  these  routines  is  that 
Lie  DSTA  outputs  every  computed  distance,  DSTN  selects 
a  best  matches  and  outputs  the  observation  number,  energy, 
3  worst  match  as  well.  Tables  4  and  5  contain  the  output 
Le  formats. 

TVo  distance  rules  are  selectable:  Ml  and  M2.  Ml  is 
ikowski  One  distance: 


N 


ere: 


distance  between  the  ith  observation  and  the 
jth  phonets 


K  =  constant  scale  factor. 

X^  =  1th  component  of  the  ith  observation. 

t  h 

=  i  component  of  the  j  phonet. 
is  the  Minkowski  Two  distance: 


ere: 


distance  between  the  it^1  observation  and  the 
,th  , 

1  phonet. 


K  =  constant  scale  factor. 

=  lfc^  component  of  the  i^ 


observation 


1 


ELEMENT  NUMBER 

ASSIGNMENT 

1 

D(i.  j) 

2 

D( i ,  j+1) 

3 

• 

D(i,  j+2) 

• 

• 

K-  j+1 

• 

• 

D( i ,  K) 

K-  j  +  2 

D( i+1,  j) 

K-  j+3 

# 

D(i+1,  j+1) 

0 

0 

2K- j+1 

• 

• 

D( i+1,  K) 

0 

(1+1)K- j+1 

• 

D( i+1 ,  K) 

where:  i  =  starting  observation  number  specified. 


j  =  starting  phonet  number  specified. 
i  +  1  =  last  observation  number  specified. 
K  =  last  phonet  specified. 

TABLE  4.  DSTA  Output  File 
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ELEMENT  NUMBER 


ASSIGNMENT 


1 

Number  of  first  observation 
specified. 

2 

Energy  of  first  observation 
specified. 

3 

Number  of  worst  phonet  choice. 

4 

Largest  distance. 

5 

Number  of  best  phonet  choice. 

6 

Minimum  distance. 

7 

Number  of  next  best  phonet  choice. 

8 

Next  minimum  distance. 

• 

• 

m 

4+2N 

0 

0 

Nfc^  minimum  distance. 

4+2N+1  1 

1 

Number  of  next  observation  specified 

0 

• 

• 

K ( 4+2N) 

0 

• 

Nfc^  minimum  distance  to 
observation. 

TABLE  5.  DSTN  Output  File 


Yj^  =  component  of  the  j**1  phonet 
Hie  observation  energy  is  not  included  in  either  distance 
computation.  Also,  the  distance  is  checked  for  overflow 
before  being  truncated  to  an  integer  and  written  to  disk. 

If  the  distance  is  greater  than  32767,  it  is  hard  limited 
to  that  value.  Hie  scale  factor,  K  in  the  computations, 
is  adjustable  to  obtain  a  satisfactory  range  of  distance 
measures.  Hie  flowchart  in  Figure  74  applies  to  both  DSTA 
and  DSTN.  Distance  file  header  assignments  are  in  Table  6. 

LIST  OF  VARIABLES: 


VARIABLE 

TYPE 

PURPOS  E 

XMEM 

Integer 

Number  of  extended  memory 
blocks  to  use. 

OUTD 

Integer 

Output  option  switch. 

FILE2 

Integer  Array 

Holds  observation  file 
name. 

FI  LEI 

Integer  Array 

Holds  phonet  file  name. 

FILE3 

Integer  Array 

Holds  output  distance  file 
name. 

STS1 

Integer 

First  observation  to  do. 

STSL 

Integer 

Last  observation  to  do. 

PTS1 

Integer 

First  phonet  to  do. 

PTSL 

Integer 

Last  phonet  to  do. 

HEADER 

Integer  Array 

Holds  header. 

NSTSTB 

Integer 

Number  of  observations 
to  do. 

NPSTS 

Integer 

Number  of  words  in 
observation. 

Ente: 


A 


% 

FIGURE  74.  Flowchart  of  DSTA  and  DSTN 
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FIGURE  74  CONTINUTED.  Flowchart  of  DSTA  and  DSTN 

> 
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ELEMENT  NUMBER 

CONTENT 

1-13 

Distance  file  name. 

14-26 

Observation  file  name. 

27-39 

Phonet  file  name. 

40 

Number  of  first  observation  time 
slice  to  do. 

41 

Number  of  last  observation  time 
slice  to  do. 

42 

Number  of  first  phonet  time  slice 
to  do. 

43 

Number  of  last  phonet  time  slice 
to  do. 

44 

Number  of  disk  block  that  holds 
first  observation  to  do. 

45 

Number  of  disk  block  that  holds 
first  phonet  to  do. 

46 

Switch:  4  =  observation  and  phonet 

files  identical/0  =  different. 

47 

NCHOICES . 

48 

Number  of  observation  time  slices 
to  do. 

49 

Number  of  phonet  time  slices  to  do. 

50 

Number  of  elements  per  time  slice. 

51 

Number  of  extended  memory  block 
used  less  those  in  window. 

52-256 

Unused. 

TABLE  6.  Distance  File  Header 
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VARIABLE 
I  BO 

FSPFSQ 

NPTSTB 

NLDTB 

JCTLD 

IBP 

FDSTR 

NPFS 

NPFP 

NBFPH 

NBFOBS 

NBFDIST 

OBSMAPCNT 

FONMAPCNT 

DISTMAPCNT 

PONBANGCNT 


TYPE 

PURPOSE 

Integer 

Holds  the  first  observa¬ 
tion  number. 

Integer 

Switch  which  selects 
phonet  file. 

Integer 

Number  of  phonets  to  do. 

Integer 

Number  of  distances  to 
compute. 

Integer 

Counts  distances  computed. 

Integer 

Holds  the  first  phonet 
number. 

Integer 

Switch  which  selects  dis¬ 
tance  rule. 

Integer 

Number  of  required  memory 
blocks  for  observations. 

Integer 

Number  of  required  memory 
blocks  for  phonets. 

Integer 

Number  of  extended  memory 
blocks  to  use  for  phonets. 

Integer 

Number  of  extended  memory 
blocks  to  use  for  obser¬ 
vations. 

Integer 

Number  of  extended  memory 
blocks  to  use  for  distances 

Integer 

First  extended  memory 
block  used  for  observa¬ 
tions. 

Integer 

First  extended  memory 
block  used  for  phonets. 

Integer 

Extended  memory  block  used 
for  distances. 

Integer 

Counts  phonets. 

Integer 

Specifies  observation 
points  left. 

OBSLEFT 


VARIABLE 

TYPE 

PURPOSE 

PONLEFT 

Integer 

Specifies  phonet  points 
left. 

NPPTS 

Integer 

Number  of  words  in  phonet. 

WKOBS 

Real  Array 

Array  Processor  memory 
array. 

WKFON 

Real  Array 

Array  Processor  memory 
array. 

DIST 

Integer  Array 

Output  distance  file. 

JDISTNBR 

Integer 

Counts  phonets  in  WKFON. 

JOUTDIST 

Integer 

Counts  distances  in  array 

DIST. 

NDBFDIST 

Integer 

Holds  disk  block  to  write 
distances  to. 

NQBFOBS 

Integer 

Number  of  quarter  blocks 
for  observations. 

NQBFPH 

Integer 

Number  of  quarter  blocks 
for  phonets. 

CNTOBS 

Integer 

First  disk  block  of  obser¬ 
vations  to  get. 

IFBC 

Integer 

First  data  point  in  disk 
block. 

I CNTOBS 

Integer 

Counts  observation  disk 
blocks  read. 

CNTFON 

Integer 

First  disk  block  of  phonets 
to  get. 

IFBD 

Integer 

First  data  point  in  disk 
block. 

I CNTFON 

Integer 

Counts  phonet  disk  blocks 
read. 

NFONLD 

Integer 

Number  of  phonets  in  next 
map. 

NPLTB 


Integer 


Number  of  phonets  left  to 
do. 


I 


VARIABLE 

TYPE 

PURPOSE 

FONLEFT 

Integer 

Number  of  points  in  WKFON 
to  fill. 

JOUTBLKS 

Integer 

Number  of  distance  quarter 
blocks  to  write  to  disk. 

JOUTDIST 

Integer 

Counts  distances. 

FON 

Real  Array 

Mapped  through  extended 
memory  to  g et  phonets. 

OBS 

Real  Array 

Mapped  through  extended 
memory  to  get  observations 

SWITCH  SETTINGS; 

SWITCH  SETTINGS 


OUTD 

FSPFSQ 

FDSTR 


0  for  disk  output  only 
10  for  screen  output  as  well 
12  for  printed  output  as  well 

4  to  run  the  observation  file  against  itself 
3  to  use  separate  phonet  file 

1  for  Ml  distance  rule 

2  for  M2  distance  rule 


RELATED  PROGRAMS: 

PLTA,  PLTO,  FITS,  PLTT,  PLTN,  DRVR,  S64,  S128,  CHOOS, 
FCTR,  GFDB,  CHKO,  CHKJ,  GRPH2. 
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ROUTINE  DSTN  NOTE:  THIS  ROUTINE  IS  FOR  FORTRAN  5  !  ! 

THIS  SUBROUTINE  CALLS  ARRAY  PROCESSOR  ROUTINES  AND  USES  EXTENDED 
NEHORY  (14  BLOCKS  OR  MORE)  TO  COHPUTE  A  SPECIFIED  NUMBER  OF  BEST 
DISTANCES,  BEST  BEING  SMALLEST ,  BETWEEN  SPECIFIED  OBSERVATIONS 
AND  SPECIFIED  PHONETS. 

BY:  CAPT  DAN  MARTIN 
DATE:  9/21/82 
SUBJ :  ACOUSTIC  ANALYSIS 

THIS  ROUTINE  IS  CALLED  BY  URVR  AND  CALLS  ROUTINES  CHOOS,  CHKO,  CHKJ, 
AND  CFDB.  ONE  OF  TUO  DISTANCE  RULES,  Ml  OR  M2,  CAN  BE  SELECTED. 
INPUT  FILES  ARE  SPECTRAL  FILES  OUTPUT  BY  DRSU.  THE  OUTPUT  FILE  IS 
AN  INTEGER  FILE  OF  CONSECUTIVE  UNITS,  EACH  UNIT  BEING  4+2N  POINTS 
LONG,  N  BEING  THE  NUMBER  OF  BEST  PHONET  CHOICES  SPECIFIED.  IN  EACH 
UNIT,  THE  FIRST  ELEMENT  IS  THE  OBSERVATION  NUMBER,  THE  SECOND 
ELEMENT  IS  THE  OBSERVATION  ENERCY,  THE  THIRD  AND  FOURTH  ARE  AH 
ORDERED  PAIR:  (PHONET  NUMBER, MAX  DISTANCE),  AND  THE  REMAINING 
POINTS  ARE  N-ORDERED  PAIRS,  ORDERED  BEST-TQ-WQRST:  (PHONET  NUMBER, 
DISTANCE). 

THIS  ROUTINE  ACCOMPLISHES  THE  FOLLOWING  TASKS: 

1)  GET  OBSERVATION  FILE 

2)  SPECIFY  OBSERVATIONS -FIRST  AND  LAST-TO  DO 

3)  GET  PHONET  FILE 

4)  SPECIFY  PHONETS-FIRST  AND  LAST-TU  DO 

5)  GET  DISTANCE  RULE 

6)  GET  FILE  NAME  FOR  OUTPUT  DISTANCES 

7)  CONSTRUCT  OUTPUT  HEADER 

8)  PARTITION  EXTENDED  MEMORY  BETWEEN  DISTANCES, 
OBSERVATIONS,  AND  PHONETS 

9)  READ  OBSERVATIONS  INTO  MEMORY 

10)  GET  ONE  OBSERVATION 

11)  READ  PHONETS  INTO  MEMORY 

12)  CET  PHONETS  TO  DO 

13)  GET  ALL  DISTANCES  FOR  CUkRENT  PHONET 

14)  CALL  CHOUS  TO  FILL  OUTPUT  FILE 

15)  GET  NEXT  OBSERVATION 

18)  GO  TO  STEP  11  UNTIL  ALL  OBSERVATIONS  ARE  DUNE 
17)  OUTPUT  DISTANCES  ACCUMULATED 
FOR  MORE  INFORMATION,  SEE  THE  USERS  MANUAL  OR  MY  THESIS. 

SUBROUTINE  DSTN 
OVERLAY  ODSTN 
INCLUDE" ARRAYP:F5 APS ,FR“ 

PARAMETER  NDP*1024 
PARAMETER  NDPHl’1023 
PARAMETER  QTRBLI =256 

REAL  WKFON,WXOBS,NPFSR,NPFPR,F ACTOR, FON,OBS,HOLDIST 
REAL  HAXD1ST  ,H1N01ST 

INTECER  FSPFH,PILE2, HEADER, STS  1 , STSL , NST STB , HPS TS ,NPFS 
INTEGER  D1ST ,FSPFSQ ,FILE1 ,PTSi ,PTSL, JQUTDIST 
INTECER  NPTSTB ,HPPTS,NPFP ,FDSTR ,FILE3 ,OUTD ,FLAG ,CB1 ,SKIP 
INTEGER  XMEM ,NBFDIST ,NBFOBS ,NBFPH ,CNTOBS ,CNTFON ,OBSMAPCNT 
INTEGER  FONHAPCNT,FONBANCCHT,OBSLEFf,NQBFOBS,NUBFPH,NGBFBIST 
INTEGER  DISTMAPCNT ,OBSBANGCNT ,FONLEFT ,ICNTOBS , IF8C 
INTEGER  NPLTB, JFONMAPCNT , JOBSMAPCNT ,DISTMAPCNTMAX ,NDBF  VIST 


INTEGER  JD1STNBR ,  JCNTOBS , JCNTFON ,  ICNTFON ,NFONLD ,  INITREH 


COHROH  /  APR  /  UKOBS(HDP) ,UKFON(NDP) 

COHROH  /  VALS  /  OBS(NDP) ,FQN(NDP)  ,DIST(NDP)  ,OBSLEFT ,FSPFN 
COHROH  /  VALS  /  IHITHEH  ,FILE3( 13 > , HE AOER ( QTRDLK )  ,NPFSR,NPFPR 
COHHQN  /  VALS  /  PTS1  ,PTSL,HPTSTB,NPPTS,NPFP,FDSTR,QUTD,NBFOIST 
COHROH  /  VALS  /  NBFOBS,NBFPH,CNTOBS,CNTFON,OBSRAPCNT,FQNHAPCNT 
CORHON  /  VALS  /  FONBANGCNT,NQBFQBS,NQBFPH,NQBFDIST,DISTRAPCNT 
COHROH  /  VALS  /  ICHTOBS, ICNTFON, 1FBC,JF0NHAPCNT,NF0NLD 
COHROH  /  VALS  /  OBSBANGCHT ,FOHLEFT  ,NPLTB ,  JQBSHAPCNT ,NDBFDI$T 
COHROH  /  VALS  /  DISTHAPCNTHAX, JOUTD1ST , JD1STNBR, JCNTOBS, JCNTFOH 
COHROH  /  VALT  /  SKIP, STSL,NSTSTB,NPSTS,NPFS,XHEH, FACTOR, FSPFSO 
COHROH  /  VALT  /  FILE1 (13> ,FILE2( 13) , FLAG ,STS1 ,ISXHEH 
COHROH  /  VALV  /  CB1(0:CBHAX) 

COHROH  /  VALU  /  H0LD1ST<HDP> ,RAXDIST,HIHDIST,HAXD1STL0C,HIH»ISTL0C 
COHROH  /  VALU  /  NCH01CES,NBRPH0NES 

TTPE“»»»»T0U  ARE  HOU  IN  ROUTINE  DSTN . FK»*«»" 

TYPE"«YOU  HAVE" ,ISXHEH,“  IK  BLOCKS  OF  EXT  HERQRY  ACCESSIBLE" 
XHEH*ISXRER 

TYPE"«ENTER  OUTPUT  DEVICE  •  “ 

ACCEPT" (0  =  NQHE/10:TfcRH INAL/ 12  =  PR INTER ) :  ",OUTD 
IFIOUTD.EQ.OtOR.OUTD.EQ. IOiOR .OUTD.Eu. 125COTO  5 
TTPE"ERROR :  INVALID  VALUE  FOR  OUTD:  “ , OUT D 
GOTO  3 

GET  OBSERVATION  FILE  HARE 

ACCEPT"*»EHTER  HAHE  OF  FILE  THAT  HOLDS  OBSERVATION  TIRE-SLICES:  “ 
READ(11 , 1 0 ) F I LE2 ( 1 ) 

FORMAT ( S 1 3 ) 

CALL  OPEN <  4 , F I LE2  ,2  ,IER) 

IF( IER.NE.l)TYPE"ERROR  OH  OPEN  OF  SPECTRAL  FILE,  IER=  "  ,1ER 

GOTO  14 
CONTINUE 

CET  NUMBERS  OF  1ST  4  LAST  TIRE-SLICES  OF  SPECTRUM 
CALL  RDBLK(4,0, HEADER,  1,IER) 

IF(OUTD.EQ.O)GOTO  19 

UR I TE ( OUT  D ,9001 ) (HEADER! 15), 15=1, 13) 

TYPE"IST  SPECTRAL  TIRE-SLICE  »  IS:* ,HEADER(55> 

TYPE’’LAST  SPECTRAL  TIRE-SLICE  »  IS , HEADER!  54) 

ACCEPT"»»«ENTER  1ST  SPECTRAL  TIRE-SLICE  •  TO  BANG:  ",STS1 
I F ( STS 1 .GE. HEADER (55) ) GOTO  22 

TYPE"ERROR:  1ST  TIRE-SLICE  TOO  SHALL!  RUST  BE  GREATER  THAN  OR" 
TYPE"EQUAL  TO" f HEADER  1 55 ) 

GOTO  20 

ACCEPT"»»*EHTER  LAST  SPECTRAL  TIME-SLICE  «  TO  BANG:  ",STSL 
IF (STSl.LE> HEADER (56) (GOTO  23 

TTPE"ERROR:  LAST  TIME-SLICE  TOO  LARGE!  HUST  BE  LESS  THAN  OR" 
TYPE"EQUAL  TO” , HEADER (56) 

GOTO  22 

I F ( STS  1  .  LT  .STSDGOTO  24 

TTPE"ERROR:  1ST  TIRE-SLICE  LARGER  THAN  LAST!" 


c* 

24 


GOTO  20 

*****  GET  «  OF  TS'S  TO  BANG  AND  AMOUNT  OF  MEMORY  REQUIRED 
NSTSTB*STSL-STS1+1  {»  TS'S  TO  BANG 

NPSTS*HEADER(57>»2  ?•  WORDS  PER  TS 

IB0»HEADER(55) 

Ctimi  GET  PHONET  FILE 

TYPE"*DO  YOU  WANT  TO  BANG  THE  SPECTRAL  FILE  AGAINST  ITSELF?" 
ACCEPT"ENTER  (4*YES/3*NQ> s  ",FSPFSG 
IF(FSPFSQ.EQ.4)G0T0  37 
35  ACCEPT"***ENTER  PHONET  FILENAME:  " 

READ (11 >10)F1LE1( 1) 

30  CALL  0PEN<3/FILE1 ,2,IER) 

CALL  CHECX(IER) 

CALL  RDBLl<3,0, HEADER, 1,IER) 

CALL  CHECK ( I ER ) 

IF (OUTD . EQ . 0) GOTO  39 
WRITE COUTD, 9001) (HEADER< 15), 15*1, 13) 

GOTO  39 

Cihm«  GET  THE  1ST  &  LAST  PHONETS 

37  DO  38  Jl*l  ,13 

38  FILE1 ( J1 ) -PILE2( J1 ) 

39  TYPE"1ST  PHONET  TS*  IS:",HEADER(55) 

TYPE"LAST  PHONET  TS*  IS:  “,HEADER(56) 

41  ACCEPT’»**ENTER  1ST  PHONET  TS  TO  DO:  ",PTS1 

IF (PTS1 .CE. HEADER (55) )GOTO  42 

TYPE"ERRORj  1ST  TIME-SLICE  IS  TOO  SHALL!  HUST  BE  GREATER” 
TYPE-THAM  OR  EQUAL  TQ“ f HEADER (55 ) 

GOTO  41 

4?  ACCEPT"*»»ENTER  LAST  PHOHET  TS  TO  DO:  ",PTSL 
IF(PTSL . LE . HEADER( 56 ) )GOTO  43 

TYPE"ERROR:  LAST  TINE-SLICE  TOO  BIG!  HUST  BE  SHALLER  THAN" 
TYPE“OR  EQUAL  TO" f HEADER ( 56 ) 

GOTO  42 

43  IF(PTS1.LE.PTSL)G0T0  44 
TYPE"ERROR:  1ST  TINE -SLICE  LARGER  THAN  LAST!" 

COTO  41 

44  ACCEPT"  ENTER  NUMBER  OF  DISTANCE  CHOICES  TO  GET:-  ",NCHOICES 

C*»***»  GET  NUMBER  OF  PHONET  TS'S  TO  BANC  AND  AMOUNT  OF  MEMORY 
C  THAT  WILL  TAKE 

NPTST8*PTSL-PTS1»1  ;»  TS'S  TO  BANG 

NLDTDsNSTSTB*(4+2*NCH01CES >  ;»  DISTANCES  TO  GET 

1BP*HEADER(55> 

C«h*i»  VERIFY  •  PNTS  IN  PHONET  TIME-SLICES  HATCH  *  PNTS  IN  SPECTRAL  TS'S 
NPPTS ‘HEADER (57 ) *2 
IF (HPPTS . EQ » NPSTS )C(JTO  60 

TYPE"ERROR:  «  PNTS  IN  PHONETIC  TIME-SLICE  MUST  MATCH  »  PNTS" 

|  >  TYPE"IN  SPECTRAL  TIME-SLICE.  AS  SPECIFIED," 

TYPE"*  PNTS  IN  PTS*",NPPTS 
TYPE"*  PNTS  IN  STS*“ , NPSTS 
GOTO  35 
60  CONTINUE 
C*m»m  GET  DISTANCE  RULE 

ACCEPT"*SELECT  DISTANCE  RULE  ( l*Hl/2--H2) :  " ,FDSTR 

C»*»»»»  GET  DISTANCE  FILE 
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GET  NAHE  OF  FILE  TO  STASH  DISTANCES  INTO 
65  ACCEPT“»**ENT£R  NAHE  OF  FILE  TO  RECIEVE  DISTANCES;  " 

READ(11 ,10)FILE3(1) 

70  CALL  0PEH(2,FILE3,2,1ER> 

IF!IER.NE.1)TY PE "ERROR  OH  OPEN  OF  DISTANCE  FILE,  IER=  “,IER 
IF(IER.NE.  1 3 ) GOTO  75 

TTPE"FILE  DOES  NOT  EXIST.  HILL  CREATE  IT  FOR  YOU." 

CALL  CFILU (F1LE3,2 , IER1 ) 

IFIIER1 .NE,1)T YPE “ERROR  ON  FILE  CREATION.  IER1 =  “,IER1 
COTO  70 

75  CONTINUE 

C>*tt««  STORE  FILENAMES  IN  HEADER  OF  DISTANCE  FILE 
DO  100  J I  *  t ,QTR8LK 
100  HEADER! J1>:0.0 

DO  110  JIM, 13 
HEADER! Jt)=FILE3(Jl ) 

J2=J1+13 

HEADER  (J2>=F1LE2(J1> 

J3= J 1 +26 

110  HEADER! J3):FILE1(J1) 

HEADER(40) =STS1 
HEADER! 41 ) =  S  T  S  L 
HEADER ( 42 ) =PTS1 
HEADER! 43) -PTSL 

HEADER! 46)  M'SPFSQ  ;SUI TCH  INDICATES  IF  OHS  AND  PHONS  ARE  SANE 

HEADER! 47 )=NCHOICES  ;•  DISTANCE  CHOICES 

HEADER(48) =NSTSTB  OBS  TIHE  SLICES 

HEADER(49)*NPTSTB  ;»  PHQN  TIHE  SLICES 

HEADER! 50) »NPPTS  ELEHENTS  PER  TIHE-SLICE 

C»**»*»  DIVVY-UP  EXTENDED  HEHURY  BETWEEN  OBSERVATION  SPECTRUM, 

C  PHOHET  SPECTRUH ,  AND  THE  DISTANCE  FILES. 

111  NPFS* 2+2*IFlX( FLOAT i NSTSTB ) *NPSTS/ (2*NDP) )  ;  HEH  BLKS  FUR  OBS  FILE 
NPFP* 2+2*1  FIX! FLOAT (NPTSTB)»NPPTS/(2»NDP) >  ;  HEH  BLKS  FOR  PHON  FILE 

114  CONTINUE 

116  Jl=IFIX((XMEN-3)/2)*2  ;NEED  EVEN  *  HEN  BLKS 

I F ( HP  FP . LT . J1 )NBFPH  =  NPFP 
IF1NPFP.GE. Jl)NBFPH:Jl 
Jl'IFIX! (XHEH-l-NBFPH)/2)*2 
IF (NPFS .LT . J1 )HBFDBS=HPFS 
IFtNPFS.GE.  JDNBFOBS-Jl 
NBFDISTM 

IF(NBFOBS.CT.O)GOTO  116 

TYPE"ERROR :  INSUFFICIENT  EXTENDED  HEHORY  AVAILABLE." 

TYPE"NEED  AT  LEAST  14  IK  BLOCKS," 

CALL  VHEn(XHEH ,  I  HR ) 

TYPE-YOU  HAVE”  ,XHEN ,"  IK  BLOCKS  OF  EXTENDED  HEHURY  AVAILABLE." 
ACC£PT“*HOW  HANY  DO  YOU  WANT?  " ,XHEH 
XHEH=XHEH-9 
COTO  114 
116  CONTINUE 

HEADER (51) =XHEH 

C  *  *  * ( *  *  INITIALIZE  TO  1ST  AVAILABLE  EXT  HEH  BLK 

OBSHAPCNT* V  ; 1ST  EXT  HEH  BLK  USED  FOR  OBSERVATIONS 

FONHAPCNT'QBSHAPCNT  +  NBFQBS  ;  1  ST  EXT  HEN  BLK 
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ROUTINE  DST A  NOTE:  THIS  ROUTINE  IS  FOR  FORTRAN  5  !  ! 

THIS  SUBROUTINE  CALLS  ARRAY  PROCESSOR  ROUTINES  AND  USES  EXTENDED 
MEMORY  <14  MEMORY  BLOCKS  OR  MORE)  TO  COMPUTE  ALL  DISTANCES 
BETWEEN  EACH  OBSERVATION  SPECIFIED  AND  EACH  PHONET  SPECIFIED. 

BY:  CAPT  DAN  MARTIN 

DATE:  9/21/42 

SUbJ :  ACOUSTIC  ANALYSIS 

THIS  ROUTINE  IS  CALLED  BY  DRVR  AND  CALLS  CHXO  AND  CFDB.  INPUT 
FILES  ARE  SPECTRAL  FILES  OUTPUT  BY  DRSQ,  THE  OUTPUT  FILE  IS  AN 
INTECER  FILE  OF  CONSECUTIVE  UNITS.  EACH  UNIT  CONTAINS  DISTANCES 
BETWEEN  AN  OBSERVATION  AND  EACH  PHONET ,  THE  K-TH  ELEMENT  OF  THE 
J-TH  UNIT  IS  THE  DISTANCE  BETWEEN  THE  J-TH  OBSERVATION  AND  THE 
K-TH  PHONET. 

THIS  ROUTINE  ACCOMPLISHES  THE  FOLLOWING  TASKS: 

1)  GET  OBSERVATION  FILE 

2)  SPECIFY  OBSERVAT1QNS-FIRST  AND  LAST-TO  DO 

3)  GET  PHONET  FILE 

4)  SPECIFY  PHONETS-FIRST  AND  LAST-TO  DO 

5)  GET  DISTANCE  RULE 

6)  GET  FILE  NAME  FOR  OUTPUT  DISTANCES 

7)  CONSTRUCT  OUTPUT  HEADER 

8)  PARTITION  EXTENDED  MEMORY  BETUEEN  DISTANCES, 
OBSERVATIONS,  AND  PHONEYS 

9)  READ  OBSERVATIONS  INTO  MEMORY 

10)  GET  ONE  OBSERVATION 

11)  READ  PHQNETS  INTO  MEMORY 

12)  GET  PHONEYS  TO  DU 

13)  GET  A  DISTANCE 

14)  OUTPUT  DISTANCES  WHEN  1024  HAVE  BEEN  ACCUHULATED 

15)  CET  MORE  PHONETS  UNTIL  ALL  ARE  DONE 
18)  CET  NEXT  OBSERVATION 

17)  GO  TO  STEP  11  UNTIL  ALL  OBSERVATIONS  ARE  DUNE 

18)  OUTPUT  DISTANCES  ACCUMULATED 

FOR  MORE  INFORMATION,  SEE  THE  USERS  MANUAL  OR  MY  THESIS. 

SUBROUTINE  DSTA 
OVERLAY  ODSTA 
lNCLUDE'ARRAYP  :F5APS . FR" 

PARAMETER  NDP= 1024 
PARAMETER  NDPH1*1023 
PARAMETER  OTRBLK’256 

REAL  WKFON, WKOBS ,NLDTD , JCTLD , NPFSR ,NPFPR , F ACTOR, FON, UBS 
INTECER  FSPFN, FI LE2, HEADER, STS1 ,STSL ,NSTST6 ,NPSTS ,NPFS 
INTEGER  DIST ,FSPFSQ ,F1LE1  PTS1  ,PTSL, JOUTDI ST 
INTECER  NPTSTB,NPPTS,NPFP,FDSTR,FILE3,0UTD, FLAG, CB1, SKIP 
INTECER  XHEM,NBFDIST,NBFOBS,NBFPH,CNTOBS,CNTFON,OBSHAPCNT 
INTEGER  FONHAPCNT  ,FONBANCCNT ,OBSLEFT , NQBF08S ,NGBFPH  ,NQBFDIST 
INTECER  DISTHAPCNT ,OBSBANGCNT ,FONLEFT , ICNTOBS , IKBC ,GSPCT 
INTEGER  NPLTB ,J FONHAPCNT ,JQBSHAPCNT ,DISTNAPCNTNAX ,NDBFDiST 
INTECER  JDISTN8R,JCNTQBS,JCHTFQN,ICNTFQN,NF0NLD,INITHEN 


COMMON  /  APM  /  WIOBS(NDP) , WK FUN < N D P ) 

COMMON  /  VALS  /  OBS(NDP) ,FON(NDP) ,D1ST(NDP)  ,OBSLEFT  , FSPFN 
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COMIOH  /  VALS  /  I HI THEN , F I LE3 ( 1 3 ) /HEADER ( QTRBLX ) , NPFSR , NPFPR 
COHHON  /  VALS  /  PTS 1 , PTSL / NPTSTB , HPPTS /HPFP #P»STft  ^OUTO #HBFD1 ST 
COHHOH  /  VALS  /  NBFOBS ,NBFPH ,CNTOBS ,CNTFON ,OBSHAPCNT ,FONHAPCNT 
COHHON  /  VALS  /  PONlANCCNT,NOBFOBS,NOBPPH,NOBFDIST,DISTHAPCHT 
COHHON  /  VALS  /  ICNTOBS,ICNTFON,IFBC,JFQNHAPCNT,HFQNLD 
COHHON  /  VALS  /  OBSDANGCNT,FONLEFT,NPLTB,JOBSRAPCNT,NDBFDIST 
COHHON  /  VALS  /  DISTRAPCNTHAX , JOUTDIST , JD1STNBR , JCNTOBS , JCNTFOH 
COHHON  /  VALT  /  SIIP,STSL,NSTSTB,HPST8,NPFS,XRER, FACTOR, FSPFSQ 
COHHON  /  VALT  /  FILE1(13),FILE2(13),FLAC,STS1,ISXHEH 
COHHON  /  VALV  /  CBl(OtCBHAX) 

TYPE“»»»aYOU  ARE  NON  IN  ROUTINE  DSTA.FR»»«*“ 


SO  TYPE”«YOU  HAVE”,ISXHEH,-  IK  BLOCKS  OF  EXT  HEHORY  ACCESSABLE” 
XHEH*ISXHEH 

3  TYPE"»EHTER  OUTPUT  BEV1CE  a” 

ACCEPT” (0«N0N£/10*TERHINAL/12*PRINTER):  “,OUTD 
IF(0UTD.E0.0.0R.QUTD.EQ.10.QR.0UTD.EQ.12)CDT0  5 
TYPE-ERROR:  INVALID  VALUE  FOR  OUTD:  -,OUTD 
GOTO  3 

CHHH  GET  OBSERVATION  FILE  NAHE 

S  ACCEPT“»*EHTER  NAHE  OF  FILE  THAT  HOLDS  OBSERVATIOH  TIRE-SLICES: 

READtll ,  10 )  FILE2  ( 1 ) 

10  FORHAT ( S13) 

12  CALL  0PEN(4,F1LE2,2,1ER) 

IF (IER. HE. 1 ) TYPE “ERROR  OH  OPEH  OF  SPECTRAL  FILE,  IER=  ”,IER 
CALL  CHIO(IER)  ;IS  UNIT  ■  IN  USE? 

COTO  10 

IS  CONTINUE 

GET  NUHBERS  OF  1ST  4  LAST  TIHE-SLICES  OF  SPECTRUH 
10  CALL  RDBLI<4,0, HEADER, 1, IER) 

CALL  CNECK(IER) 

1F10UTD • NE i 0)HRITE( OUTD, 9001) (HEADER (15) ,15*1,13) 

19  TYPE”  1ST  SPECTRAL  TIHE-SLICE  •  IS:",HEADER(55) 

TYPE“LAST  SPECTRAL  TIHE-SLICE  a  IS:” ,HEADER(56> 

20  ACCEPT” ataENTER  1ST  SPECTRAL  TIHE-SLICE  •  TO  BANG:  ”,STS1 
IF (STS1 »GE*HEADER(5S) )COTU  22 

TYPE-ERROR:  1ST  TIHE-SLICE  TOO  SHALLI  HUST  BE  GREATER  THAN  OR” 
TYPE-EOUAL  T0”,HEADER(55) 

GOTO  20 

22  ACCEPT”a«iEHYER  LAST  SPECTRAL  TIHE-SLICE  a  TO  BANG:  ”,STSL 
IF(STSL.LE.HEADER(SO))GOTO  23 

TYPE-ERROR:  LAST  TIHE-SLICE  TOO  LARGE!  RUST  BE  LESS  THAN  OR" 
TYPE-EOUAL  TO” ,HEADER(50) 

GOTO  22 

23  1F1STS1  iLBiSTSDCOTO  24 
TYPE”ERR0Ri  1ST  TIRE-SLICE  LARGER  THAN  LAST! " 

GOTO  20 

aaaaa  GET  •  OF  TS'S  TO  BANG  AND  AHOUHT  OF  HEHORY  REQUIRED 
HST8TD*8TSL-STSW1  {§  TS'S  TO  BANG 

NPSTS*2>NEADER(S7)  \9  SORBS  PER  TS 


C* 

24 


IB0«HIADBR(55) 

£>«•!•■  6BT  PHONET  PUB 

TYPB"«D0  YOU  WANT  TO  BANC  THE  SPECTRAL  PILE  AGAINST  ITSELF?" 
ACCBPT’ENTER  (4«YES/3«N0) :  ",FSPFSQ 
IP ( PSPPSO • BO • 4 ) GOTO  37 
35  ACCEPT" •••ENTER  PHOHET  FILENAME:  " 

REAKll;  lO)FILEl(l) 

34  CALL  0PEH(3,FILE1,2,IER) 

CALL  CHECK  I ER) 

CALL  RDBLK3,0, HEADER, 1,1ER) 

CALL  CHECK ( I ER > 

IF(0UTD .HE <0)HRITE(0UTD,9001> (HEADER (15) ,15*1,13) 

60T0  39 

C»m»«»  GET  THE  1ST  t  LAST  PHOHETS 
37  DO  30  Jl»l,13 

30  FILE1  ( J1 )  SFILE2(«I1  > 

39  TYPE" 1ST  PHOHET  TS«  IS:" ,HEADER(55> 

TYPE"LAST  PHOHET  TSa  IS:  ",HEADER(54) 

41  ACCEPT"«MENTER  1ST  PHONEY  TS  TO  DO:  ",PTS1 
IF (PTS1 .CE. HEADER (55) )COTO  42 

TYPE"ERROR:  1ST  TIHE-SL1CE  IS  TOO  SHALL!  HUST  BE  GREATER* 
TYPE-TNAN  OR  EQUAL  T0",HEADER(55> 

GOTO  41 

42  ACCEPT"«*aENTER  LAST  PHOHET  TS  TO  DO:  “ ,PTSL 
IF(PTSL.LE.HEADER(54) )GOTO  43 

TYPE"ERROR:  LAST  TIHE-SLICE  TOO  BIG!  HUST  BE  SHALLER  THAN" 
TYPE-OR  EQUAL  TO’ , HEADER (54) 

GOTO  42 

43  IFIPTS1  .LE.PTSDGOTO  44 

TYPE"ERROR:  1ST  TIHE-SLICE  LARGER  THAN  LAST! ~ 

GOTO  41 

Caaaaaa  GET  HUHBER  OF  PHONET  TS'S  TO  BANG  AND  AHOUNT  OF  HEHORY 
C  THAT  HILL  TAKE 

44  NPTSTB’PTSL-PTSWl  ;•  TS'S  TO  BANC 

NLDTD«HSTSTB*HPTSTB  *a  DISTANCE  TO  COMPUTE 

JCTLD*0  ; SET  COUNTER  FOR  a  DISTANCES 

IBP*HEADER(55) 

C««a«««  VERIFY  a  PNTS  IN  PHONET  TINE-SLICES  HATCH  a  PNTS  IN  SPECTRAL  TS'S 
NPPTS*2*HEADER(57) 

IF(NPPTS»EO»NPSTS)GOTO  40 

TYPE-ERROR:  a  PNTS  IH  PHONETIC  TIHE-SLICE  HUST  HATCH  •  PNTS" 
TYPE"IN  SPECTRAL  TIHE-SLICE.  AS  SPECIFIED," 

TYPE’a  PNTS  IN  PTS*",NPPTS 
TYPE"*  PNTS  IN  STS*";NP8TS 
GOTO  35 

40  CONTINUE 
Chiih  GET  DISTANCE  RULE 

ACCEPT"«SELECT  DISTANCE  RULE  (1*H1/2*H2):  “,FDSTR 
C»*«»«*  GET  DISTANCE  FILE 

C»««»»»  GET  NAHE  OF  FILE  TO  STASH  DISTANCES  INTO 

45  ACCEPT"*»»EHTER  NAHE  OF  FILE  TO  REC1EVE  DISTANCES:  " 

READdl  ,10>FILE3(1) 

70  CALL  0PEN(2,F1LE3^,1ER) 

IF(IER.NE.1)TYPE"ERR0R  ON  OPEN  OF  DISTANCE  FILE,  I ER  *  “,IER 
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CALL  CHKO(IER) 

IP(ICR*HE»13)COTO  73 

TYPE"FILE  DOBS  NOT  EXIST.  WILL  CREATE  IT  FORYOU." 

CALL  CFILN!FILE3,2,IERl) 

IF(IBRl.HE.l)TYPE“ERROR  ON  FILE  CREATION.  IER1*  “,IERl 
GOTO  70 
73  CONTINUE 

C>*M«t  STORE  FILENANES  IN  HEADER  OF  DISTANCE  FILE 
DO  100  Jl*l jQTRBLX 
100  HEADER! J1)>0.0 

DO  110  J1«1,13 
HEADER! J1)*FILE3!J1) 

J2*JU13 

HEADER  ( J2)*FILE2(J1) 

J3«JU24 

110  HEADER! J3)*FILB1!J1) 

HEADER(40)*STS1 
HEADER! 41 ) <STSL 
HEADER(42>»PTS1 
HEADER ( 43 ) =PTSL 

HEADER! 46 )*FSPFSQ  ?SH1CH  INDICATING  IF  DBS  AND  PHONS  ARE  SANE 

HEADER!46)«NSTSTB  i*  OBS  TIKE  SLICES 

HEADER!49)*NPTSTB  ?•  PHON  TIKE  SLICES 

HEADER!50)*NPPTS  ;«  WORDS  PER  TIKE-SLICE 

!••»»««  DIVVY-UP  EXTENDED  HEHORY  BETWEEN  OBSERVATION  SPECTRUH, 

C  PHONET  SPECTRUH,  AND  THE  DISTANCE  FILES. 

111  NPFS* 2*2«IFIX! FLOAT ! NSTSTB ) oNPSTS/ ! 2*NDP ) )  ;*  HEN  BLKS  FOR  OBS  FILE 

NPFP * 2+2 tIFIX ! FLOAT ( NPTSTB ) *NPPTS/ ! 2*HDP ) )  }  HER  BLXS  FOR  PHOH  FILE 

114  CONTINUE 

116  Jl»IFIX!!XHEH-3>/2)«2  ; NEED  EVEN  ■  HEN  BLXS 

IFINPFP.LT. J1)NBFPH*NPFP 
IF!HPFP.CE» J1)HBFPH*J1 
Jl«IPIX!<XHEH-l-HBFPH>/2>*2 
1F!NPFS.LT. J1)NBF0BS*NPFS 
IF(NPFS.GE.J1)NBF0BS«J1 
NBFD1ST*! 

IF ( NBFOBS . CT . 0 ) GOTO  116 

TYPE- ERROR z  INSUFFICIENT  EXTENDED  KEKOKY  AVAILABLE." 

TYPE“NEED  AT  LEAST  14  IX  BLOCKS," 

CALL  VHEN!XNEN,IER) 

TYPE-YOU  HAVE“,XKEN,“  IX  BLOCKS  OF  EXTENDED  HEHORY  AVAILABLE." 
ACCEPT-iHOW  HANY  DO  YOU  WANT?  “ ,XREH 
XNEH*XHEH-9 
GOTO  114 
116  CONTINUE 

HEADER!51)«XHEH 

Cmmm  INITIALIZE  TO  1ST  AVAILABLE  EXT  HEN  BLI 

OBSRAPCNT*?  {1ST  EXT  NEK  BLK  USED  FOR  OBSERVATIONS 

FONHAPCNT*OBSNAPCNT«NBFOBS  ;1ST  EXT  HEN  BLK 

{USED  FOR  PHONETS 

DISTRAPCNWONHAPCNT+NBFPH  ?1ST  EXT  REK  BLK 

lUSED  FUR  DISTANCES 

FONBANGCNT'O 

OBSLEFT'NDP 


FILE:  CHOOS 

LANGUAGE:  FORTRAN  5 

DATE:  September  21,  1982 

AUTHOR:  D.  Martin 

SUBJECT:  Acoustic  Analysis 

CALLING  SEQUENCE:  CHOOS 

DATE  OF  LAST  REVISION:  September  21,  1982 

PURPOSE: 

Ibis  program  uses  the  Array  Processor  to  choose  the 
N-best  phonet  choices  from  a  file  of  up  to  512  distance 
measures*  The  number  of  choices,  N,  is  selected  in  the 
routine  which  calls  this,  DSTN. 

DESCRIPTION: 

Location:  DP4: BRATCHET 

Size:  CHOOS. FR 

CHOOS. RB 

PROGRAM  USE: 

This  routine  calls  CHKJ  and  is  called  by  DSTN.  It 
uses  the  Array  Processor  to  choose  a  selectable  number  of 
best  distances,  best  being  minimum.  It  fills  the  output 
file  with  the  observation  number,  the  observation  energy, 
the  phonet  number  of  the  maximum  distance,  the  maximum 
distance,  and  consecutive  ordered  pairs  for  the  choices: 
the  first  number  being  the  phonet  number,  the  second  is  the 
distance;  the  first  pair  is  the  best  choice  and  the  N  pair 
is  the  Nth  choice.  Each  file  segment  of  4+2N  contains  these 


3994  bytes 
2870  bytes 


assignments  in  ths  order  given  here.  The  output  file  is 
an  integer  file. 

Just  prior  to  assignment,  the  energy  value  and  each 
distance  is  checked,  while  real,  for  integer  overflow. 

Any  value  larger  than  32767  is  assigned  the  value  32767, 

which  is  the  largest  integer  the  Eclipse  S/250  recognizes. 

» 

The  number  of  choices  is  selected  in  routine  DSTN. 

A  flowchart  of  this  routine  is  in  Figure  75.  When 
called,  the  distances  are  loaded  into  Array  Processor 
memory.  After  the  energy  and  observation  number  are  loaded 
into  the  output  array  DIST,  the  Array  Processor  is  used  to 
first  get  the  input  distance  file  maximum  and  correspond¬ 
ing  phonet  number.  Then  the  routine  enters  a  loop  which 
successively  gets  the  minimum  distance  and  phonet  number, 
replaces  the  minimum  distance  with  the  max,  then  gets  the 
next  minimum.  This  continues  until  the  selected  number  of 
choices  is  obtained. 

LIST  OF  VARIABLES; 


VARIABLE 

TYPE 

PURPOSE 

NPFSR 

Integer 

Energy  in  observation. 

NPTSTB 

Integer 

Number  of  phonets. 

HOLDIST 

Real  Array 

Holds  distances. 

WKOBS 

Real  Array 

Holds  distances  in  APM. 

JOUTDIST 

Integer 

Counts  output  array 
elements. 

* 

OBSBANGCNT 

Integer 

Counts  observations  done 

Load 

distance 
array  into 
APM 


I 

Loab  obser¬ 
vation  num¬ 
ber  and 
ene 


^T^ntoJ 


Get  max 
distance 
land  locatioH 

Load  max- 
distance  and 
location 
into  OIST 


Get  min 
distance 
land  location 

X 


Load  min 
distance  and 
location  in-) 
to  DIST 


kepi ace  min 
distance 
with  max 
distance 


Yes 


No 


FIGURE  75.  Flowchart  of  CHOOS 


VARIABLE 

TOPE 

PURPOSE 

DIST 

Integer 

Output  array. 

MAXDIST 

Integer 

Maximum  distance  in 
HOLDIST. 

MAXDISTLOC 

Integer 

Phonet  number  correspond 
ing  to  MAXDIST. 

MINDIST 

Integer 

Minimum  distance  in 
HOLDIST. 

MINDISTLOC 

Integer 

Phonet  number  correspond¬ 
ing  to  MINDIST. 

SWITCH  SETTINGS;  None. 

RELATED  PROGRAMS; 

PLTO,  PLTA,  PLTS ,  PLTN,  PLTT,  GRPH2,  DRVR,  DSTA,  DSTN, 


S128,  S64,  CHKO,  CHKJ,  PCTR,  GFDB 
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»•■•••  ROUTINE  CH008*  THIS  ROUTINE  IS  FOR  FORTRAN  S  !  ! 

THIS  ROUTINE  USES  THE  ARRAY  PROCESSOR  TO  CHOOSE  THE  N-BEST  PHONETS 
'  FRON  AN  ARRAY  OF  UP  TO  512  DISTANCE  HEASURES.  THE  NURIER  OF 

CHOICES,  M,  IS  SELECTED  IN  THE  ROUTINE  WHICH  CALLS  THIS  ONE,  DSTN. 

DYt  CAPT  DAN  NARTIN 

DATE)  9/21/62 

SUDJt  ACOUSTIC  ANALYSIS 

THIS  ROUTINE  CALL8  CNEJ  AND  IS  CALLED  SY  DSTN.  IT  FILLS  AN  OUTPUT 
FILE  WITH  THE  OBSIRVATION  NURDER,  THE  OBSERVATION  EHERCY ,  THE 
PHONET  NURDER  OF  THE  RAX I HUH  DISTANCE,  THE  RAXIHIHUH  DISTANCE, 

AND  CONSECUTIVE  ORDERED  PAIRS  FOR  EACH  CHOICE:  (PHONET  NURDER, 

NEXT  SHALLEST  DISTANCE).  THE  OUTPUT  FILE  CONSISTS  OF  A  HEADER 
OF  ONE  RISE  BLOCK  AND  CONT1CUOUS  UNITS  OF  4*2N  POINTS  IH  LENGTH, 

H  BEING  THE  NURDER  OF  CHOICES  RADE.  EACH  UNIT  CONTAINS  THE  ABOVE 
LISTED  VALUES  FOR  EACH  OBSERVATION. 

THIS  ROUTINE  ACCOHPL1SHES  THE  FOLLOWING  TASKS: 

1)  LOAD  DISTANCE  ARRAY  INTO  ARRAY  PROCESSOR  REHORY 

2)  LOAD  OBSERVATION  NURDER  AND  ENERGY  INTO  OUTPUT 
ARRAY  BIST 

3)  GET  RAX  DISTANCE  AND  LOCATION 

4)  LOAD  RAX  DISTANCE  AND  LOCATION  INTO  BIST 

5)  GET  HIN  DISTANCE  AND  LOCATION 

6)  LOAD  HIN  DISTANCE  AND  LOCATION  INTO  BIST 

7)  REPLACE  HIN  DISTANCE  WITH  RAX  DISTANCE 
S)  CO  TO  TASK  5  UNTIL  ALL  CHOICES  ARE  HADE 

FOR  NORE  INFORMATION  SEE  THE  USERS  MANUAL  OR  NY  THESIS. 

SUBROUTINE  CHOOS 

OVERLAY  OCHOO 

INCLUDE "ARR A YP:F5 APS. FR’ 

PARAMETER  HDP*1024 
PARAMETER  NDPR1-1023 
PARAMETER  QTRBLK*254 

REAL  WXFON ,HXODS ,HPFSR ,NPFPR , FACTOR ,FON , ODS , HOLD  1ST 
REAL  HAXDlSTfHINDlST 

INTEGER  FSPFN, FILED, HEADER, STS1,STSL,NSTSTB,NPSTS,NPFS 
INTEGER  BIST , FSPFSQ , F I LE 1 ,PTS1 ,PTSL ,JOUTDIST 
INTEGER  NPTSTB,NPPTS,NPFP,FDSTR,FIIB3,0UTD, FLAG, CD1, SKIP 
INTEGER  XREH ,NBFDIST ,MBFQBS ,NDFPH ,CHTQBS ,CHTFON , OISMAPCNT 
INTEGER  FONRAPCNT ,FOHBANCCN? , ODSLEFT ,NQDFODS ,NQDFPH ,NOBFDIST 
INTEGER  DISTHAPCNT ,ODSDANCCNT ,FONLEFT , ICNTODS , IFDC ,GSPCT 
INTEGER  HPLTD,JFONNAPCNT,JQDSRAPCNT,DISTHAPCNTRAX,NDDFDIST 
INTEGER  JDISTNDR , JCNTODS , JCNTFON , ICNTFON ,NFONLD , INI THEM 
INTEGER  HAXDI8TL0C,HINDISTL0C,NCH0ICES,NDRPH0NE8 

CORRON  /  APR  /  WKODS(NDP) ,WKFON(HDP) 

COMMON  /  VALS  /  OBSi'HDP) , FOH ( NDP )  ,DIST ( NDP ) , ODSLEFT ,FSPFH 
CORRON  /  VALS  /  1NITHEH ,F1 LE3 ( 13) /HEADER (QTRDLK ) ,NPFSR,NPFPR 
CORRON  /  VALS  /  PTS1,PTSL,NPTSTB,NPPTS,NPFP,FD8TR,0UTD,HBFDI$T 
CORRON  /  VALS  /  HDFOBS,NDPPH,CMTOBS,CHTFON,OBSRAPCNT, FONRAPCNT 
COMMON  /  VALS  /  PONIANCCHT,NQIPODS,NOIFPH,NQBFDIST, DISTHAPCNT 
COMMON  /  VALS  /  ICNTODS , ICNTFON , IFDC , JFONHAPCNT ,NFQNLD 
CORRON  /  VALS  /  0D8DAN6CNT,F0NLEFT,NPLTD,J0DSHAPCNT , NDBFDI ST 
COMMON  /  VALS  /  D1STMAPCNTMAX, JOUTD1ST , JblSTHDR ,  JCHTOBS , JCHTFOM 
CONNON  /  VALT  l  SI I P , STSL , NSTSTD ,NPSTS , HPFS , XHEtl , FACTOR , FSPFSQ 
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COh«W  /  VALT  /  P1LEK13)  ,FILE2(13),FLA6,STS1 
COMMON  /  VALV  /  CBl(OtCBNAX) 

COItHOH  /  VALU  /  HOLD1ST ( MDP ) ,HAXDIST , HINDIS! ,HAXBISTLOC ,H1NDISTL0C 
COMMON  /  VALU  /  HCN0ICE8 ,NBRPHONES 

SEIP-0 

Ciihh  INITIALIZE  APS  ANB  APN  NAP  FOR  ERBB 
IF ( SKIP • EQ • 0 ) GOTO  502 
SI1P>0 

CALL  APINIT<NIL,DTARAY,M*lfHEH,IER) 

CALL  APHAPt DTARAY ,0 ,4 , IER ) 

502  CONTINUE 

Caa*«»*  LOAB  ARRAY  HOLBIST  INTO  UXOBS 
CALL  APSETL<NPTSTB,IER) 

CALL  CBSET <CB1 ,CBAAHH , HOLD I ST ,CBAXR , UXOBS ,  IER ) 

CALL  VLDR(CBl) 

Ci*»«H  LOAB  OBSERVATION  «  INTO  OUTPUT  ARRAY  BIST 
DIST<J0UTDIST)«0B6BANCCNT-1 

CALL  CHXJ  {SEE  IF  BIST  IS  FULL. 

J0UTDISKJ0UTDISK1 
Ctatatt  LOAB  OBSERVATION  ENERGY  NEXT 

NPFSR-HPFSR/100  {SCALE  ENERGY 

IF (NPFSR.GT. 32767 )NPFST*32767 
DISK  JOUTD  1ST)  *HPFSR 
CALL  CHXJ 

JOUTDISKJOUTDISKl 

C»»»»»»  GET  MAX  BISTANCE  ANB  ITS  LOCATION 

CALL  CBSEKCBl,CBAXR;NXQBS,CBAAMR,HAXDIST,CBAIHH,nAXDISTLDC,IER) 
CALL  VHXR(CBl) 

HAXDlSTLOCsNAXDISTLOC+l 
DISK JOUTD 1ST) *HAXDISTLOC 
CALL  CHXJ 

JOUTDISKJOUTDISKl 
DISK JOUTD 1ST) *MAXD 1ST 
CALL  CHXJ 

JOUTDISKJOUTDISKl 

Caaa«a»  GET  HINIHUN  DISTANCES  AND  THEIR  LOCATIONS,  I.E.,  PHONET  aS 
DO  2000  J2*1,NCH0ICES 

CALL  CBSET(CB1,CBAXR,HK0BS,CBAARN,HINBIST,CBAINN,NINDISTL0C,IER) 
CALL  VHNRICBl)  {GET  HINIHUN 

DISKJOUTDISTKHINDISTLOKl  {LOCATION  OF  HINIHUN 

MKOBS(DISKJOUTDIST)  )*NAXDIST 
CALL  CHXJ 

JOUTDISKJOUTDISKl 
BIST( JOUTB 1ST) * HINDIS! 

CALL  CHXJ 

2000  JOUTB  1ST*  J0UTB1  SKI 

RETURN 

END 
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FILE: 

LANGUAGE: 

AUTHOR: 

DATE: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OF  LAST  REVISION: 


FCTR 

FORTRAN  5 
D.  Martin 

September  21,  1982 
Acoustic  Analysis 
FCTR  (tf,  M,  H,  Q,  A,  B,  C,  D) 
September  21,  1982 


PURPOSE: 

Ttiis  subroutine  is  called  by  GFDB  and  calls  no  sub¬ 
routines.  It  factors  an  integer  into  a  linear  combination 
of  four  terms.  GFDB  uses  it  to  determine  the  unit  (disk 
block,  etc. )  and  starting  data  point  in  that  unit  at  which 
specified  data  reside.  For  example,  if  the  N  element 
of  an  integer  disk  file  is  to  be  accessed,  then  - GFDB  will 
call  FCTR  to  specify  the  disk  block  and  location  in  that 
block  where  that  Nth  element  resides. 


DESCRIPTION: 

Location:  DP4: BRATCHET 

Size:  FCin.FR  605  bytes 

FCTR. RB  218  bytes 


ARGUMENT  STRUCTURE: 


ARGUMENTS 

TYPE 

PURPOSE 

W 

Integer 

Integer  to  be  factored 

M 

Integer 

Factor. 

A 

Integer 

Coefficient  of  M. 

H 

Integer 

Factor. 
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ARGUMENTS 

TYPE 

PURPOSE 

B 

Integer 

Coefficient  of  H. 

Q 

Integer 

Factor# 

C 

Integer 

Coefficient  of  Q. 

D 

Integer 

Coefficient  of  one 

PROGRAM  USE: 

Values  of  W,  M,  H,  and  Q  are  passed  to  FCTR.  FCTR 
returns  values  of  A*  B,  c,  and  D  for  which: 

W  *  AM  +  BH  +  CQ  +  D. 

RELATED  PROGRAMS: 

PLTO,  PLTS ,  PLTA,  PLTN#  PLTT,  DRVR,  DRSQ,  DSTA,  DSTN, 
S64,  S128,  GFDB,  GRPH2,  CHKO,  CHKJ. 
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Chiih  ROUTINE  FCTR.  THIS  ROUTINE  POR  FORTRAN  S  !  ! 

C  THIS  SUBROUTINE  IS  CALLED  BY  CPDB  AND  CALLS  NO  SUBROUTINES.  IT 

C  FACTORS  AN  INTEGER  INTU  A  LINEAR  COHBINATIQH  OF  FOUR  1-ST  ORDER 

C  TERNS.  GFDB  USES  IT  TO  DETERMINE  THE  DISK  BLOCK  AND  STARTING 

C  ELEMENT  IN  A  DISK  FILE  GIVEN  A  SPECIFIED  TIRE  SLICE  OR  SPECTRAL 

C  SLICE  NUMBER. 

C  BY x  CAPT  DAN  MARTIN 

C  DATE;  9/21/52 

C  SUB J s  ACOUSTIC  ANALYSIS 

C  SEE  THE  USERS  MANUAL  OR  MY  THESIS  FOR  MORE  INFORMATION. 

SUBROUTINE  FCTR<M,H,QjA,B,C,D> 

INTEGER  U,H,H,Q,A,B,C,D,X,Y 
Chiih  U*A«H«B»H*CiQ*D 
A*U/H 
X*U-A»M 
B*X/H 
Y*X-B«H 
C«Y/Q 
D»Y-C«0 
RETURN 
END 
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PILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OP  LAST  REVISION 


GFDB 

FORTRAN  5 

September  21,  1982 
D.  Martin 
Acoustic  Analysis 
GFDB  (IA,  IB) 
September  21,  1982 


PURPOSE; 

This  routine  is  called  by  PLTA,  PLTS,  PLTN,  PLTT,  DSTA, 
and  DSTO.  It  calls  FCTR  to  compute  the  location  of  a 
specified  data  point. 


DESCRIPTION: 

Location:  DP  4 : BRATCH2T 

Size:  GFDB. FR  644  bytes 

GFDB. RB  278  bytes 


ARGUMENT  STRUCTURE: 
ARGUMENT  TYPE 

IA  Integer 

IB  Integer 


PURPOSE 

Passes  file  slice  number. 

Returns  first  disk  block. 

Passes  number  of  integer 
words  per  slice. 

Returns  first  data  word. 


PROGRAM  USE: 

The  plot  routines  PLTA,  PLTS,  PLTN,  PLTT,  DSTA,  and 
DSTN  use  this  routine  to  locate  specified  data. 
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ROUTINE  CFDB.  THIS  ROUTINE  IS  FOR  FORTRAN  5  I  ! 

THIS  ROUTINE  6ETS  THE  1ST  OISE  BLOCK  (254  INTEGER 
U0RI8)  AND  THE  1ST  ELEMENT  OF  THE  ARRAY  IIHICH  HAS  LOADED 
AT  THAT  DISK  BLOCK.  THE  FIRST  ELEHENT  OF  THE  ARRAY 
IS  NUNBER  ONE. 

BY  s  CART  DAN  HARTIN 

DATE:  9/21/02 

SUBJ:  ACOUSTIC  ANALYSIS 

HHEN  PASSED:  1A*FILE  SLICE  NUHBER 

IB'NUHBER  OF  INTEGER  WORDS  PER  SLICE 
WHEN  RETURNED:  IA*1ST  DISK  BLOCK 
IB* 1ST  DATA  WORD 

FOR  HORE  INFORMATION,  SEE  THE  USERS  MANUAL  OR  MY  THESIS. 

SUBROUTINE  GFDB(1A,1B) 

INTEGER  H,A,B,C,D 

W*IB»(IA-1> 

CALL  FCTR(U,1024,512,256,A,B,C,D) 

IA«4*A+2»B+C 
1 8  s  D«  1 
RETURN 
END 
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FILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OF  LAST  REVISION: 


CHKJ 

FORTRAN  5 

September  21,  1982 
D.  Martin 
Acoustic  Analysis 
CHKJ 

September  21,  1982 


PURPOSE: 

This  routine  is  called  by  DSTN  and  calls  no  subroutines. 
It  checks  the  number  of  distances  accumulated,  writes  them 
if  1024  distances  are  accumulated,  resets  the  counter,  and 
increments  a  pointer  to  the  next  disk  block  to  write. 


DESCRIPTION: 


Location:  DP4: BRATCHET 

Size:  CHKJ. FR  2324  bytes 

CHKJ. RB  1012  bytes 


LIST  OF  VARIABLES: 
VARIABLE  TYPE 

JOUTDIST  Integer 


DIST 


Integer  Array 


NDBFDIST  Integer 


PURPOSE 

Number  of  distances 
accumulated. 

Holds  distances  accumu¬ 
lated. 

Disk  block  to  write  out 
put  to. 


PROGRAM  USE: 

Ibis  routine  is  designed  for  use  by  DSTN,  to  be  called 
each  time  JOUTDIST  is  incremented. 
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RELATED  PROGRAMS: 

PLTO,  PLTA,  PLTS,  PLTN,  PLTT,  PCTR,  GRPH2,  CHOOS,  DRVR, 
DRSQ,  S128,  S64,  DSTA,  DSTN,  CKHJ,  CHKO. 
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ROUTINE  CNXJ.  THIS  ROUTINE  FOR  FORTRAN  5  I  ! 

THIS  ROUTINE  IS  CAILEI  BY  BSTN  ANB  CALLS  NO  SUBROUTINES.  IT 
CNECES  THE  NUMBER  OF  I1STAHCES  ACCUMULATES ,  NRITES  THEM  IF 
1024  HAVE  BEEN  ACCUMULATED ,  RESETS  THE  COUNTER,  AND  INCREMENTS 
A  FDINTER  TO  TNI  NEXT  D1SI  BLOCK  TO  URITE  TO. 

BYt  CAPT  DAN  MARTIN 

DATIi  9/21/12 

SUBJi  ACOUSTIC  ANALYSIS 

THIS  ROUTINE  IS  DESICNEI  FOR  USE  BY  DSTN  TO  BE  CALLED  EACH  TINE 
JOUTBIST  IS  INCREMENTED.  FOR  MORE  INFORMATION,  SEE  THE 
USERS  MANUAL  FOR  DSTN  OR  THIS  ROUTINE,  OR  SEE  MY  THESIS. 


SUBROUTINE  CHXJ 
OVERLAY  OCNXJ 
PARAMETER  NDP‘1024 
PARAMETER  NDPH1>1023 
PARAMETER  QTRBLK*236 

REAL  HXFON,HXOB8,NPFSR,NPFPR, FACTOR, FON,OBS,HOLDIST 
REAL  HAXDIST,M1ND1ST 

INTECER  FSPFN,FILE2, HEADER, STS1,STSL,NSTSTR,NPSTS,NPFS 
INTECER  DIST ,FSPFSQ,FILE1 ,PTS1 ,PTSL, JOUTBIST 
INTECER  NPTSTB,NPPTS,NPFP,FDSTR,FILE3,0UTD,FLAC ,CB1  ,SKIP 
INTEGER  XHEH ,NBFDIST ,HBFOBS ,NBFPH,CNTOBS ,CNTFON ,OBSHAPCHT 
INTEGER  FONHAPCNT ,FOHBAHGCHT ,DBSLEFT ,NQBFOBS , NQBFPH , NQBFO 1 ST 
INTECER  DISTHAPCNT ,OBSBANGCNT ,FONLEFT , ICHTOBS , IFBC,GSPCT 
INTECER  HPLTB,JFOHNAPCNT ,JOBSHAPCNT,DISTHAPCNTHAX,HDBFDIST 
INTECER  JDISTNBR , JCNTOBS , JCNTFON , ICNTFOH ,NFOHLD , IH1TMEH 
INTECER  MAXDISTLOC,HINDISTLOC,NCHOICES,HBRPHOHES 

COMMON  /  APH  /  UXOBS(NDP) ,UXFON(NDP) 

COMMON  /  VALS  /  OBS ( NDP >, FON ( HOP ), DISK NDP) ,OBSLEFT,FSPFN 
COMMON  /  VALS  /  INITHEH >F ILE3 ( 13 ) ,HEADER(QTRBLX) ,NPFSR,HPFPR 
COMMON  /  VALS  /  PTS1,PTSL,NPTSTB,KPPTS,NPFP,FDSTR,0UTD,NBFDIST 
COMMON  /  VALS  /  NBFOBS ,NBFPH ,CHTOBS ,CHTFOH ,DBSMAPCNT , FONHAPCNT 
COMMON  /  VALS  /  FONBAMCCNT,NOBFQBS,NBBFPM,NQBFDIST, DISTHAPCNT 
COMMON  /  VALS  /  ICHTOBS, 1CNTFON,IFBC,JFONHAPCNT,NFONLD 
COMMON  /  VALS  /  OBSBANCCNT , FOMLEFT ,HPLTB ,JQBSRAPCHT  ,HDBFDIST 
COMMON  /  VALS  /  B1STMAPCNTHAX, JOUTBIST, JDISTNBR, JCNTOBS, JCNTFON 
COMMON  /  VALT  /  SXIP,STSL,NSTSTB,HPSTS,NPFS,XMEM, FACTOR, FSPFSO 
COMMON  /  VALT  /  FILEK13)  ,FILE2(  13)  ,FLAC,STS1 
COMMON  /  VALU  /  HOLDISKNDP),HAXDIST,HINDIST,NAXDISTLOC,NINDISTLOC 
COMMON  /  VALU  /  NCH01CES ,NBRPHQNES 

IF( JOUTBIST .LE.HDP) GOTO  365  {BIST  ARRAY  HOT  FULL  YET 
JOUTDISTM  {NEED  URITE  TO  BISK  IF  FULL 

IF(0UTB.NE.0)HRITE(OUTB,90O0)(BIST( J2) ,J2*1 ,NBP) 

Cimii  URITE  4  QTR  BLKS  OF  DISTANCES 

CALL  EURB(2,NDBFDIST,36+N0BF0BB«N0BFPH,4,IER) 

IF ( IER . NE •  i ) TYPE1* ERROR  ON  EURB  OF  DISTANCES!  I8R«",1ER 
HDBFDIST»NDIFBIST+4 
365  CONTINUE 

9000  F0RRAKSC12.4) 

RETURN 
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FILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OF  LAST  REVISION 


PL  TO 

FORTRAN  5 

September  21,  1982 
D.  Martin 
Graphics 
PLTO 

September  21,  1982 


PURPOSE: 

This  routine  calls  one  of  four  plot  routines  selected 
from  a  displayed  menu.  The  graph  is  displayed  on  a  Tektronix 
4010-1  Graphics  Terminal  using  the  routine  GRPH2  by  G.  Shaw 
as  modified  by  L.  Kizer  and  D.  Zambon. 

DESCRIPTION: 

Location:  DP4: BRATCHET 

Size:  PLTO. FR  1491  bytes 

PLTO.RB  1650  bytes 


PROGRAM  USE: 

This  program  is  called  by  DRVR  and  calls  PLTS,  PLTA, 
PLTN,  and  PLTT.  When  called,  this  routine  presents  a  menu 
from  which  the  operator  is  to  choose  one  of  five  options, 
one  being  to  return  to  DRVR.  If  the  option  to  plot  spectrum 
is  chosen,  PLTS  is  called  to  plot  the  spectrum  of  a  time 
slice  from  a  disk  file  computed  by  DRSQ.  If  the  option  to 
plot  all  distances  is  selected,  PLTA  is  called  to  plot  dis¬ 
tances  between  one  observation  and  selected  phonets  from  a 
disk  file  computed  by  DSTA.  If  the  option  to  plot  N-best 

distances  is  selected,  PLTN  is  called  to  plot  observation 


194 


energy  and  a el acted  distances  from  a  disk  file  computed  by 
DSTO.  Finally*  if  the  option  to  plot  an  integer  file  is 
selected,  then  one  unit  of  up  to  512  elements  from  any 
integer  disk  file  is  plotted. 

PURPOSE 

Currently  unused,  but  can 
be  used  as  unit  number  for 
screen  output. 

Switch  to  select  plot 
option. 

SWITCH  SETTINGS: 

SETTINGS 

6  to  call  PLTS  to  plot  spectrum 

7  to  call  PLTA  to  plot  distances 

8  to  call  PLTN  to  plot  N-best  distances, 
worst  distance,  and  observation  energy 

9  to  call  PLTT  to  plot  a  segment  of  an 
integer  file 

10  to  return  to  DRVR 

RELATED  PROGRAMS: 

PLTA,  PLTS,  PLTN,  PLTT,  DRVR,  DSTA,  DSTN.  DRSO,  CHKJ, 
CHKO,  CHOOS,  FCTR,  GFDB,  S128,  S64,  GRPH2. 


SWITCH 

SKIP 


LIST  OF  VARIABLES: 
VARIABLE  TVPE 

outd  Integer 

SKIP  Integer 
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ROUTINE  PLTO.  THIS  ROUTINE  IS  FOR  FORTRAN  S  i  ! 

THIS  ROUTINE  CALLS  PLTA,  PLTS,  PLTN,  AND  PLTT  ,  AND  IS  CALLED 
IT  DRVR.  IT  SELECTS  THE  APPROPRIATE  PLOT  ROUTINE  TO  PLOT  THE 
DATA  AS  REQUESTED. 

BY «  CAPT  DAN  HARTIN 
DATE)  9/21/02 
SUBJx  ACOUSTIC  ANALYSIS 

FOR  HORS  INFORMATION  SEE  THE  USERS  MANUAL  OR  MY  THESIS. 

SUIROUTINE  PLTO 
OVERLAY  OPLTO 
INCLUDE  "ARRAYP*F5APS.FR" 

PARAMETER  HNDP«1024 
PARAMETER  0TRBLX>256 
REAL  DATA,XHIN,XMAX 

INTEGER  GSPCT,FILE2, HEADER, FSPFN,OUTD,OBSNBR,PHONl,PHQHl 
INTEGER  DIST ,FONHDR,FRSTDP,FRSTDBLK,ICNT ,D1,IER,SKIP 

COMMON  /  VALS  /  DATA<HNDP),DIST(1024),OBSNBR,FSPFH,OUTD,PHON1 
COMMON  /  VALS  /  PHONL , 1CNT , IER ,FONNBR ,FRSTDP ,FRSTDBLX ,FILE2 ( 13 ) 
COHHON  /  VALS  /  HE AD ER ( OTRBLX > , SKIP 
COMMON  /  VALT  /  Dl< 37) ,ISXHEH 

TYPE“»*»*YOU  ARE  NON  IN  ROUTINE  PLTO.FR****" 

OUTD'IO 

40  TYPE"**NHAT  HILL  YOU  HAVE  HE  PLOT?" 

TYPE"»**»CHOOSE  OPTION****" 

TYPE"*" 

TYPE" A  TO  PLOT  SPECTRUM." 

TYPE"7  TO  PLOT  ALL  DISTANCES, " 

TYPE" A  TO  PLOT  H-IEST  DISTANCES." 

TYPE”?  TO  PLOT  INTEGER  FILE." 

TYPE" 10  TO  RETURN  TO  THE  PROGRAM  WHICH  CALLED  THIS  ONE." 
ACCEPT'aaENTER  DESIRED  OPTION)  “ , SUP 

100  IF(SXIP.NE.*)COTO  110 

CALL  PLTS 
GO  TO  40 

110  IF4SKIP.NE. 7 )GOTO  120 

CALL  PLTA 
GOTO  40 

120  IF(SE1P.NE.S)G0T0  130 

CALL  PLTN 
GOTO  40 

130  IFtSIlP.NE.TlGOTO  140 

CALL  PLTT 
GOTO  40 

140  IPCSXIP.EQ. 10DG0T0  1000 

TYPE"ERROR)  COULDN'T  FIND  YOUR  OPTION." 

TYPE"YOU  SELECTED  OPTION)  ",SKIP 
GOTO  40 

1000  CALL  RESET 

CALL  OVEXIT(OPLTO,IER) 
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Em 


GALL  CHECK (1BR) 

RETURN 

ENI 
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PILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OP  LAST  REVISION 


PLTA 

FORTRAN  5 

September  21,  1982 
D.  Martin 
Graphics 
PLTA 

September  21,  1982 


PURPOSE: 

Ibis  routine  plots  distances  from  a  disk  file  computed 
by  DSTA.  Distances  between  an  observation  and  all  phonets 
specified  by  the  operator  are  displayed  on  a  Tektronix 
4010-1  Graphics  terminal  using  the  plot  routine  GRPH2  by 
G.  Shaw  as  modified  by  L.  Kizer  and  D.  Zambon. 


DESCRIPTION: 

Locat ion :  DP4 : BRATCH  ET 

Size:  PLTA.FR  4371  bytes 

PLTA.RB  G712  bytes 


PROGRAM  USE: 

Ihis  program  is  called  by  PLTO  and  calls  GFDB  and 
GRPH2.  It  plots  distances  from  a  disk  file  between  a 
specified  observation  and  the  first  and  last  phonets 
specified. 

The  operator  is  first  prompted  for  an  option  to  print 
numerical  values  of  data  to  be  plotted.  Ttie  data  values 
plotted  can  be  displayed  on  the  operator's  terminal,  written 
to  the  line  printer,  or  not  displayed  at  all.  Next,  the 
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operator  is  prompted  for  the  name  of  the  disk  file  holding 
the  data  to  be  plotted.  The  file  header  is  read  and 
identifying  information  from  it  displayed  for  the  operator. 
Then  the  desired  observation  number,  and  the  numbers  of  the 
first  and  last  phonet  to  be  plotted,  are  requested.  When 
signalled  to  proceed,  the  program  will  display  the  plot. 

LIST  OF  VARIABLES; 


VARIABLE 

TYPE 

PURPOSE 

OUT 

Integer 

Switch  for  data  outpht 
option. 

OUTD 

Integer 

Channel  number. 

ETLE3 

Integer  Array 

Holds  disk  file  name. 

OBSNBR 

Integer 

Observation  number. 

PHONl 

Integer 

First  phonet. 

PHONL 

Integer 

Last  phonet. 

PONNBR 

Integer 

Number  of  distances  to 
plot. 

FRSTDBLK 

Integer 

First  disk  block  to  read 

PRSTDP 

Integer 

First  data  point  in  disk 
block. 

DATA 

Real  Array 

Holds  data  to  plot. 

DIST 

Integer  Array 

Data  read  from  disk. 

SWITCH  SETTINGS: 

SWITCH 

SETTING 

OUT  » 

10  for  screen  output  of  plotted  data 

» 

12  for  line  printer 

output 

0  for  no  output  of 

numbers . 

OUT 


-  =  ="l 

RELATED  PROGRAMS: 

PLTO,  PLTS,  PLTN,  PLTT,  DRSQ,  S128,  S64,  DSTA,  DSTN, 

CHOOS,  CHKJ,  CHKO,  PCTR,  GFDB,  DRVR,  GRPH2. 
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ROUTINE  PLTA.  THIS  ROUTINE  IS  FOR  FORTRAN  S  !  ! 

THIS  ROUTINE  IS  CALLED  BY  PLTO  AND  CALLS  THE  PLOT  ROUTINES  GRPH2 
AND  CFBD.  CRPH2  HAS  WRITTEN  BY  C.SHAU  AND  HODIFIED  BY  L.KIZER 
AND  D.ZAHBON.  IT  IS  IN  THE  SPEECH  LAB  PROCRAN  LIBRARY. 

THIS  ROUTINE  PLOTS  DISTANCES  COMPUTED  BY  DSTA.  IT  HILL  PLOT 
DISTANCES  BETHEEN  A  SPECIFIED  OBSERVATION  AND  SPECIFIED  PHONETS. 
THE  USER  IS  PROMPTED  FOR  THIS  INFORMATION  AS  NELL  AS  FOR  THE 
FILE  NAME. 

BY:  CAPT  DAN  MARTIN 

DATE:  9/21/&2 

SUBJ:  ACOUSTIC  ANALYSIS 

PLTA  ACCOMPLISHES  THE  FOLLOWING  TASKS: 

1)  GET  DISTANCE  FILE 

2)  DISPLAY  INFORMATION  FROM  HEADER  TO  IDENTIFY  THE  FILE 

3)  GET  OBSERVATION  NUMBER  AND  PHONET  NUMBERS 

4)  LOCATE  DATA  IN  DISK  FILE 

5)  PREPARE  DATA  FOR  PLOT  ROUTINE  GRPH2 
it  CALL  THE  PLOT  ROUTINE 

FOR  MORE  INFORMATION  SEE  THE  USERS  MANUAL  OR  MY  THESIS. 

SUBROUTINE  PLTA 
OVERLAY  OPLTA 
INCLUDE" ARRAYP:F5APS.FR“ 

PARAMETER  HNDP*512 
PARAMETER  QTRBLK:254 
REAL  DATA,XMIN,XMAX 

INTEGER  GSPCT,FILE3, HEADER, FSPFN,QUTD,OBSNBR,PHONl,PHOHL 
INTEGER  DI ST , FONNBR , FRSTDP ,FRSTDBLl , ICNT ,D l , I ER ,OUT 

COMMON  /  VALS  /  DATA < HNDP > ,D 1ST ( 1024)  ,0BSNBR,FSPFN,0UTD,PH0N1 
COMMON  /  VALS  /  PHONL, ICNT, IER, FONNBR, FRSTDP, FRSTDBLK ,F I LE3 ( 1 3 ) 
COMMON  /  VALS  /  HEADER(QTRBLK) 

COMMON  /  VALT  /  Dl(37) 


TYPE"»*»»YOU  ARE  NON  IN  ROUTINE  PLTD.FR»*««" 

30  TYPE“DO  YOU  WANT  OUTPUT  OF  DATA  PLOTTED?" 

TYPE“ENTER:  10  FOR  SCREEN  OUPUT  OF  NUMBERS" 

TYPE*  12  FOR  PRINTED  NUMERICAL  OUTPUT" 

ACCEPT"  0  FOR  NO  NUMERICAL  OUTPUT:  ",OUT 

IF (OUT .EQ .0. OR. OUT ,EQ. 10 .OR.QUT.EQ. 12) GOTO  40 
TYPE"ERROR:  ILLEGAL  INPUT." 

GOTO  30 
40  OUTO 3 1 0 

CiHHi  CET  DISTANCE  FILE  NAME 

SO  ACCEPT"«»ENTER  NAME  OF  FILE  HOLDING  DISTANCES:  " 

READ (11 , 100 ) FILE3 ( 1 ) 

100  FORMAT ( SI 3 ) 

120  CALL  CLOSE! 2 ^ IER ) 

CALL  0PEN(2,FItE3,2,IER) 

CALL  CHECK  IER) 

IF (IER. NE.l )TYPE“ ERROR  ON  OPEN  OF  SPECTRAL  FILE,  IER3  ",IER 
GOTO  ISO 
ISO  CONTINUE 
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T**’ '  :  :  *fr.%  'flip* ,  ••• 


CET  HUHBERS  OF  1ST  4  LAST  TINE-SLICES  OF  SFECTRUH 
ISO  CALL  RDtiX(2,0, HEADER, 1 , IER) 

CALL  CHECK  IER) 

WRITE(0UTD/9110) 

WRITE (OUTD/9 100) (HEADER! 15) ,15*1,45) 

WRITE (0UTD/9000H HEADER (15) ,15=1,13) 

WRITE(0UTD;9001)( HEADER! 15)  ,15*14,24) 

WRITE < OUTD ( 9004) HEADER! 40) ^ HEADER ( 4 1 > 

WRITE (OUTD, 9002) (HEADER( 15), 15*27, 39) 

WRITE(0UTD/9005) HEADER (42) /HEADER (43) 

ACCEPT** I S  THIS  THE  FILE  YOU'RE  AFTER?  (1*YES/0*N0> I  ",FSPFN 
IF(FSPFN.EQ.0)G0T0  50 

190  ACCEPT"ENTER  THE  OBSERVATION  HUNBER:  "/OBSHBR 

IF(OBSHBR,CE.HEADER(40).AND.OBSHBR.LE.HEADER(41))COTO  200 
TYPE**ERRORs  THIS  NUHBER  HUST  LIE  IH  THE  RANCE" 

TYPE" ["/HEADER (40) /" /"/HEADER (41)  /“]" 

GOTO  190 

200  ACCEPT'ENTER  THE  FIRST  PHOHET  HUHBER:  “,PH0N1 

IF (PH0N1.GE. HEADER (42) . AND.PH0N1 ■ LE » HEADER ( 43 ) )GOTO  220 
TYPE“ERR0R:  THIS  NUHBER  HUST  LIE  IH  THE  RANGE” 

TYPE  "C" /HEADER (42) /"," /HEADER (43)  /  ”]” 

GOTO  200 

220  ACCEPT'ENTER  THE  LAST  PHOHET  NUHBER:  ”/PHOHL 

IF(PH0NL.GE.PH0N1 . AND. PHQNL.LE. HEADER (43) (GOTO  240 
TYPE"ERROR:  THIS  HUHBER  HUST  BE  IN  THE  RANGE” 

TYPE“C" /PH0N1 /"/" /HEADER! 43) 

GOTO  200 

240  FONNBR  *PHQNL -PHQN 1*1 

IF(FONNBR*LE*HNDP)COTO  240 

TYPE”ERROR:  YOU  REQUESTED  HORE  THAN  512  DISTANCES.” 

GOTO  200 
240  CONTINUE 

FONHBR*HEADER (43) -HEADER (42) +  1 
C*«*»»*  CET  DISK  BLOCK  AND  1ST  DATA  POINT 
IA*OBSHBR -HEADER (40) tl 
IB«FOHNBR 
CALL  GFDB(IA/IB) 

FRSTDBLK*IA+1 

FRSTDP*IB 

CALL  RDBLX(2/FRSTDBLK,DIST,3/ICHT/IER) 

IF(1ER.NE.9)C0T0  300 

TYPE”READ  EOF!  SUCCESSFULLY  TRANSFERRED" /ICNT/”  QTR  BLKS” 
TYPE'PROCEEDINC  WITH  "/ICNT /"QTR  BLOCKS  TRANSFERRED." 

IER*1 

300  CALL  CHECX(IER) 

CALL  CLOSE ( 2 /IER > 

CALL  CHECK! IER) 

IF(OUT.EQ.O) GOTO  310 
IF ( OUT • EQ .12) OUTDs 12 

Ciiim  HOVE  THE  DATA  UP-FROHT  AND  INTO  THE  ARRAY  TO  PLOT 
310  IL«PH0N1-HEADER(42>+1 

IU*PH0NL-HEADER(42)+1 
FRSTDP*FRSTDPfIL-l 
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10  320  Jl«l,IU-IL4l 
DATA ( J I > *DIST(Jl4FR8TDP-l ) 

IF(0UT. BO. 0)6010  322 
NR1TE(0UTD, 9100) (DATA (15) ,I5*1,IU-IL+1) 

322  ACCEPT'OE  10  CONTINUE?  (EN1ER  ANY  NUMBER)  ",15 

IPdS.NE. 0)6010  324 
324  CONUNUE 

F0NHBR«PH0HL*PH0Nl4l  . 

COLL  6RPH2("  BIS10NCE  PLOT", 1, DATA, U,F0NNBR,0,XMIN,XMAX,0) 
ACCEPTS!  10  CONUNUE?  (ENTER  ANY  NUMBER)  “ ,15 

iras.Eo. 0)6010  sso 

350  CONUNUE 

9100  FORMA! (5612.4) 

9110  P0RN0K15, "HERE'S  THE  FIRST  45  ELEMENTS  OP  THE  HEADER:") 

9000  P0RM0K13, "DISTANCE  FILE  NOME  IS:  ",13S2) 

9001  P0RHAKT5, "OBSERVATION  FILE  NONE  HAS:  ",13S2) 

9002  F0RH0K15, "PHONE!  FILE  NOME  NAS:  ",13S2) 

9004  F0RHAKT2, ‘FIRST  OBSERVATION  TS«  *  ",I5,"  AND  LAST  *  ",I4> 

9005  F0RH0KT2, "FIRST  PHONE!  TS«  *  ",I5,"  AND  LAST  *  ",I6> 

RETURN 

END 
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PILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OF  LAST  REVISION 


PLTS 

FORTRAN  5 

September  21,  1982 
D.  Martin 
Graphics 
PLTS 

September  21,  1982 


PURPOSE: 

This  routine  plots  spectrum  from  a  dish  file  computed 
by  DRSQ.  The  plot  is  displayed  on  a  Tektronix  4010-1  Graphics 
Terminal  using  the  routine  GRPH2  by  G.  Shaw  as  modified  by 
L.  Kizer  and  D.  Zambon. 


DESCRIPTION: 

Location:  DP4 : BRATCHET 

Size:  PLTS. PR  3785  bytes 

PLTS.RB  5302  bytes 


PROGRAM  USE: 

This  program  is  called  by  PLTO  and  calls  GFDB  and  GRPH2. 
It  plots  spectrum  from  a  disk  file  computed  by  DRSQ. 

The  operator  is  first  prompted  for  an  option  to  print 
numerical  values  of  data  to  be  plotted.  The  data  values 
plotted  can  be  displayed  on  the  operator's  terminal,  written 
to  the  lire  printer,  or  not  displayed  at  all.  Next,  the 
operator  is  prompted  for  the  name  of  the  disk  file  holding 
the  data  to  be  plotted.  The  file  header  is  read  and 
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identifying  information  from  it  is  displayed.  The  time 
slice  number  corresponding  to  the  desired  spectrum  is 
requested  and  the  dB  or  magnitude  option  prompted.  When 
the  program  is  signalled  by  the  operator  to  proceed,  it 
will  display  the  plot.  The  value  of  the  first  component 
will  be  the  observation  energy. 

LIST  OF  VARIABLES ; 


VARIABLE 

TYPE 

PURPOSE 

OUT 

Integer 

Switch  for  data  output 
option. 

OUTD 

Integer 

Channel  number. 

PILE2 

Integer  Array 

Holds  disk  file  name. 

OBSNBR 

Integer 

Time  slice  number. 

PONNBR 

Integer 

Number  of  spectral  com¬ 
ponents. 

FRSTDBLK 

Integer 

First  disk  block  to  read. 

FRSTDP 

Integer 

First  data  point  in  disk 
block. 

HEADSR(  57 ) 

Integer 

Number  of  spectral  com¬ 
ponents. 

DATA 

Real  Array 

Data  to  be  plotted  - 
passed  to  GRPH2. 

DIST 

Real  Array 

Data  read  from  disk. 

I  DP 

Integer 

Switch  for  dB  or  magnitude 
plot. 

SWITCH  SETTINGS; 

SWITCH  SETTING 

10  for  screen  output  of  plotted  data 
12  for  line  printer  output 
0  for  no  output  of  numbers 


OUT 


-wynli.i 


SWITCH  SETTING 

I OP  =  1  for  dB  versus  component  number 

»  0  for  magnitude  versus  component  number 

RELATED  PROGRAMS; 

PLTA,  PLTO,  PLTN,  PLTT,  DRVR,  DRSQ,  DSTA.  DSTN,  CHKJ, 
CHKO,  CHOOS,  FCTR,  GFDB,  S128,  S64,  GRPH2. 
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ROUTINE  PITS.  THIS  ROUTINE  18  FOR  FORTRAN  S  I  ! 

THIS  ROUTINE  18  CALLED  BY  PLTO  ANO  CALLS  THE  PLOT  ROUTINES  6RPN2 
ANO  CF1D.  6RPH2  HAS  HRITTEN  IT  C.SHAH  ANI  ROIIFIES  SY  L.KIZER 
ANO  D.ZANBON.  IT  18  IN  THE  SPEECH  LAS  PR06RAN  LIBRARY. 

THIS  ROUTINE  PLOTS  8PBCTRUH  OF  TIRE  SLICES  CONFUTE!  SY  SRSO. 

THE  USER  IS  PRONPTES  FOR  THE  TINE  SLICE  HUMBER,  FILE  NAME, 

ANB  OTHER  NECESSARY  INFORNATIOH. 

BY*  CAPT  SAN  MARTIN 
BATE*  9/21/82 
SUBJs  ACOUSTIC  ANALYSIS 
PLTA  ACCONPLISNES  THE  POLLOUINC  TASES* 

1)  CET  SPECTRAL  FILE 

2)  SI6PLAY  INFORMATION  FROM  HEABER  TO  IBENTIFY  THE  FILE 

3)  CET  OBSERVATION  NUMBER  (TIRE  SLICE  NUMBER) 

8)  LOCATE  DATA  IN  BISK  FILE 
5)  PREPARE  BATA  FOR  PLOT  ROUTINE  CRPH2 
8)  CALL  THE  PLOT  ROUTINE 

FOR  MORE  INFORMATION  SEE  THE  USERS  MANUAL  OR  NY  THESIS. 

SUBROUTINE  PLTS 
OVERLAY  OPLTS 
INCLUDE "ARRAYPjF5APS.fr" 

PARAMETER  HNDP»512 
PARAMETER  0TRBLK*258 
REAL  BATA ,XH1N ,XHAX  ,DIST 

INTEGER  GSPCT,FILE2, HEADER, FSPFN, OUTS, OBSNBR,PHONl,PHONL 
INTEGER  FONNBR,FRSTDP ,FRSTDBLX,ICNT ,D1 , IER,OUT 

COMMON  /  VALS  /  DATA(HHBP) , BIST (5 12) , OBSNBR, FSPFN , OUT D ,PHONl 
COMMON  /  VALS  /  PHONL , ICHT , IER/FONNBR ,FRSTDP ,FRSTDBLX ,FILE2 ( 13) 
COMMON  /  VALS  /  HEABER (OTRBLI) 

COMMON  /  VALT  /  Bl(37) 

TYPE“«»*»Y0U  are  non  IN  ROUTINE  PITS. from" 

30  TYPE"DO  YOU  UANT  OUTPUT  OF  DATA  PLOTED?" 

TYPE"EHTER*  10  FOR  SCREEN  OUTPUT  OF  NUMBERS" 

TYPE"  12  FOR  PRINTED  NUMERICAL  OUTPUT" 

ACCEPT"  0  FOR  NO  NUMBERS:  ",OUT 

lF(OUT.EQ.O.OR.OUT.EQilO.OR.OUT.EQ.12>GOTO  40 
TYPE'ERROR*  ILLEGAL  INPUT." 

GOTO  30 
40  0UT0*10 

GET  SPECTRAL  FILE  NAME 

50  ACCEPT"»»ENTER  NAME  OF  FILE  HOLDING  SPECTRUM:  " 

READS  Ilf 100) FILE2C i ) 

100  FORMAT ( 813 ) 

120  CALL  CL0SE(4,1ER) 

CALL  0PEN(4,FILE2,2,IER) 

CALL  CHECK (1ER) 

IF( IER.Hk. 1 )TYPE"ERROR  ON  OPEN  OF  SPECTRAL  FILE,  IER*  ",IER 
GOTO  180 
130  CONTINUE 

Cihim  GET  NUMBERS  OF  1ST  (  LAST  TIME-SLICES  OF  SPECTRUM 
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110  CALL  RDDLK 4,0, HEADER, 1,IIR> 

CALL  CHICK  II«) 

WITHOUT!, 0110) 

WITHOUT!,  mOMHEAlER(I5>,I5*l, 45) 

WITHOUT!, fOOO)  (HEA!ER(I5)  ,15*14,24) 
NRlTH0UT!,f00l><HIAIBRd5>,I5«i,13> 

WITH  OUT! , T004 ) HEARER  ( 55 ) , HE A!ER  ( 54 ) 

ACCBPT‘II  THIS  TNE  PILE  YOU'RE  AFTER?  (1*YES/0«N0> :  ",FSPFN 
IP(FSPPN«EQ.O)COTQ  50 

190  ACCEPT"INTER  THE  TINE  SLICE  NUMBER:  ",OBSNBR 

IK 0!SM!R . 6B . HEARER 1 55  > . AN! , 0!SH!R . LE. HEADER (54) ) GOTO  200 
TYPI'ERRORt  THIS  HURIER  RUST  LIE  IH  THE  RANGE" 

TYPE" C ‘ , HEARER ( 55 ) , V , HEARER ( 54 ) , 1 "  1 " 

SOTO  1 90 
200  COHTIHUB 

240  F0NNBR«HRA!SR(37> 

I F ( PONHBR •  LE  < HHRP ) SOTO  240 

TYPE"ERROR«  YOU  REOUESTE!  HORE  THAN  512  SPECTRAL  COMPONENTS. 
SOTO  200 
240  CONTINUE 

Cmmm  DISI  ILOCE  AN!  ARRAY  LOCATION  OF  1ST  DATA  POINT  IS 
IA*0!SNBR-HEA!ER(53)*1 
I!«HEADER(57)i2 
CALL  6FDI(IA,I!) 

FRSTDBLI*  WIA 
FRST1P* ( IB-1 )/2il 

CALL  RDBLK(4,FRSTDBLK,DIST,4,ICNT,IER) 

IFdER.HE.TIGQTQ  300 

TYPE"READ  EOF!  SUCCESSFULLY  TRANSFERRED", ICHT,“  QTR  BUS" 
TYPE-PROCEEDING  U1TH  -,1CHT,“QTR  BLOCKS  TRANSFERRED." 

IER*1 

300  CALL  CHECK(IER) 

CALL  CLOSE ( 4 , 1 ER > 

CALL  CHECK  IER) 

I F ( OUT • EQ  *12) OUT! *  1 2 

C*»i«M  HOVE  THE  DATA  UP-FROHT  AND  IHTO  THE  ARRAY  TO  PLOT 
ACCEPT-ENTER  1  FOR  DB  PLOT/  0  FOR  HAGNITUDE:  ",IDP 
IFUDP.EQ.DCOTO  310 
DO  305  J1*1,HEADER<57> 

305  DATA( J1 )  *DIST(  JUFRSTDP  - 1 ) 

GOTO  320 

310  DO  315  J1*1,HEADER(57) 

315  DATA(Jl)«20iAL0G10(DISI(JUFRSTDP-i)> 

320  1F(0UT.EQ.0)G0T0  321 

URITE(OUTD, 0100)  (BATAdS), 15*1, FOHHBR) 

321  ACCEPT-OI  TO  CONTINUE?  (ENTER  ANY  NUMBER)  ",I5 

322  1F(15»E0.0)C0T0  324 

324  CONTINUE 

CALL  SRPH2<”  SPECTRAL  PLOT", 1 , DATA, U, PONHBR, 0,XRIN,XNAX,0> 
ACCEPT-OI  TO  CONTINUE?  (ENTER  ANY  NUHBER)  “ ,15 
I F ( 1 5  *  EQ  *  0 ) COTO  350 
CONTINUE 


350 


mo  miwmsii.o 

till  FltMT(1S#*NIKl< I  THE  FIRST  45  ELEMENTS  OF  THE  NEARER*-) 
fOOO  FRRHATCT3, “SPEECH  FILE  NAM  18*  *,1382) 

tOll  FORRAT ( T5 ,“ SPECTRAL  FILE  NARE  HAS*  *,1382) 

8004  FORRAT (T2, “FIRST  8FICTRAL  TSa  *  “,I5,“  AHO  LAST  »  “,I4> 

RETURN 

EN» 
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PILE:  PLTN 

LANGUAGE:  FORTRAN  5 

DATE:  September  21,  1982 

AUTHOR:  D.  Martin 

SUBJECT:  Graphics 

CALLING  SEQUENCE:  PLTN 

DATE  OP  LAST  REVISION:  September  21,  1982 

PURPOSE: 

This  routine  plots  distances  from  a  disk  file  computed 
by  DSTN.  The  plot  is  displayed  on  a  Tektronix  4010-1 
Graphics  Terminal  using  the  routine  GRPH2  by  G.  Shaw  as 
modified  by  L.  Kizer  and  D.  Z am bon. 

DESCRIPTION: 

Location:  DP4 : BRATCHET 

Size:  PLTN. FR 

PLTN.RB 

PROGRAM  USE: 

This  program  is  called  by  PLTO  and  calls  GFUB  and  GRPH2. 
For  a  specified  observation,  this  routine  plots  the  observa¬ 
tion  energy  as  the  first  data  point  and  the  distances  at 
locations  corresponding  to  their  phonet  numbers.  In  this 
way,  a  sort  of  scatter  plot  is  obtained. 

The  operator  is  first  prompted  for  an  option  to  print 
numerical  values  of  the  data  to  be  plotted.  The  values 
plotted  can  be  displayed  on  the  operator's  terminal,  written 
to  the  line  printer,  or  not  displayed  at  all.  Next,  the 


4132  bytes 
5970  bytes 
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operator  ia  prompted  for  the  name  of  the  disk  file  holding 
the  data  to  be  plotted.  The  file  header  is  read  and  identi 
fying  information  displayed.  Then,  the  observation  number 
is  requested.  The  program  will  display  the  plot  when 
signalled  by  the  operator  to  proceed. 

LIST  OF  VARIABLES: 


VARIABLE 

TYPE 

PURPOSE 

OUT 

Integer 

Switch  for  data  output 
option. 

OUTD 

Integer 

Channel  number. 

FILE3 

Integer  Array 

Holds  disk  file  name. 

OBSNBR 

Integer 

Observation  number. 

PONNBR 

Integer 

Number  of  data  points 
(one  plus  number  of 
phonets ) . 

FRSTDBLK 

Integer 

First  disk  block  to  read 

FRSTDP 

Integer 

First  data  point  in  disk 
block. 

DATA 

Real  Array 

Holds  data  to  be  plotted 

DIST 

Integer  Array 

Data  read  from  disk. 

SWITCH  SETTINGS: 

SWITCH  SETTING 

OUT  =  10  for  screen  output  of  plotted  data 

=  12  for  line  printer  output 

=  0  for  no  output  of  numbers 

RELATED  PROGRAMS; 

PLTO,  PLTS,  PLTA,  PLTT,  DRVR,  DRSQ,  S64,  S128,  DSTA, 
DSTN,  CHOOS,  CHKJ,  CHKO,  FCTR,  GFDB,  GRPH2. 
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C»»«««»  ROUTINE  PLTN.  THIS  ROUTINE  IS  FOR  FORTRAN  5  !  ! 

C  THIS  ROUTINE  IS  CALLED  BY  PLTO  AND  CALLS  THE  PLOT  ROUTINES 

C  GRPH2  AND  CFBD.  CRPH2  HAS  HRITTEH  BY  C.SHAH  AND  HQDIFIED  BY 

C  L.IIZER  AND  D.ZAHBON.  IT  IS  IN  THE  SPEECH  LAB  PR06RAH  LIBRARY, 

C  PLTN  PLOTS  DISTANCES  COMPUTED  BY  DSTN.  FOR  A  SPECIFIED 

C  OBSERVATION,  IT  HILL  PLOT  THE  OBSERVATION  EHERCY  AS  THE  FIRST 

C  DATA  POINT,  AND  THE  DISTANCES  AT  LOCATIONS  CORRESPONDING  TO 

C  THIER  PHONET  NUMBERS,  IN  THIS  HAY  A  SORT  OF  SCATTER  PLOT  OF 

C  N-BEST  DISTANCES  IS  PLOTTED. 

C  BY:  CAPT  DAN  MARTIN 

C  DATE:  9/21/B2 

C  SUBJ :  ACOUSTIC  ANALYSIS 

C  PLTN  ACCOMPLISHES  THE  FOLLOHING  TASKS: 


C 

1) 

GET  DISTANCE  FILE 

c 

2) 

DISPLAY  INFORMATION  FROM  HEADER  TO  IDENTIFY  THE  FILE 

c 

3) 

GET  OBSERVATION  NUMBER  AND  PHONET  NUMBERS 

c 

A) 

LOCATE  DATA  IN  DISK  FILE 

c 

5) 

PREPARE  DATA  FOR  PLOT  ROUTINE 

c 

6) 

CALL  THE  PLOT  ROUTINE 

C  FOR  MORE  INFORMATION  SEE  MY  THESIS  OR  THE  USERS  MANUAL. 

SUBROUTINE  PLTN 
OVERLAY  OPLTN 
INCLUDE" ARR A YP:F5 APS , FR" 

PARAMETER  HNDP«512 
PARAMETER  0TRBLK254 
REAL  DATA,XMIH,XRAX 

INTEGER  CSPCT,FILE3, HEADER, FSPFN,OUTD,OBSNBR,PHON1,PHONL 
INTECER  D1ST,F0NNBR,FRSTDP,FRSTDBLX,ICNT,D1,IER,D2,0UT 

COMMON  /  VALS  /  DATA ( HNDP ) ,DIST ( 1024 ) ,0BSNBR,FSPFN,0UTD,PH0N1 
COMMON  /  VALS  /  PHOHL , 1 CNT ,IER,FONNBR,FRSTDP ,FRSTDBLK,FILE3(13) 
COMMON  /  VALS  /  HEADER(QTRBLK) 

COMMON  /  VALT  /  D1I37) 

COMMON  /  VALU  /  D2I2054) , NCHO I CES , NBRPHQNE3 

TYPE"*»»*YQU  ARE  NOU  IN  ROUTINE  PLTD,FR»»*»” 

30  TYPE“DO  YOU  HANT  OUTPUT  OF  DATA  PLOTTED?” 

TYPE”ENTER:  10  FOR  SCREEN  OUTPUT  OF  NUMBERS- 

TYPE”  12  FOR  PRINTED  NUMERICAL  OUTPUT" 

ACCEPT"  0  FOR  NO  NUMBERS:  ",OUT 

IFIOUT.EQ.O.OR.OUT.EQ, iO.OR.QUT.EQ. 12)G0T0  AO 
TYPE-ERROR:  ILLEGAL  INPUT.” 

COTO  30 
AO  0UTD*10 

Ci*«h  GET  DISTANCE  FILE  NAME 

SO  ACCEPT-MENTER  NAME  OF  FILE  HOLDING  DISTANCES:  " 

READ (1 1,100 ) FILE3( 1) 

100  FORMAT I S 13 ) 

120  CALL  CLOSE ( 2 , IER) 

CALL  0PEN(2,F1LE3,2,1ER) 

CALL  CHECK  IER) 

Chhh  GET  NUMBERS  OF  1ST  I  LAST  TIME-SLICES  OF  SPECTRUM 
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160  CALL  RDBLK (2,0, HEADER, 1,IER) 

CALL  CHECK  C 1 ER ) 

WRITE(0UTD,9110) 

WRITE (OUTD, 9100) (HEADER (IS) ,15:1,65) 

WRITE (OUTD ,9000) (HEADER ( 15) ,15s 1,13) 

UR1TE ( OUTD ,9001 ) (HEADER (15) 

UR I TE ( OUTD ,9004) HEADER (40) /HEADER  <  4 1 ) 

URITEtOUTD ,9002) <HEADER( 15) ,15*27,39) 

WRITE (OUTD, 9005) HEADER (42) /HEADER (43) 

ACCEPT" IS  THIS  THE  FILE  YOU’RE  AFTER?  <1--YES/0*H0) :  “,FSPFN 
IF(FSPFH.EQ.O)COTO  50 

190  ACCEPT’ENTER  THE  OBSERVATION  NUMBER :  ",OBSNBR 

IF (OBSHBR.CE. HEADER (40) . AND . OBSNBR »LE . HEADER ( 4 1 ! ) GOTO  200 
TYPE"ERRORt  THIS  NUMBER  RUST  LIE  IN  THE  RANCE" 

TYPE** C ** /HEADER ( 40 )  ,**,** /HEADER i 4 1 ) 

GOTO  190 
200  CONTINUE 

240  FONNBR* HEADER (43) -HEADER (42) +2 

IF(FONNBR.LE.HNDP)COTO  260 

TYPE"ERROR:  YOU  REQUESTED  MORE  THAN  512  DISTANCES." 

COTO  1000 
260  CONTINUE 

CHHII  GET  DISK  BLOCK  AND  1ST  DATA  POINT  IS 
IA*QBSNBR-HEADER(40)*1 
IB=4*2«HEADER(47) 

CALL  GFDB(1A,IB) 

FRSTDBLK-'IAtl 

FRSTDP*1B 

CALL  RDBLX(2,FRSTDBLK,DIST,3,ICNT,IER) 

IF(IER.NE.9)CQT0  300 

TYPE"READ  EOF!  SUCCESSFULLY  TRANSFERRED", ICNT,"  QTR  BLKS" 
TYPE"PR0CEED1NC  WITH  ", ICNT, "QTR  BLOCKS  TRANSFERRED." 

I ER  *  1 

300  CALL  CHECK(IER) 

IF ( OUT. EQ. 12) OUTD =12 

C*m»m  ZERO-OUT  ARRAY  DATA 
DO  400  JJ*l,FONNBR+l 
400  DATA(JJ)*0.0 

FIRST  DATA  POINT  IS  OBSERVATION  ENERGY 
DATA( 1)*D1ST(FRSTDP»1) 

TYPE-OBSERVATION  NUMBER*", DISK FRSTDP) 

ACCEPT" IS  THIS  THE  ONE  YOU  WANT?  (1«YES/0*N0) :  ",FSPFN 
IF(FSPFN.EQ.O)COTO  120 
C»»*»»»  THEN  FILL  DATA  WITH  DISTANCES 

DO  410  JJ*FRSTDP+2,3+FRSTDP*2»HEADER(47>,2 
410  DATA(DIST( JJ)tl )*D1SK JJtl) 

IF (OUT »NE.O) WRITE (OUTD, 9100 )(DATA( 15) ,15*1, HEADER (43 )-KEADER(42)tl) 
310  ACCEPT-OI  TO  CONTINUE?  (ENTER  ANY  NUMBER)  ",15 

IF(I5.EQ.0)C0T0  312 
312  CONTINUE 

CALL  CRPH2("  DISTANCE  PLOT" ,1 ,DATA,U,FONNBR,0,XMIN,XHAX,0) 
ACCEPT"OI  TO  CONTINUE?  (ENTER  ANY  NUMBER)  ",15 
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350 


IF(I5.EB.0)CQT0  350 
CONTINUE 
CALL  RESET 
9100  FORMAT ( 5G 1 2 • 4  > 

9110  F0RHAT<T5, “HERE'S  THE  FIR6T  «5  ELEHEHTS  OF  THE  HEADER:") 

9000  F0RHAT(T3, “DISTANCE  FILE  NAME  IS:  ",13S2) 

9001  F0RHANT5, “OBSERVATION  FILE  NAHE  HAS:  ",13S2) 

9002  F0RHAT(T5,“PH0NET  FILE  NAHE  HAS:  ",13S2) 

9004  FQRHAT(T2, “FIRST  OBSERVATION  TS«  *  ",15,"  AND  LAST  *  " ,16) 

9005  F0RHAKT2, "FIRST  PHONET  TS«  *  ",I5,"  AND  LAST  *  ‘,14) 

1000  RETURN 

END 
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FILE: 

LANGUAGE: 

DATE: 

AUTHOR: 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OF  LAST  REVISION: 


PLTT 

FORTRAN  5 

September  21,  1982 
D.  Martin 
Graphics 
PLTT 

September  21,  1982 


PURPOSE: 

This  routine  plots  an  integer  disk  file  in  units  of 
from  one  to  512  integer  words  on  a  Tektronix  4010-1 
Graphics  Terminal  using  the  routine  GRPH2  by  G.  Shaw  as 
modified  by  L.  Kizer  and  D.  Zambon. 

DESCRIPTION: 

Location:  DP4 : BRATCHET 

Size:  PLTT. FR  2840  bytes 

PLTT. RB  3290  bytes 


PROGRAM  USE: 

This  program  is  called  by  PLTO  and  calls  GFDB  and 
GRPH2.  It  will  plot  segments  of  an  integer  disk  file: 
the  segment  length  can  be  from  one  to  512  integer  words. 

Hie  routine  was  designed  to  be  used  to  plot  segments  of 
speech  files  and  to  plot  distances  between  several  con¬ 
secutive  observations  and  phonets.  The  operator  is  prompted 
for  the  disk  file  name,  the  segment  number,  and  the  segment 
size.  To  illustrate,  say  the  header  of  a  file  computed  by 


either  DRSQ,  DSTA,  or  DSTN  is  to  be  examined.  Using  this 


routine,  the  operator  could  specify  the  first  segment  of 
256  words  along  with  a  print  option  to  examine  the  header. 
The  plot  could  be  discarded.  On  the  other  hand,  suppose  a 
curious  segment  of  a  speech  file  is  to  be  examined.  One 
determines  which  disk  block  contains  the  curious  segment, 
and  uses  PLTT  to  plot  that  disk  block. 

LIST  OF  VARIABLES: 


VARIABLE 

TYPE 

PURPOSE 

OUT 

Integer 

Switch  for  data  output 
option. 

OUTD 

Integer 

Unit  number. 

FILE3 

Integer  Array 

Holds  disk  file  name. 

OBSNBR 

Integer 

Disk  file  segment. 

PHONL 

Integer 

Disk  file  segment  length 

FRSTDBLK 

Integer 

First  disk  block  to  read 

FRSTDP 

Integer 

First  data  point  in  disk 
block. 

DATA 

Real  Array 

Holds  data  to  be  plotted 

DIST 

Integer  Array 

Data  read  from  disk. 

SWITCH  SETTINGS: 

SWITCH  SETTING 

OUT  =  10  for  screen  output  of  numbers  to  be  plotted 

=  12  for  line  printer  output 

=  0  for  no  numerical  output 

RELATED  PROGRAMS: 

PLTO,  PLTA,  PLTN,  PLTS,  DRVR,  DRSQ,  S128,  S64,  DSTA, 
DSTN,  CHOOS,  CHKJ,  CHKO,  FCTR,  GFDB,  GRPH2. 
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ROUTINE  PLTT.  THIS  ROUTINE  IS  FOR  FORTRAN  511 
THIS  ROUTINE  IS  CALLED  BY  PLTO  AMD  CALLS  THE  PLOT  ROUTINES  6RPH2 
AND  CFBD.  SRPH2  NAS  NRITTEN  BY  C.SKAU  ANB  MODIFIED  BY  L.KIZER 
AND  D.ZANBON.  IT  IS  IN  THE  SPEECH  LAB  PROCRAH  LIBRARY. 

THIS  ROUTINE  PLOTS  INTE6ER  DISK  FILES.  IT  NUL  PLOT  SE6HENTS 
OF  THE  FILE}  EACH  SB6HENT  IS  SPECIFIED  BY  ITS  NUMBER  AMD  POINT 
LENCTH.  THE  POINT  LENCTH  CAN  BE  FROM  1  TO  512  INTECER  HQRBS. 

THE  USER  IS  PROMPTED  FOR  THIS  INFORMATION  AS  HELL  AS  FOR  THE 
FILE  NAME. 

BY :  CAPT  DAN  MARTIN 
DATEi  9/21/62 
SUBJ:  ACOUSTIC  ANALYSIS 
PLTA  ACCOMPLISHES  THE  FOLLONINC  TASKS: 

1)  GET  FILE 

3)  GET  UNIT  NUMBER  AND  UNIT  SIZE 

4)  LOCATE  DATA  IN  DISK  FILE 

5)  PREPARE  DATA  FOR  PLOT  ROUTINE  GRPH2 
4)  CALL  THE  PLOT  ROUTINE 

FOR  MORE  INFORMATION  SEE  THE  USERS  MANUAL  OR  MY  THESIS. 

SUBROUTINE  PLTT 
OVERLAY  OPLTT 
INCLUDE* ARRAYP:F5 APS. FR* 

PARAMETER  HNDP*512 
PARAMETER  QTRBLK*254 
REAL  DATA,XHIN,XHAX 

INTEGER  GSPCT,FILE3, HEADER, FSPFN,OUTD,OBSNBR,PHONl,PHONL 
INTEGER  D1ST,F0NNBR ,FRSTDP ,FRSTDBLK ,ICHT ,Dl , 1ER ,QUT 

COMMON  /  VALS  /  DATA(HNDP),DIST(1024),OBSNBR,FSPFN,OUTD,PHON1 
COMMON  /  VALS  /  PHOHL , ICNT , IER ,FONHBR , FRSTDP , FRSTDBLK , F I LE3 (13) 
COMMON  /  VALS  /  HEADER(QTRBLK) 

COMMON  /  VALT  /  Dl(37) 

TYPE“»»»«YOU  ARE  NON  IN  ROUTINE  PLTT.FR»»«*“ 

30  TYPE*DO  YOU  WANT  OUTPUT  OF  DATA  PLOTTED?* 

TYPE*ENTER:  10  FOR  SCREEN  OUPUT  OF  NUMBERS" 

TYPE"  12  FOR  PRINTED  NUMERICAL  OUTPUT" 

ACCEPT*  0  FOR  NO  NUMERICAL  OUTPUT:  ",OUT 

1F(OUT.EO.O.OR.OUT.EQ. 10. OR » OUT «EQ. 12) GOTO  40 
TYPE"ERROR:  ILLEGAL  INPUT." 

GOTO  30 
40  OUTD’IO 

60  CONTINUE 

C«m«m  GET  INTEGER  FILE  NAME 

ACCEPT"»ENTER  NAME  OF  FILE  HOLDING  INTECER  DATA:  " 

READ (11 ,100)F1LE3(1) 

100  FORMAT ( S13 ) 

120  CALL  0PEN(5, FILES, 1,IER) 

CALL  CHECK(IER) 

190  ACCEPT“ENTER  THE  TIME  SLICE  NUMBER:  ",OBSHBR 

220  ACCEPT"ENTER  THE  NUMBER  OF  POINTS  IN  TIME  SLICE:  ",PHONL 

IF(PHONL.LE.HNDP)COTO  260 
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TYPE** ERRORt  YOU  REQUESTED  NORE  THAN  512  POINTS." 

GOTO  220 
200  CONTINUE 

Oiiiii  USE  ILOCK  THAT  HOLDS  1ST  DATA  PQIHT  IS 
1A-0BSNBR 
IS*  PHONL 

CALL  CPDI(IA|1B) 

FRSTDILK* I A 
FRSTDP>IB 

CALL  RDBLX(S;FRSTDDLXjDIST;4f ICNT,IER) 

XF < IER • HE • 9 )COTO  300 

TYPE"READ  EOF!  SUCCESSFULLY  TRAHSFERRED" ,ICHT,“  OTR  BUS 
TYPE** PROCEEDING  WITH  ",ICNT,“QTR  SLOCKS  TRANSFERRED." 

I ER  *  1 

300  CALL  CHECK (IER) 

CALL  CLOSE(S)IER) 

CALL  CHECK (IER) 

IF(OUT.EO.O)COTO  310 
IF(0UT.EQ.12)0UTD*12 

Ciiiiit  HOVE  THE  DATA  UP-FROHT  AND  INTO  THE  ARRAY  TO  PLOT 
310  lU’PHONL 

DO  320  J1*1,IU 

320  DATA(  Jl)*DISTt  JWFRSTDP-I) 

IF(0UT.EQ.0)C0T0  322 
TYPE“OUTD*“ >OUTD 

URITE(OUTD, 9100) (DATA( 15), 15*1, IU) 

322  ACCEPT-OX  TO  CONTINUE?  (ENTER  ANY  NUHDER)  ",I5 

IF(I5.NE.0)C0T0  324 
324  CONTINUE 

CALL  CRPH2("  TINE  SLICE  PLOT ~ , 1 , DAT A , U ,  1 U , 0 , XH IN , XH AX > 0 > 
ACCEPT'OX  TO  CONTINUE?  (ENTER  ANY  NUNBER)  ",I5 
IF( I5.EQ.0)C0T0  350 
350  CONTINUE 

UOO  F0RNAT(5C12.4) 


RETURN 

END 


PILE; 

LANGUAGE: 

DATE: 

AUTHOR; 

SUBJECT: 

CALLING  SEQUENCE: 

DATE  OP  LAST  REVISION 


CHKO 

FORTRAN  5 

September  21,  1982 
D.  Martin 
Acoustic  Analysis 
CHKO  (IER) 
September  21,  1982 


PURPOSE: 

This  routine  is  called  by  DSTA  and  DSTN,  and  calls  no 
subroutine.  It  checks  the  value  of  IER  returned  from  an 
OPEN  file  attempt.  It  returns  a  caution  statement  if  the 
unit  number  is  in  use. 


DESCRIPTION: 

Location:  DP4:BRATCHET 

Size:  CHKO. FR 

CHKO. RB 

LIST  OF  VARIABLES:  None. 

ARGUMENT  STRUCTURE: 

ARGUMENT  TYPE 

IER  Integer 


586  bytes 
288  bytes 


PURPOSE 

Holds  error  code  from  OPEN 
system  call. 


PROGRAM  USE: 

Tto  be  called  by  DSTA  and  DSTN  to  provide  caution 
statement. 


223 


RELATED  PROGRAMS; 

PLTO,  PLTA,  PLTS,  PLTT,  PLTN,  FCTR,  GRPH2 ,  GRDB,  DRVR, 
DRSQ,  DSTA.  DSTN,  CHKJ,  CHOOS,  S128,  S64. 
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ROUTINE  CHKO.  THIS  ROUTINE  IS  FOR  FORTRAN  3  !  ! 

THIS  ROUTINE  IS  CALLED  BY  DSTA  AHD  DSTN,  AND  CALLS  NO  SUBROUTINES. 
IT  CHECKS  THE  VALUE  OF  IER  RETURNED  FROH  AN  OPEN  FILE  ATTEHPT. 

IT  RETURNS  A  CAUTION  STATEMENT  TO  THE  SCREEN  IF  THE  UNIT  HUHBER 
IS  IN  USE. 

BYt  CAPT  DAN  MARTIN 

DATE:  9/21/82 

SUB J  t  ACOUSTIC  ANALYSIS 

FOR  MORE  INFORMATION  SEE  THE  USERS  MANUAL  FOR  DSTA,  DSTN,  OR  THIS 
ROUTINE.  OR  SEE  MY  THESIS. 

SUBROUTINE  CHKO(IER) 

OVERLAY  OCHKO 
INTEGER  IER 

I F ( I ER . NE . 3098 ) GOTO  100 
IER*1 

TYPE” •••••■CAUTION:  UNIT  NUMBER  IS  IN  USE!»«»«»»* 

RETURN 

END 
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